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Two dljncntlonal tubionic wind tunnol t«itt navt boon conducud on a 
20% thlcknaaai chord ratio circulation controlled alliptic aerofoil 
aection equipped with forward and reverae blowing alots. Overall 
performance measurementa were made over a range of trailing edge 
blowing roomentun coefficients from 0 to 0.04; some included the effect 
of leading edge blowing. The effective incidence was determined 
experimentally and lift augmentations, 3C./3C , of 70 were obtained at 

li VJ 

low blowing rates. 

A detailed Investigation of the trailing edge wall jet, using split 
film probes, hot wire probes and total head tubes, provided inoasure- 
ments of mean velocity components, Reynolds normal and shear stresses, 
and radial static pressure. Corrections for the effects of ambient 
temperature variation, flow angle and shear flow gradient upon the 
various probes were examined and some corrections for the low band- 
width of the split film probes proposed. 

In some cases, the effects of slot height and slot lip thickness were 
investigated. The results were mostly taken at a geometric incidence 
of 0°. 


The closure of the two dimensional angular momentum and continuity 
equations was examined using the measured data, with and without 
CO section, and the difficulty of obtaining a satisfactory solution 
illustrated. 

The experimental results have led to some suggestions regarding the 
nature of the flow field which should aid the understanding of Coanda 
effect and the theoretical solution of highly curved wall jet flows. 
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1. INTRODUCTION 
If 1 Th« Coanda Eff»ot 

Th» affact by which a fluid jat attachaa itaalf to an adjacant aurfaca 
and ramalna attachad waa Initially obaarvad by M. Hanri Coanda, after 
vftom the affect waa namad. 'Coanda affect* ia capable not only of 
attaching a free jet to a aurfaca but can alao enable a tangential jet 
to negotiate and remain attached to a highly curved walli The effect 
producea very strong entrainment of the aurroundlng fluid, independ- 
ently of whether the external fluid ii moving or stationary, and 
significantly reduces the surface static pressure under the jet. The 
point at which the flow separates from a curved surface in a two 
c.lmensional case can be controlled by the jet blowing momentum. The 
detailed physics of the effect are still not wholly understood. 

M. Coanda investigated many applications of the effect including 
leading edge blown aerofoils and thrust augmentors. Reference 1 is an 
example of this work, reproduced after the war. Further app.icatlons 
have been investigated by other researchers, for example, the blown 
cylinder, upper surface wing/flap blowing (as on the Buccaneer and 
YC-14 aircraft) , fluidic amplifiers and wall blowing to improve wind 
tunnel diffuser operation (see Figure 1). More recently interest has 
been focused upon circulation control aerofoils and their potential 
benefits for helicopter rotors. 


If aerofoils are considered, then a conventional sharptrailing edged 
aerofoil exhibits the well known 'Kutta* condition. This states that 
the rear stagnation streamline must emerge from the trailing edge in 
order to avoid discontinuities in velocity and pressure. Hence the 
circulation around the aerofoil is uniquely defined by the aerofoil 
geometry, incidence and free stream velocity, if however# the trailing 
edge of the aerofoil ia rounded, then the rear stagnation point is 
free to move, depending upon the other parameters. If the aerofoil is 
an ellipse aligned at aero incidence to the free stream direction, then the 
upper and lower surface separation points should be located at the same 
chordwise station and the nett circulation will be sexo. If now a jet 
of fluid is injected tangentially into the upper surface boundary 
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layor, near the trailing edge, the Coanda effect will entrain the 

boundary layer and delay the eeparatlon of the upper aurface flow. 

Thle cauiea a nett increaee in the circulation around the aerofoil. 

The momenturo of the blowing Jet now controla the position of the tear 

stagnation pointj the aerofoil la subject to 'Circulation Control' by 

blowing. If the blowing jet is si'.rong enough to discharge excess 

momentum into the wake/ then the aerofoil performs in a similar manner 

to one fitted with a jet flap, The lift is no longer produced solely 

by delaying the upper surface separation but has a jet reaction 

thrust component which reduces the nett lift augmentations 3C. /aCp, 

L 

( see Figure 2 ) . 


1 • 2 The Applications of Circulation Control 

The application of circulation control-by-blowing to helicopter rotors 
and more recently to stopped rotor aircraft (X-wing, see Figure 3) is 
shown in References 2 and 3. 


Conventional helicopter rotor systems suffer from a variety of 
problems, including vibration, mechanical complexity, high drag and 
retreating blade stall at high forward speed. It has been shown 
(References 4, 5 and 6) that a circulation control rotor is capable of 
reducing many of these problems. 

Circulation control aerofoils are able to develop lift independently 
of incidence, and to a first approximation, at low blowing rates, 
velocity. Hence the rotor may develop lift purely as a function of 
blowing rate. At low advance ratios - w ■ helicopter forward speed/ 
rotor tip speed less than 0.5 - the lift may be increased at the fore 
and aft asinuth positions of the rotor disc to produce a more 
efficient value of the thrust coefficient/solidity ratio (C^/o) . In 
particular, the hover performance may be greatly inproved since a 
circulation control rotor is capable of producing high lift 
coefficients compared with conventional rotors. At advance ratios in 
excess of 0.7 the area of reversed flow on the retreating blade side 
tends to increase (see Figure 4) . It has been shown by Ottensosar^ 
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that incorporation of a aacond blowing alot on tha loading adga of tha 
blada. ajacting air tangantial to tha aurfaca and in an oppoaita 
aanaa to Uia normal fraa atraam diraction, anablaa tha blada to 
produce lift in a ravaraad flow region. Tha dual alota may ba blown 
independently aa required or continuoualy with little aignificant loaa 
of lift. Typical flowa around a dual slotted aerofoil are ahown in 
Figure 5. Uaeful lift has been shown on test rigs at advance ratios 
greater than 2 (approximately 400 knots forward speed) . 


Three distinct flight regimes of Improved VTOL performance arise as a 

consequence of circulation control: 

i) High forward speed (advance ratios of the order of 2) by 
incorporating blowing slots at the leading and trailing edges. 

ii) High forward speed at highei advance ratios permitted by slowing 
the rotor and incorporating leading edge blowing only. 

iii) Low speed performance as a conventional helicopter and trans- 
lation to a 'stopped rotor' for high speed flight (X-wing mode). 
Transition would be between 100 - 180 Kts. In this case dual 
blowing would be available for the helicopter mode and 'trailing 
edge' only blowing would be used to augment the lift whilst 
operating as a stopped rotor aircraft. 


All of these techniques would use the blowing system as an attitude 
control system offering improved control response at low forward speed. 


Each of these systems offers many advantages over the current 
technology improvements in helicopter performance being made available 
by tip geometry, blade construction and airframa improvements. 


Inherent in circulation control aerofoils is the increased blade 
stiffness which will enable the rotor to be smaller and lighter. A 
circulation control blade may be as much as lo times stiffer in 
bending c<»pared with a conventional 'o' spar blade. 
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Sines collsctlve and cyclic blowing can ba aubatitutad for collactiva 
and cyclic pitch, tha complaxity, aiaa and waight of a circulation 
control hub can ba significantly laaa than ita convantional countar- 



part. 

This offers three main advantages; 

j 

i) 

reduced hub aerodynamic drag 

I 

. f 

11) 

reduced hub weight 

f ' 

iii) 

reduced rotor system maintenance 


Tha third advantage being moat significant in terms of operating costs 
and time on the ground. A circulation control rotor also has the 
advantage of having fewer dynamic force compaients and hence a higher 
harmonic control system is possible giving leas vibration and improved 
gust response. 

Reference 8 describes the first practical application of circulation 
control to a helicopter rotor (Kaman XH-2/CCR) and indicates the 
simplicity and effectiveness of installing a circulation control rotor 
on an existing airframe. The Kaman flight demonstrator is essentially 
a low speed test bed; a more realistic stopped rotor has been tested 
in the NASA-Ames 40 ft x 60 ft wind tunnel. 

1.3 Previous Experimental Work 

With tha advent of these advanced technology applications, the existing 
knowledge of the Coanda effect and circulation control aerofoils has 
had to be reviewed. The basic aim of tha research was, and still is, 
to produce a reliable prediction method for estimating the performance 
of practical applications of the Coanda effect, in particular, 
circulation control aerofoils. 

Accepting that the heart of the problem is to solve the highly curved 
wall jet flow, experimental research over the last two decades can be 
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grouped In three categories i 

I) plane jets 

II) curved jets 

ill) circulation control aerofoils 


An excellent recent summary of wall jet experiments is given in 
Reference 9. 

1*3.1 Plane jet experiments 

References 10 - 18 are given as examples of experimental work available 
on plane wall jets. These reports cover a variety of effects including 
longitudinal pressure gradient, external flow, jet turbulence level 
and slot geometry. 

The self preserving, simple plane wall jet flow has been shown to 
exhibit a linear half velocity growth rate given by 


^ V2 

dx 


0.073 ± 0.002 


( 1 . 1 ) 


however, this region may not begin until some 50 slot widths down- 
stream. The growth rate is some 30% leas than the corresponding free 
jet value. Paw useful measurements exist within the important starting 
region of the flow, although the study of lrwin^° is particularly 
thorough in the self preserving region. Also, good agreement with 
genaral mosientum balances is somewhat lacking and it has been 
suggested in Reference 9 that this is due mainly to unsatisfactory 
turbulence measurements. The problem of matching results from 
different test apparatus will be further illustrated in the later 
sections of this report. 
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S«v«ral axamples for the caae of a plana turbulant wall jat in tha 
praaanca of an axtarnal atraani hava baan givan. Thaaa hava a mora 
diraot ralationahip to tha blown boundary layar application, but tha 
addition of a furthar variabla, tha axtarnal flow, raaraly appaara to 
coraplicata tha test situation. The work of Kruka and Eskinaai^^ la 
considered to be one of the most reliable in this field, indicated 
in this work are the variation of maximum velocity decay, jet and 
free shear layer growth. Additionally the difficulty of providing a 
single similarity profile which applies across the whole flow is 
demonstrated. Gartshore and Newman^® provide some interesting 
comparisons between various test results and a proposed theory. The 
theory was used to investigate some assumptions generally made 
regarding the shear stress profile across a plane wall jet. This 
work emphasises how any prediction method has to be based on well 
documented, reliable teat results. 

1*3.2 Curved jet experiments 

References 19 - 25 are a selection of works concerned with wall jets 
around curved surfaces. The radial pressure gradient and extra strain 
on the fluid due to curvature make a convex curved flow more 
susceptible to spanwise irregularities and hence many of the experi- 
mental results must be viewed with suspicion. True two-dimensional 
flow around a cylinder is considerably more difficult to obtain than 
with a plane jet. The interaction effects at the extremities of the 
blowing jet may cause large longitudinal vortices to be shed into the 
jet flow. The slot geometry and, in particular, small spanwise 
irregularities in the slot lip, may also cause severe three- 
dimensional effects. The growth rate of the curved jet is considerably 
higher than that of a plane jet due to the extra strain imposed by the 
curvature upon the flow. A further problem involved with the testing 
of curved, rather than plane, wall jets is the provision of 
satisfactory instrunentation for measurements in the highly curved 
flow. In general, to avoid excessive blowing requirements, a small 
cylinder and slot height is chosen and this reduces the width of the 
flow and hence demands smaller spatial resolution of any probes 



posltiontd within th« flow. With the higher growth rate, th, flow 
«l.o ceaeee to be effectively parallel to the local lurface. Thie, 
coupled with the high turbulence levels .xperienced at the edge of 
the jet and the increased entrainment, make determination of any 
turbulence parameters exceedingly difficult. Reference 21, by Wilson 
and Goldstein, may be considered to give a good indication of the 
turbulence properties within a highly curved wall jet in still air. 

In an attempt to produce a self preserving curved wall jet some 
researchers have investigated the log spiral surface, in which the 
curvature is a function of the distance from the slot. 


R - K 
o 


( 1 . 2 ) 


where R ■ surface radius 
o 

S ■ distance from slot. 

Giles et al^^ is an example of this work. An approximately self 
preserving wall jet is created, although the growth rate is larger 
than for circular cylinders of comparable surface radius. Little 
reliable turbulence data is available for this type of flow. 

. 3 . 3 circulation cortrol aerofoil experiments 
There are two main areas of interest in the aerofoil experiment.! 

1 ) overall force and pitching moment coefficients 

ii) mean values and turbulence parameter, in the curved wall jet. 

Early work at N.G.T.E. and N.P.L. provided much information upon the 
high lift capabilities of cylinders and ellipses with circulation 
control by blowing. Ihe majority of the work however was performed at 


high blowing momantusi coefficients, C^, in excess of 0.5, 
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*nd hanct waa not of aijnlflcant uaa in the current appUcatlona. 

More recently a aeriea of atudlea have bean performed; References 1 , 

26 - 34 are good examples of this work. Typically, lift coefficients 
in excess of 6.5 for < 0.25, Lift augmentations 3Cj^/3c of greater 
than 70, and section nffective lift/drag ratios of more than 30 have 
been demonstrated. 

Only references 28, 32, 33 indicate any detailed trailing edge measure 
roenta. Englar used a hot film surface mounted shear stress meter to 
measure wall shear stress beneath the wall jet. Some indications of 
radial static pressure distributions were also obtained from static 
tappings located on a machined flange mounted normal to the trailing 
edge surface. The aerofoil was a modified 20% ellipse of 11 inch 
chord, tested up to high subsonic speeds. Kind^^ produced velocity 
profiles and longitudinal turbulence profiles at various angular 
positions around the trailing edge using a single hot wire probe. 

This aerofoil was again a 20% uncambered ellipse of inch chord 
and was tested at low subsonic speeds. Jones produced a large 
chord ( -V 4 ft) modified elliptic model in an attempt to increase the 
size of the measuring zone. Velocity and shear stress profiles were 
produced, again at low subsonic speeds, using a single slant hot wire 
probe. Examples of these results, where of interest, will be used in 
Chapter 5. 

The main reason for the lack of detailed investigations around the 
trailing edge of a circulation control aerofoil arises from the 
cesaplexity of the testing environmant. The small scale of the flow 
(slot heights of 0.5 mm are typical), the high curvature of the 
trailing edge and the problem of establishing two-dimensional flow 
around the aerofoil have restricted many researchers to simple overall 
pressure measurements. 


1 • 4 Prtviout Theoretical Work 


Th« th«or»tlc«l work can b* dlvidad into two areas i 

i) Ths wall jet without a co-flowing stream 

11) The wall jet with a co-flowlng stream/aerofoll trailing edge 
wall jet flow. 


The theoretical prediction of the first type of flow is greatly 
simplified since the boundary conditions are invariant from case to 
case. This suggests that both integral and differential methods 
should succeed since the velocity profile can be described simply, by 
a single equation. The bases of the two methods are summarised by 
Newman and Irwin and their relative qualities discussed. In 
general the technique used is to split the flow into a series of 
strips and to solve the assumed equations by establishing the boundary 
conditions at the edges of the strips. Unfortunately the extension of 
these types of technique to a trailing edge wall jet flow appears 
impractical. 

A great variety of wall jet velocity profiles may exist around the 
highly curved trailing edge of a circulation control aerofoil (see 
Figure 6) and this precludes the use of a simple universal set of 
equations to describe the velocity profiles. For this reason, the 
Integral technique will be only briefly discussed (see Section 1.4.1) 
and the reader is referred to the wor)c of Kind^^ for further details. 

In general, the solution of the flow around a circulation control 
aerofoil follows a similar pattern regardless of closure technique. 

i) determine the pressure distribution around the aerofoil using a 
potential flow solution. 

li) calculate the lower surface boundary layer development to 

separation using a suitable method, visually of the integral type. 
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ill) calculate the upper aurfaca flow davelopmant aa far aa the 
blowing alot, again uaually by an Integral method. 

Iv) cloaa the flow by calculating the curved wall jet development, 

from the atartlng condltlona of the upper aurface boundary layer 
at the alot and Cy, to separation. It la usually assumed that 
the wall jet separates at the same pressure as the lower surface 
boundary layer as suggested by Thwaitea^®. 

1.4.1 The Integral method 
37 

Dunham modified Spalding's plane wall jet method by Incorporating 
corrections for curvature and entrainment to represent a wall jet on a 
circular cylinder. Kind^^ then extended the method to an elliptic 
aerofoil; a block diagram of his method is given as Figure 7. However 
as previously mentioned, these Integral type calculations have limited 
applicability. 

1*4.2 The finite-difference method of Dvorak and Kind^^ 

38 

Dvorak and Kind used a similar overall calculation scheme which 
includes allowance for viscous effects, and uses a finite difference 
calculation to determine the wall jet development (see Figures 8, 9) . 
The finite difference scheme is initiated by merging the final 
boundary layer profile at the slot with a nearly uniform slot flow 
profile. The calculation mesh distribution is arremged with points 
concentrated in the regions of high shear. 

The static pressure distribution both normal and tangential to the 
surface is obtained from thn potential flow calculation. A correction 
to take account of the excess momentum flux in the jet (J ) is 
applied to the deduced radial static pressure distribution. A 
variation of «t the surface to sero at the velocity minimum 

is superimposed upon the distribution. 


The equations of motion are then closed using an eddy viscosity model 


for th« Reynold's ihtar itress tarn. Tha modal in Ita final form la 


I 

i' 



i 

! 


-u' V* ■ V 


|li - c 

t ay ''1 1+Hy 
eddy viacoaity 


f(«) . i 


(1.4 


Ihia ia largely baaed on aimpla mixing length argumenta. Valuaa for 
the additional empirical curvature correction, have been quoted 
from 1 to 25 by varioua reaearchera. In thia model, the value of 
has largely been tailored to fit the small amount of exiating 
experimental data (Kind^^) and varies thus; 


For y 


For 


m 


- i82n - 32n^ + i9on^ 

- JL 


■m 


> y. 


m 


1 

n - 


- 33 - 32n 


‘y - 

^^max “ 


(1.5) 


This type of variation is assumed since, at some small distance from 
the wall (y < y^) , where 3u/3y is still large and positive, must 
also be large to allow u*v* to become positive, as has been shown by 
Jones and Wilson and Goldstein^^. Thia is a possible indication 
that the present eddy viscosity model ia inappropriate for this type 
of flow. 


Ihe finite difference scheme *diich devel^s the wall jet flow to 
separation is checked against the previously stated Thwaites equal 
pressure criterion, and the aerofoil calculations are then iterated to 
a satisfactory closure. 
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Tht diBcr»t« vorfx wthod of Swlth^^ 

R«c*ntly, how«v«r, « n«w ttchniqu* hat baan davalopad by at 

al to datarmlna tha davalopmant of tha curvad wall Jat on an aaro- 
foil. At praaant, tha mathod uaaa a almpla potantial flow aolutlon 
around a cylindar with tha wall jat rapraaentad by a aariaa of 
diacrata point vortlcaa (and thair Imagaa) ahad from tha alot lip at 
diacrata time intarvala. Tha vorticaa rotata in tha sanae which 
deacribaa tha outar free ahaar layer of the wall jat. Tha davalopmant 
of the wall jet ia then calculated uaing the mutual influence of the 
vorticaa to predict ateady state separation. Tha smmnation of the 
induced effects of the vorticaa upon each other is calculated and tha 
locus of the vortex stream is predicted (see Figures 10, 11). Wie 
initial strength of the vortices is based upon the velocity difference 
at the slot, representing C , and the vortux strength then decays 
exponentially representing the dissipation of the vortex energy. An 

artificial viscosity has also been found necessary to stabilise the 
vortex motion. 

This method, simple in technique, has shown reasonable agreement with 
experiment, as shown in Figure 12, but is very dependent upon the 
prescribed decay rate of the vortex strength and the vortex shedding 
frequency (i.e. the time interval between each successive vortex). At 
present, the technique is undergoing development in order to remove 
the empirical decay rate and improve prediction. Further discussion 
of this technique and its relevance to the current study will be given 
in Chapters 5 and 6. 

1*5 The Present Investigation 

As has baan stated, there exists a requirement to be able to predict 
tha performance of Ooanda flows, in particular for circulation control 
aerofoils of various geometries. The existing theories lack universal 
applicability and so it was decided to undertake a detailed 
experimental study in order to provide more data and a better under- 
standing of the mechanics of a Coanda flow. 
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Th« tait «*rofoll wai to ba mounted in th* 7 ft k 5 ft high tpaad 

working Motion of the Bath Univaraity wind tunnal (aaa Figure 13 and 

Baction 2.2). The blowing aupply waa to ba delivered from a 

compraaaor capable of aupplying 1300 c.f.m. (r.A.D.) at a nominal 

80 p.B.l.g. 

The following design guidellnei were adopted in the light o. prevloua 

work and the available facilitiee. 

i) The chosen aerofoil section should be an ellipse as this is a 
section which is easily transformed mathematically, simplifying 
the evaluation of the theoretical pressure distribution around 
the aerofoil. 

ii) The thickness : chord ratio of the ellipse was chosen as 20% to 
enable direct comparison with previous research, particularly 
the work of Kind^^. 

iil) The chord should be as large as possible within the restraints 
of the working section in order to enlarge the scale of the 
boundary layer at the trailing edge. 

iv) It should be possible to adjust the incidence over the range 
±10 from the tunnel centreline. 

V) The aerofoil should Incorporate a plenum chamber and the skin of 
the aerofoil should not deflect significantly when blowing 
pressure was applied. 

Vi) Slots should be Incorporated at the leading and trailing edgea 
in order to examine the dual /reverse blowing cases. The slots 
should exhaust as near tangentially to the local surface as 
poasibla and should cause a minimum of discontinuity of the 
surface profile. 
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vii) The slot height should bs adjustsbls within ths optlnun range 
of Blot height i chord ratio as auggeatad by previous work 
(0.001 < h/c 5 0.002 Bnglar and Willlama^®) . 

vlll) The flow should be as near to ths two-dimensional case as 

poBsibla« either by using a large aspect ratio or by Incorp- 
orating a secondary blowing eye tom. 

lx) A travoralng system wae to be mounted on, or in, the model, to 
enable hot wire and pressure probes to be accurately positioned 
around the trailing edge of the aerofoil. 

It was proposed to examine the performance of this circulation 
control aerofoil within the range of blowing momentum coefficients of 
0 to 0.1, this being the most efficient regime in terms of maximum 
lift augmentation. 

1.6 A Brief Guide to the Contents of the Remainder of this Thesis 

The design considerations raised in Section 1.5 will be discussed in 
more detail in Chapter 2. Also in Chapter 2, a full explanation of 
the precautions taken to control the two-dimensionality of the flow, a 
description of the probe positioning devices, the wind tunnel and the 
blowing air supply will be given. 

Chapter 3 contains all the information on the instrumentation used in 
this experiment. The calibration, operation and sources of errors of 
both wire and split film anemometer probes are considered. 

Chapter 4 is a presentation of the results obtained both for the 
overall performance and the detailed trailing edge investigation. 

In Chapter 5 soma suggestions are made regarding the mechanics of a 
Coanda flow field in the light of the experimental results. The 
ia^llcatlons of the suggestions and their effect upon the inter- 
pretation of the results is discussed , Justification for the 


li 


■uggcttsd flow fl«ld ii given and the llnltationa of the 
Inatrumantation are llluitratad. 


Chapter 6 review^ the incorporation of the euggeeted flow field into 

a full aerofoil perfomnee leheme and aleo dieouaaee the deficienciee 
of current nethode. 

Where poaaible, the parametera ueed have been made non-dlmenelonal. 
Some reference! to both the Imperial and Metric systema hae been 
unavoidable; herever the unite are clearly stated in all caiee and 
the author apologieea for any inconvenience thia may cause. 


2 , APPARATUS 

2.1 Aerofoil D>»iqn 

2.1.1 Deeiqn congideratlone 
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Many of the points made in the following sectlone have been gleaned 

from the problems encountered by previous workers, notably, Kind^^, 

33 28 * * 

Jones , Englar and in particular finglar and Williams^®. 

2. 1 . 1.1 Aerodynamic 

Three-dimension li effects dominate the aerodynamic considerations when 
designino a circulation control aerofoil model. They arise from three 
sources: 

i) high lift testing of large models in closed tunnels. 

ii) effects of low, undefined aspect ratio (of the order of 1) and 
endplate boundary layer interactions due to the large adverse 
pressure gradients at the aerofoil trailing edge. 

iii) spanwise irregularities in the slot flow. 

The recommendations of Englar and Williams'*'^ cover many of these 
points and suggest that the cherd to tunnel height ratio should be 
less than 0.3 and that some form of secondary blowing should be 
incorporated to re-energise the spanwise extremities of the flow. 

The secondary blowing could take the forms (shown in Figure 14) of 
either 

I) blowing on the end plates aft of the half chord point 

II) blowing between the leading edge and the half chord tangentially 
over the aerofoil surface. 
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111) blowing with excess momentum compared to the main jet at the 

spanwise extremities of the main trailing edge slot. The' c are 
called 'tip jets', (see Figure 14). 

For reasons of blowing air economy and manufacturing simplicity, the 
tip jet option was chosen. 

To reduce the boundary layer growth on the end plates, they were 
positioned with less than a half chord projecting in front of the 
leading edge of the aerofoil. The position and dimensions of the end 
plates is discussed further in Section 2. 1.1. 2. 

The reason for the low aspect ratio is primarily an attempt to 
increase the chord and hence the size of the trailing edge wall jet. 
The bigger the wall jet, the fewer the problems encountered when 
velocity profiles and turbulence surveys are required. A compromise 
must be drawn between this and the overall tunnel constraints. 

Jones outlines qualitatively the problems involved in obtaining a 
reasonedjle solution. 

The thickness of the boundary layer just upstream of the blowing slot 

and hence the chord, also indicates the range of slot widths which 

may be used. For example, if the boundary layer is too thin, then 

manufacture of narrow, accurate slots (h < 0.25 mm (0.010 inch)) 

becomes extremely difficult, while if the boundary layer is too thick, 

there is a chance of causing shear layer separation in the wall jet, 

41 

as found by McGahan . It was proposed in this work to vary the slot 

height : chord ratio (h/ ) around the optimum value of 0.002 proposed 
40 ^ 

by Englar et al . The slot Itself should be free frcxn obstruction 
for at least 2.5 cm (1 inch) from the slot exit (this limit allowed 
the majority of the contraction to be 'clean') and t-'.o slot flow 
shovild exhaust as near tangentially to the outer surface as 
physically possible, with the minimum cross section occurring at the 
exit Itself. 
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Th« slot chordwli* poiition was chosan to ba batwaan 96 and 96.5% x/c 
In ordar that tha slot axit wax cloia to tha point at which tha 
unblown boundary layar would normally aaparata* this wax axpactad to 
be the moat efficient location. Tha final poxltlon wax calculated ax 
a function of the geometry of the trailing edge and the optimum xlot 
height « and wax located at 96.45% x/c, ax xhown In Figure 15. 

As a further aerodynamic consideration the model geometric Incidence, 
Oq, was made adjustable by +10° relative to the tunnel centreline. 

2. 1.1. 2 Mechanical 

The two main problem areas concerned the blowing air and the model 
mounting In the tunnel. 

32 

Following the experiences of Kind , It was noted that the aerofoil 
should not be allowed to deform under blowing pressure. This placed 
limitations on the model construction and in particular, on the choice 
of materials. Other researchers have found that if the plenum air 
pressure is allowed to act on retaining screws in the outer skin, the 
screws tend to loosen and leakage and loss of structural Integrity 
may occur. 

The vertical mounting of the model In the wind tunnel working section 
(to reduce chord : height ratio) and the size and poxltlon of the 
access panels (see Section 2.2) Imposed restrictions on the end plate 
design. In particular. In allowing the end plates to extend less 
than a half chord In front of the leading edge, the axis of the model 
and the end plate size were fixed by their positions relative to the 
roof and floor access panels. It was also a requirement that the 
model should be easily removable from the working section to allow 
the tunnel to be available for other testing and to permit easy model 
maintenance. For this reason the model was mounted on Its own 
'floor' and could be lifted Into place as a complete unit Incl'jdlng 
control InstruBMntatlon, on a simple hoist system. 
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2. 1.1. 3 Conitructional 

The modal overall dlnanaions ware restrictad by tha capacity of tha 
available workshop machinery. This was a major factor in arriving at 
the final design scheme. 


Care was also taken over the choice of materials, with particular 
reference to corrosion. Wherever possible aluminium or brass were 
used, but where steel could not be avoided, the components were 
cadmium plated. 


To further simplify production it was decided that the leading and 
trailing edges should be constructed from d,. t steel tube of a 
suitable diameter. 

2. 1.1. 4 Operational 

Consideration was given to shortening the actual tunnel running time, 
by optimising the model operation at the design stage. It was 
assumed that, since the investigation was concerned with detailed 
trailing edge surveys, variations of geometric Incidence were of 
secondary importance compared with the accurate positioning of a 
probe (wire, film or pressure) around the trailing edge, without 
stopping the tunnel or the blowing air supply. Por this reason, any 
probe traversing system had to be operable externally from the 
working section and should not change its effect upon the flow field. 
Further dicussion of these points is included in Sections 2.1.4, 3.2, 
3.3. 


Following consideration of the above requirements, certain aspects of 
the aerofoil were finalised irrespective of span. 

i) An elliptic 20% section with 24 inch chord was chosen to match 
previous work, provide a relatively thick section for ease of 
manufacture and to reduce the curvature at the leading and 
trailing edges. 
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II) The leading and trailing edges should be made from 2 Inch 
(50.6 mm) O.D., 10 G wall, drawn steel tube. 

III) The positioning of the tubes should be such that all external 
surfaces should join tangentially. 

Iv) The leading and trailing edge .ilots should be symmetrically 

placed at 3.55% and 96.45% chord. These figures were defined 
from the assumed geometry, an optimum slot height and slot lip 
thic)(ne88 and the condition of surface tangency. 

v) Ihe slot width h, should be variable by ±0.015 inch (0.375 mm) 
about h = 0.035 inch (0.875 mm) giving h/c - 0.CXH5 and 

" 28.5, vdiere is the radius of the trailing edge. 

vi) The radial traversing gear should be housed within the trailing 
edge tube, rotatable from outside the tunnel and have a 
resolution of better than ±0.001 inch (±0.025 mm) . 

vii) Blowing air supply to the model plenum should be from both ends 
of the model. 

vlli) No pressurised surface should deflect more than 0.005 inch 
(O. l25 mm) . 

ix) The slot contraction should be free from obstruction for at 
least 1 inch (25.4 mm) ahead of the slot exit. 

The first model design was of 54 inch (1.37 m) span, aspect ratio =2.2, 
mounted vertically in the wor)clng section. Three inches (7.S cm) were 
allowed on Jooth sides for removal of the working section wall boundary 
layer and the model was to be of rib and skin aluniniun construction. 
This design was abandoned mainly due to production difficulties (for 
long lengths) and a smaller span model (27 inch, 0.685 m) of simpler 
construction was adopted. This consisted of an inner steel box 
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Pl.n^ with wocd.„ c»tour bitch, to .« th. out.r profll. voo. It. 
m.oy .dv.nt.,.., th. pl.„o„ could b. built .nd ...l.d, .p , clllc.il, 

M .ll«ln.t. led..,., without th. con.tr.lot of h.,l„, . d..lr.d 
..rodynamlc prolll., hl„ th. load upon th. block, w.. ,mi„l..d 
•Inc. th.t. w.r. only ..rodynanlc lore, upon them. Thl. mod.l had to 
b. n.odlfl«i bacauB. of warping dlfflcultl.. In manulacturlng th. 
wooden contour blocks. Th. final production ood.l ».. . conprool.. 

b.tw..n th... two d..lgn. and 1. th. l.ttcr pl.„u„ a..l,n with a rib 
and skin outer profile (Section 2.1.2). 


2.1.2 Final scheme 

Th. final production modal (a.. Pig^:., 15 , ig, con.lstwi of a 
box planum constructad from h Inch < 6.34 mm) aluminium alloy plat, 
with a rib and skin outar proflla. Th. planum g.omatry was maln- 
talnad by 16 small mtarnal riba, which also located th. Lading and 
trailing adg. cyllndars. Aluminium spaers, >. Inch (12.7 mm) O.D 
war. also used to add rigidity to th. planum surf. can. Th. pl.num 
was divided in two by a spanwls. main spar, separating th. Lading and 
trailing mig. blowing plenums, six V Inch ( 9.53 mm) aluminium rib. on 
the top and bottom plates provided th. contour tor th. <, Inch (3.18 mm) 
aluminium alloy outer skin which was held In place by 6 BA countersunk 
head bolts, at k inch (19 mm) pitch. In th. ribs. Th. Inner surface. 
Of th. plenum w«r. coated with non setting adhesive (Uon llguld 
jointing) before aesmnbly to provide ...ling .g«i„.t .ir L,k,. m. 
plenum stagnation pressure and tmaperatur. were m...ur«l via a pitot 
tub. and a thermocouple mounted in th. middle of each pl.num. 

Th. slots were made by two numerically machlnml slot lips (s.e 
section 2.1.3) which loc.tmj directly onto th. plenum dumber plates. 
The slot gap was adjustable «>d could be lockml by 24 locating bolu 
In each Up. The lower surface was ...Ld using two strip, elmll.r 
to th. .lot Up., machln«l to fit flu.h agalnat th. cylinder., 
sealing between th. Ups, s„Ung strip, and th. pLnu. wall. ... 
achl.v«l by 1/16 inch ( 1.6 m) thick, .oft rubb.r strip. ..cursd 
bstween th«. "Ssalastlk’ ... .l„ ussd betw.„ th. waling strip. 
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and tha adga cyllndars to provide additional aealing. The elot lipe, 
aealing etrlpe and leading and trailing edge cylinders were cadaiun 
plated to protect the slot contraction against corrosion. 

The aerofoil was supported by two spigots attached to the end ribs 
and rotating in thrust bearings housed in the support assemblies. 

The support assemblies and leading and trailing edge cylinders which 
extended to the tunnel walls, were left unfaired. This ensured that 
the unfaired cylinders would not produce any lift and therefore could 
not produce the asymmetric interference encountered by Kind^^. 

End plates were fitted and they located into the end outer ribs and 
rotated with the model. The end plates measured 34 in x 54 in 
(0.06 m X 1.37 m), the longest side parallel to the aerofoil chord, 
and projected approximately half a chord in front of the model 
leading edge. The end plates were wire braced to each other and the 
tunnel walls to provide additional stiffness. 

The final overall model chord was 23.425 in (0.595 m) . 

There were 44 static pressure tappings located around the model 
centreline (see Figure 17) . One of the tappings proved unservice~ 
able due to being drilled incorrectly. Tubes from the static 
tappings were routed through the outer ribs to avoid contact with the 
high pressure plenum air. 

Details of the blowing air supplies are given in Section 2.3 and of 
the pressure data reduction system in Chapter 3. 

The aerofoil incidence was sat by rigidly mounted incidence locking 
arms located near the tunnel walls. 

2.1.3 Slot design 

Tha importance of good slot design was highlighted by the variable 
performance obtained by other researchers. Mills the main sndel was 
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bting d««l 9 nad a imaU alot geometry modal waa built, and thie ia 
ahown in Figure 18. It conaiated of a 2 in (50.8 mm) ariuare aection 
tube, with optical glaaa aidea and contained adjuatable linera. The 
tube waa conatricted by an offaet 2 in (50.8 mm) diameter cylinder, 
creating a variable height alot. The linera could be interchanged 
to represent various slot geometries. A vertical splitter plate was 
inserted to chec)t that slot flow Interference was minimal and that no 
significant wake was present over the cylinder. Liquid film, smoke 
flow visualisation and schlieren were used to check on the flow within 
the slot. This simple model was used to establish the production 
slot geometry. 

The production slot geometry is shown in Figure 19; details of the 
location and sealing are also shown. The geometry was initially 
drawn by eye and dimensioned later. 

Both the slot lips and the lower surface sealing strips were 
produced on a numerically controlled milling machine, this being the 
only complex machining involved in the final model design scheme. 

The slot lips were designed to have an optimum geometric position at 
h » 0.035 in (0.875 mm) and contraction ratios in excess of bo were 
typical. The actual slot lip thickness was 0.010 in (0.25 mm), this 
being a precaution against lip deflection under aerodynamic or 
blowing load. 

2.1. 3.1 Two-dimensionality of slot flow 

Many researchers, notably Gartshore and Newman^® have indicated 
severe problems arising from slot flow which is not two-dimensional. 
Figures 20, 21 show the alot flow for this model with a slot gap of 
0.021 in (0.52 mm) measured with a flattened pitot tube and a multi- 
tub# rake. Apart from tho various surface and plenum irregularities, 
the uniformity of the flow is excellent. The slot gap was constant 
to within ±0.001 in (0.025 mm) across the span and the slot lip had 
been machined to a sharp 90® edge to ensure a strong fixed 
separation from the slot lip. The excess tip jet blowing is clsarly 
seen (in Figure 20) at the edge of the main jet flow. The 
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infliction* in th« total praaiur* diatribution at th* boundary of th* 
main and tip j*t* may alio indicat* th* pr***nc* of a sicondary 
longitudinal vortan causad by *xc*siiv* tip J*t blowing a* diacuiiad 
in Saction 2.2. 3.1. Apart from th* inflactiona, th* tip j.t blowing 
ia ahown to have little harmful effect upon the main jet blowing 
diatrlbution. 


2.1.4 Traversing gear 

It waa required to be able to position a probe at any position within 
1 inch <25.4 mm) of the surface of the trailing edge cylinder, 
mainly on the centreline, although some spanwlse movement was 
desirable. The accuracy was required to be within 0.001 in (0.025 mm) 
radially and on 0 (angular position from slot). The 0 position 
was to be adjustable from outside the working section and the probe 
supports should offer a minimum of flow disturbance; 

There were three possible arrangements for the traversing gear: 

1) hou.«i fr^ th. moa.1 with the prcb. brought to th. 

trailing edge, (Kind^^) . 

ill hcu..d lnt.rn.lly within th. pl„„ or tr.lllng edge cylinder 

«.d projecting th. probe redlelly through th. .urf.ce, (Jon..^^), 

Ill) hou.«J within th. telling edge cylinder with th. prob. returning 
to th. .urf.c. .om. dl.t.nc. .w.y fron th. tr.ver.lng ge.r. 

m. third option w.. cho..n. In view of th. probUm. .ncountered on th. 
flr.t two .y.t«.. by pr.wlou. r....rch.ra, wd .!«, b.c.u>. of It. 
I»pro».d rigidity «d lnt.rf.r.nc. propertl... The on. probl.. with 
thl. option wu thet th. g..r .y.t«. ... to b. houMd within • tub. of 
1% in (94 ibb) internal diameter. 
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2 - 1 < 4 . 1 Mechanical 

» d.t .11 of tho final layout of tha travaralnj gaar lyatam 1 . .h 
Fi,ur. Th. driv. unit choaan ... I 

tfP. 9904 111 0400 With a .tap , 03 ,. ^ 

-L ~tl~ - — 

-3^.at.d inoiudlno. tact and Piuionr^U^rtH^^^ 
c«na, in,antu.Uy a . 0 ™ gear ay, ten. .aa choaan aa thla offered 
■aver, advantagea, including,- aoall aiae for large hunt in tar 
^a.u«ion, the teeth are in ccnatant oeah, cheap, leaa auaceptihle to 


TO alginate bachlaah in the ayaten, t«, featurea were incorporated. 
Firstly, the mesh of the gears was ad-tu«»<-*Ki 

contact with th. u , adjuatable to keep the wcrn, in firm 

ntact .1th the .heel. Secondly, a aprung collar ..a fitted to the 

drive shaft inside the .omi .heel to keep the threads 1 

way contact irrespective of the direction of ootlon. 

The two probe supports were *'ecriilr«a,s 4-m 4 j 

rrroK . equired to give adequate rigidity to the 

Lnle " ^ -) Pitch. 


The rotation of the traversing gear aa a unit .as achieved via a fi.ed 

eft on the .,1. of the trailing edge cylinder, e.tendlng beteeen the 

aver, ng gear body and the .tapper motor mounting block, since the 
Stepper motor was at a fixed rei*M«.a i 

unit did not affect the a,” a^u^g eyT: 

Stepper motor mounting plate was held i The 

V4 . position by grub screws 

which clamped onto the trailing edge cylinder. 


The resolution of the system was 


ss follows) 


1 pulse of motor ■ 7^30' 

50)1 worm gear reduction - 9* 

32.6 t.p.i. on .haft . 0.0307 inch., per revolution 
.. ipula. - 1.279 a 10-9 inch.. 

•• approklmately 75 .tep, - o.ool in (0.025 .) 
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Th« motor could b« driven either romotoly by the PDP-11 comput«r or via 
a manual control. Tha position of tht probs was rscordsd on a countsr 
as a numbtr of stsps from a zsro datum and was also availabla as a BCD 
output from the stepper drive. 

To obtain maximum advantage of the accuracy of the traversing gear a 
reliable zero setting procedure was required. This is described in the 
next section. 

2.1.5 Position control system 

Previous researchers have used either a mechanical or optical, wind off, 
zero setting device for the. positioning of hot wire probes. These 
systems have many disadvantages, the main one being the deflection of 
the probe under airflow, and hence can only be considered accurate to 
within ±0.001 in (0.025 mm) at best. 

To obtain improved accuracy, the capacitance position control system, 
shown in Figure 23, was developed. 

Simply, the probe was used to detect the capacitance between Itself and 
the Surface of the model. The model as a whole was energised with a 
5 volt r.m.s., 5 KHz a.c. signal; the charge detected was amplified and 
passed through a phase sensitive detector set at the same frequency. 

The final output voltage was compared with a preset reference voltage 
which enabled inhibition of the stepper drive system. A typical output 
for measured voltage against distance from the surface is shown in 
Figure 24. The voltage produced was a function of probe orientation 
relative to the local surface, but for any run the orientation was 
constant. Using this system, the probe could be positioned optically 
wind off and the zero returned to, reliably, wind on, even if the 
drive system miscounted the number of step pulses or the probe was 
moved to a new angular location ai-ound the trailing edge. There were 
four main advantages with this system t- 

1) the electric field is unaffected by airflow end thm tha probe 
could be positioned whilst the tunnel was running. 
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li) wh«n cloM to th« surface, the repeatable resolution was better 
than O.OOOl in (0,C»25 mm). 

lil) any wire or film probe can use this system. Any insulated 

metallic pressure probe could also use the system, although this 
was not investigated in the current work. 

iv) a probe can be positioned over any metallic part of the model. 

using this system and a travelling microscope, a single sensor probe 
can easily and repeatably be positioned at its minimum distance from 
the surface dictated by the prong tip diameter. 

The potential of this system for probe positioning over the entire 
surface area of the model was not realised at the initial design stage. 
Hence a thin brass strip was Inserted around the trailing edge, on 
the centreline, to enable probe positioning on the model centreline. 
During the experimental work however, the anemometry probes were 
positioned over the steel cylinder away from the pressure tappings to 
avoid the majority of interference effects from the surface 
discontinuities . 


It should also be made clear that the anemometer probes were electrically 
isolated from the model. This was achieved by supporting the probe 
holders in non-conducting nylon bushes and Isolating the model from a 
common earth. Th.e probe connecting cable was also isolated from the 
fluctuating field by the addition of a further braided wire sheath 
around the cable. The field was automatically inhibited when the 
probes were activated. 

The system has been shown to be accurate enough to be used to test for 
vibration of the probe mounting and probe supports. Thm probe was 
driven into and out from the surface and the output monitored on a 
storage oscilloscope, (see Figure 25). The results show little 
detectable vibration, wind on or off. 


Th« following ■•ctlon lUuitratta th« um of th« poiitlon control/ 

travarsing gaar/anamonwtar lystami to perform a limpla boundary layer 

traveree. 

i) ewitch stepper control to 'Manual' and 'Ref. pt' , 

11) mount probe and align with surface. 

ill) move probe towards the surface, monitoring the gap with a 
travelling microscope. 

iv) when the probe Is at the desired zero position, adjust the 
reference voltage to match the probe output voltage. Mote, all 
metallic equipment other than the model and probe should be 
removed from the wor)tlng section as they affect the field 
around the probe and hence the reference voltage would be 
incorrect. 

v) set step counter to zero. 

vi) move probe away from the surface emd switch to 'DiSA'. This 
enables the probe to Ise heated and the zero flow voltage 
measured. It should be noted that the proximity of the model 
may significantly affect the measured zero flow voltage. It 
was found necessary to determine the variation of zero flow 
voltage against distance from the surface so that a true zero 
flow voltage could be determined at any known position from the 
surface. 

vii) switch back to 'Ref. pt.* having switched off the probei start 
up the tunnel. 

viil) whilst monitoring the probe output voltage, drive the prcAie 

towards the surface. If the probe stopa due to reference point 
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voltafl., than the probe he* deflected toward* the *urfac*r lero 
the counter. If the probe *top* due to aero on the counter, the 
probe ha* deflected away from the aurface, allow the counter to 
90 through zero until the reference point voltage i* encountered 
and then zero the counter. Typical deflection* due to airflow 
were of the order of 0.001 in (0.02S mm). 

lx) the probe is now positioned at the original set point, switch 
to 'DISA' and activate the probe. The probe can now be driven 
away from the surface and the anemometer output voltages 
recorded. 

X) having completed a traverse, switch off the anemometer (note: 
tunnel still running) and switch the stepper control to 'Ref. 

pt.'. The traversing gear can now be moved to its new angular 
position and locked in place. 

xi) repeat viii) , ix; , x) 

Using this procedure the probe could be accurately reset close to the 
surface, taking account of small surface irregularities and deflection 
of the probe due to airflow. 


2.2 The 7 ft X 5 ft Wind Tunnel 


The wind tunnel used was the Bath University large dual purpose tunnel, 
shown in Figure 13. The model was mounted vertically in the 7 ft x 
5 ft high speed working section a* shown in Figure 26. The tunnel i* 
powered by a 170 h.p. motor which drives a 4 bladed fan, 10 ft ( 3 . o'? m) 


in diameter. The high speed working section has a 4tl contraction 
upstream, fitted with a gauze, and downstream there is a multicell 
diffuser. The maximum continuous centreline velocity in the high 
spend section is 160 f.p.s. (49 me"^) . Vhm return duct of the tunnel 
is a 12 ft X 10 ft industrial working section with a maximum centre- 
line velocity of 40 f.p.s. (12.3 m* M. A three component balance 1* 
mounted above the high speed working section but this was not used as 

the model was mounted verUcally to obtain the maximum tunnel height: 
chord ratio. 
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2.2.1 Calibration of th» 7 ft x 5 ft working ■tetli>n 

Prior to tasting tho aurofoil, tha working faction was racalibratad to 
taka into account tha affact of tha aarofoil upon tha rafarance atatic 
praaaura diffaranca tapping locatad on tha aida wall at tha antranca 
to tha working aection. A aariaa of new tappinga waa poaitionad in 
the ceiling at the end of the contraction and the atatic praaaure 
difference across the contraction was recorded with and without the 
model in place. The static pressure difference was measured on a Betz 
micromanometer and referenced to the dynamic pressure at the gauze 
measured by a pitot static tube and an alcohol micrenanometer. The 
static pressure tapping which suffered least model interference for 
maximum contraction static pressure difference was chosen and a 
calibration of Betz manometer readings against working section dynamic 
pressure obtained {Figure 27) . The third tapping from the start of 
the working section was found to best satisfy these conditions. 

The tunnel was operated with the low speed working section vented to 
atmosphere to avoid air inflow through the alternative vents in the 
high speed test section. 

During initial commissioning runs of the aerofoil it was apparent that 
the flow in the working section was severely disturbed and an initial 
rake survey showed the flow field to be asymmetric as shown in 
Figure 28. On examination of the tunnel, the contraction gauze was 
found to be partially blocked with dirt causing a recirculating 
separation bubble in the contraction. The gauze was cleaned and the 
calibration and rake survey repeated. The flow was shown to be 
greatly improved (Figure 29) . 

2.2.2 Investigation of free stream turbulence 

During the calibration of the high epeed working section, the free 
stream turbulence levels were investigated using a dual sensor DXSA 
hot wire probe. The results (Figures 30, 31} clearly show the 
disturbed nature of the lower flow prior to the cleaning of the gauze. 
1*he free stream turbulence in the undisturbed region is less than 1% 
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and can ba conaldarad laotroplc. Tha anamomatar waa oparatad as 
datcrlbad In aactlon 3.2.4, 

2.2.3 Propoaad taatlnq conditions 

All expariroents ware carriad out at a constant fraa stream velocity of 
31.9 ms 1 given by a reference pr ensure difference of 60 mm of water 
on the Betz micromanometer. The velocity was kept at this modest 
value to enable a better comparison with the work of Kind^^, and also 
to ease the noise problem for the operator during the anticipated long 
run times. These conditions gave a test Reynolds number of 
approximately 1.3 x 10® based on the aerofoil chord. 

The slot blowing rate was seen to have a slight effect upon the 
reference pressure difference but it was small compared with its 
effect upon the downwash corrections. At each test point, the tunnel 
r.p.m. was adjusted to give a constant reference pressure difference. 

The low speed return section of the wind tunnel was vented to atmos- 
phere, providing a free stream stagnation pressure slightly less than 
atmospheric . 

^**atmos ■ < 1 cm of water 


This proved to be of advantage in the determination of the radial 
static pressure distributions (see section 3.4). 

2.2.3. 1 Effect of tip jets 

Following the recommendations of Englar and Williams^® discussed in 
Section 2. 1.1.1, small tip jets were installed at the extreme spanwlse 
positions of the trailing edge blowing slot. The basic idea was to 
blow the tip jets two or three times harder than the main jet and 
effectively disrupt any vorticity which might be shed at the junction 
of the aerofoil and the and plates. The vorticity present can be 
reduced by careful filletting between the end plate and the model skin 
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In addition, tha larga adveraa praaauro gradlanta which axlat at the 
trailing edge of the aerofoil have a very pronounced effect upon the 
end plate inner surface boundary layer. 

An attempt was made to ascertain the effectiveness of the tip jets by 
measuring the flatness of the spanwlse, upper surface, half chord 
pressure distribution. Eight static pressure tappings were incorporated 
at the 50% chord station, equlspaced to within l«s in (76 mm) of either 
end plate. It was found however that little non-linearity could be 
detected with or without tip jet blowing. 

The overall effect of the tip jets, as shown in Figure 32, was to 
increase the value of by an amount which increased in proportion to 
the amount of tip jet blowing. Since no plateau or peak was indicated 
in the production of extra it was decided to attempt an 
optimisation using flow visualisation. It was found that a single 
wool tuft positioned over the splitter between the tip and main jets, 
extending over the trailing edge, proved most successful. 

When no tip jet blow was used, the main jet caused the tuft to deflect 
towards the centreline of the model. The amount of tip jet blow could 
then be adjusted until the tuft was parallel to the endplates, 
indicating no crosbilow in the stream. This is suggested as one of 
the best methods of ensuring two-dimensionality in circulation control 
testing, particularly at low free stream velocities. It should be 
noted that the tip jet blowing system became less effective at 
increased slot heights due primarily to the nature of the non- 
adjustable sealing strips between the main and tip plenum chambers in 
the slot contraction. Also, if the tip jets were blown too hard, a 
secondary vortex was created between the tip and main jet flows, 
ve-introducing three-dimensional effects. 


The optimum amount of tip jet blowing was not obvious but is thought 
to k)e a function of the main jet momentum, the overall circulation and 
the geometry of the end plates. 
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2.2.4 Blockage correctiona 

A simple solid blockage correction was applied to the free stream 

dynamic pressure since the effects of the high lift coef f icients » 

model mounting « air supply pipes and the impingement of the strong 

wake on the multicell diffuser were difficult to determine individually. 
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The correction applied was from Pankhurst and Holder and it was 
assumed that the model spanned the tunnel, to take account of the end 
plates and mounting tubes. 

The correction is given by 


( 2 - 1 ) 

where - t 2 (j) 

X “ 0.822 for a closed tunnel 

t _ maximum aerofoil thickness 
h txinnel height 

X • *j(l + ^) for an ellipse 

c _ aerofoil chord 
t aerofoil thickness 


This gave a blockage factor of 1.0275 for the free stream dynamic 
pressure. The teat conditions were then 

- 613 Nm"^ 

throughout the experiments. 





2 . 3 Th« Blowing Air Supply 


2.3.1 Df cription of syitam 


A ■chomatic diagram of th« blowing ayitem i> shown in Figure 33. Ths 
blowing supply system was based on a Beilis and Morcom, oil-less, 

3 cylinder compressor delivering 1300 c.f.m. P.A.D. up to BO p.s.l.g. 
±10%. The compressor delivered the air direct to a settling tank. A 
self regulating Area valve installed at the outlet from the settling 
tank, provided initial regulation of the air sv^ply. A galvanised 
pipe was installed, external to the buildings, between the settling 
tank emd the wind tunnel, and a Hattersly stop valve was provided at 
the entry to the wind tunnel building. To improve hot wire/film life 
and reliability a Dcminick Hunter 0.1 um filter was fitted iiranediately 
downstream of the stop valve . 

The flow then passed through a manifold into three separate channels, 
each controlled by its own Hale Haunilton RL6D dome valve. Each RL6D 
was operated by an L15 controller, powered from a tapping just down- 
stream of the filter. This provided a system with independent control 
for the leading and trailing edge slots and for the tip jets. Each 
flow, regulated by an RL6D, was then split into two to supply each aide 
of the model. Shortly after splitting, the flow for the main slots 
was transferred to 2 in (SO. 8 mm) I.D. braid reinforced nylon tube of 
40 p.s.l.g. maximum working pressure connected to the model manifolds. 
Orifice plates (see Section 2.3.2} were positioned in the flexible 
hose for the main leading and trailing edge flows only. The tip jet 
flow passed to the model through \ in (9.5 mm) bore nylon tube. 

Ihe maximum pipe flow Mach number was kept below 0.45 to reduce losses i 
the worst possible conditions were evaluated thus: 




norpt Maximum 

flow caaa(c.f.m) Mach number 


6" pipe 
4" pipe 
)" pipe 
2" pipe 


settling tank to atop valve 
stc^ valve to 3-way manifold 
3-way manifold to 2-way manifold 
final supply to model 


1300 

1300 

1300 

650 


0.1 

0.22 

0.395 

0.44 


The final case was equivalent to nearly three times the flow required 
to choke the main trailing edge slot. 


It was not considered necessary to accurately monitor the tip jet flow 
and the dome regulation pressure was found to be sufficient for control, 


Calibration of orifice plates 


The orifice plates fitted in the upper and lower, leading and trailing 
edge main jet supply lines were designed to BS1042^^. since the 


required operating flow range was large and the connecting tubing was 
unable to comply with the recommendations of BS1042^^, in terms of 


minimum straight lengths either side of the orifice plates, a 
calibration was made. 


A series of nozzles was manufactured to the general design shown in 
Figure 34. These could be fitted in turn to the upper leading edge air 
supply pipeline, the lower leading edge hose being blanked off. By 
ensuring that the flow through the nozzle was choked, the fo3 lowing 
relationship could be assumed: 


m 


"nozzle “ 0-585^ (Ij** ^g 
where a - eree of exit (m^) 


( 2 . 2 ) 


- itegnetlon teny>erature (®K) measured 
by a thermooouple upstream of the noisle 

prwtsure in Bara, measured 
by e atethem pressure transducer upstream 
of the orifice plate 

Y»B ere the usual gas constants, for air 
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Thr.. nozzlM W.r. m.d. with ilff.tznt zxlt .r..z, .nzblln, .id. 
coverage of the flow range. 


With th. m... flow, „d obtalnzd .. .bov, ,„d the przz.wr. 
dlff.r.n=. zoroz. th. put.. ™..„.d on . =o..truct.d w.t.r 

manometer, the calibration was obtained asi- 


m 



■ 0.146 



(figure 35) 


(2.3) 


wher^=! Ah is the pressure difference 
in cm of water. 


across the orifice plate measured 


It w.. Ut.r found that ther. wa. a diff.tenc. betw..n th. flow, of th. 
upp.r and lowar .uppli.., but. zinc. th. calibration appli.a to all flow 
rates, it applies equally well to all four individual .uppli... 


BI°" I-nq supply control and in.truo«nt.Ho„ 

With th. blowing sy.t™ coupled to th. model, it wa. found pozzlble, 1„ 
th. mas, flow calibration, to replace th. stagnation pr...ur. up.tr.ma 
of th. orlflc. pl.t, by th. model plenum stagnation pressure, with 
negligible error. Th. .y.tem wa. thus controllml primarily by 
monitoring th. plenum pressure on a Statham 0-25 p.a.i.d. pre.sur. 
transducer. Although .mall variation, in Ah and T^ occurted from day 
to day, in using constant p^ settings, th. ov.r .11 repeatability of th. 
value. Of CU was found to be better thm. 2*, th. transducer sensitivity 

being adjustmj accordingly. Pull detail, of th. calculation of Cu are 
given In Section 2.3.4. 

Alternatlvaly, the upper leading edge blowing supply could also be uamS 
a. th. supply tor the «,.nomet.r cellbratlon rig (see Section 3.2 mid 
Plgur. 36) . Dependent upon the nozzle .1.., , .imii.r accuracy for 
nozzle exit velocity az for Cp could be attained. 

A. mantloiiad in Section 2.3.2 a dlff.renc. ewlstad betwam th. upper 
and lowar modal pi.nizt suppllaa, due to the ditfarance in length of 
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pipe run*. Thi* w*s alio th* ca*« on th* tip j*t lupply. Hh*r« th* 
main ■uppllt* axhauitad into a naar stagnant planum, tha tip jata 
maintatnad tha asymmatry through the blowing alot. For thi* reason, a 
small restrictor valve was placed in tha shorter tip jet supply line 
enabling equalisation of the blowing rates. Thi* valve had to be 
adjusted according to the slot height and unfortunately slightly 

reduced tha capability of the tip jets to control three-dimensional 
effects. 


Serious instabilities were encountered on the RL6D dome valves in the 
form of a resonance between the valve lip and seat just as the valve 
opened. As the valve 'cracked' , a region of high suction was formed 
between the lip and seat causing the valve to close and then open, 
once more, under the dome pressure. The resultant ' trun?>etlng' was so 
severe as to prevent testing of the aerofoil. After extended testa 
the only satisfactory solution was to increase the valve spindle 
damping by filling the spindle support cavity with Tate and Lyle 
'Golden Syrup'. All oils tested lacked sufficient viscosity and all 
greases failed to flow back to the closed position. The application 
of syrup not only stopped the resonance but also reduced the low 
frequency variations in blowing supply from sundry causes such as 
interactions with the Area valve. 


2.3.4 Evaluation of the blowing momentum coefficient. Cu 
The blowing momentum coefficient is given by 


m V. 


»|PV^S 

w 


(2.4) 


now for two-dimensional tests 




*m._ v_ 
2D J 

J|pvi C 


where 


2D 


mess flow per unit spen 


(2.r>) 


Vj ■ Jet velocity 
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Th« valu* of tha jet velocity wae obtained from a amall iterative 
procedure baaed on an iaentropic expanaion in the nozzle. 


Since the atagnation temperature and preaaure in the plenum are known, 

a value for o can be found, 
o 



( 2 . 6 ) 


Aaeuming that the flow is incompressible, an Initial value of the slot 
Mach number can be found. 

( 2 . 7 , 

Pq ^ ^ 

Due allowance was made for the reduction in slot span due to the 
presence of the tip jets, when relating the orifice plate calibrated 
flowrate to the aerofoil slot flow. 


This now allows a first estimate for p, since 



( 2 . 8 ) 


The procedure can be repeated and the value found for which eolvee 
the loop. This takes no account of any losses or boundary layer 
growth within the slot. 


The value of the plenum stagnation pressure is found relative to the 
unblown slot exit static pressure. The error caused by not referring 
to the actual slot exit static pressure in the presence of the jet was 
small since low blowing rates were used and the free stream dynamic 
preasura was small cowparud with the plenum stagnation pressure. 


Iinitf biiiiiiiiiiLiiiuii i: i i ! 


M A check on the celculatlon procedure, the trailing edge jet 
velocity wae determined by three independent methoda. 

i) Aeeuming ieentropic conditions at the nozzle givea, 


D \Y - 1 


' Vd- 


(2.9) 


^ ^ 2Y f / <» 

“lean " (Y - DRT^ iUjj 


( 2 . 10 ) 


where ■ duct temperature (®K) 

-2 

■ duct preeaure (Nm ) 

Pjiij ■ free stream static pressure (Nm”^) 


The duct variables were assumed to be the plenum stagnation quantities 
The slot contraction ratio for this test was 70 : 1. 


A single sensor hot wire probe was placed in the middle of the 
jet and the velocities calculated from its calibration. The 
probe was moved 0.001 in (0.025 mm) either side of the centre- 
line of the slot and no velocity change was apparent. The heat 
loss to the model had a large affect on the performance of the 






ill) The Cu calculation program could b* uaad f.o datarmina tha Hot 
valoclty. 

The reaulti of the three methoda are ehown In Figure 37. The agreement 
la aurpriaingly good. The deviation of the Cp calculation could be due 
to loeaea in the a lot, while the apparent increaae maaaured by the hot 
wire could be due to boundary layer build-up in the alot and probe 

blockage at high speeds, .»nd to model Interference effects at low 
speeds. 


The slot height for this test was 0.021 in (0.52 mm) and the agreement 
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3. INaTRUMBWTATIQW 
3 . 1 PrMiurt Data Aequlgition Svi f m 

A block diagram of tha ayatam ia ihown in Pl;fura 28, Two 46 port *j' 
larlaa Scanivalvaa wara uaad, allowing altarnata tapping! on altarnata 
Scanivalvaa to provida tha aaquanca for tha atatic tappinga. Thia 
anabled each tranaducer to return to ambient praaaura before the 
experimental preaaura waa applied, thua avoiding praaaure lag problema. 

The Scanivalves were controlled by an online PDP-0 computer which 
stepped them in sequence and analysed the data. At each measuring 
point the transducer was allowed a settling time of 250 and then 
the pressure at the tranaducer face was found as the average of 
50 samples over a 50 ns period. 

The transducers were a Statham +2.5 p.s.i.d, and a Setra +5 p.s.l.d., 
they were identically calibrated over the same range using the low 
pressure transducer calibration system (see Section 3.1,1). The range 
of the calibration was determined to give the maximum signal 
permissable for the tests envisaged within the 10 volt A/D converter 
range of the PDP-8 computer. The derived transducer sensitivity was 
then input as a program constant. The maximum pressure coefficients 
allowed were ±5.5 based on a free stream dynamic pressure of 613 Nm”^ . 
This range was only possible because of the low blowing rates being 
usedj previous researchers had encountered pressure coefficients in 
excess of -20 in the trailing edge suction peak. 

The analysis of tha data was performed, as shown in Figure 39, by a 
Scanlvalve control program on the PDP-8, discussed further in Appendix 
1. Briefly, points were interpolated midway between adjacent data 
points using a Langrangian Interpolation technique. At the leading 
and trailing edges, where high curvature of the dietributlons was 
present, certain raatrictions ware placed on the interpolation routine 
to reduce instabilitiea. These had a small effect upon the overall 
P*i^^®>^®anoe results. The lift coefficients and pressure drag 
coefficients were evaluated using a simple numerical integration 
technique . 
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It !■ usual to avaluAts ths normal fore# cotfflclants including 
allowanoaa for tha raactlon thruat of tha jat. 


s ■ * s*‘"“ 
s " s ■ 

m ^ 


(3.1) 


where M Is the angle of the jet exit to the horizontal and the sub- 
Bcript m indicates the initial measured quantity. Since this work was 
mainly concerned with the detailed trailing edge investigation, these 
corrections were not included in the computer program. The effect on 
the lift coefficient would be small since Cp was generally less than 
0.03 and p was small. The effect on the drag force Is much larger) 
however the initial accuracy of the drag integration was questionable 
due to the lack of static pressure tappings around the trailing edge. 

The results, including the interpola-.ed values, were displayed on a 
CRT display at the tunnel console, enabling an immediate visual check 
on the flow. 

The aerofoil pressure distributions were used to evaluate the effective 
incidence of the aerofoil and this will be described in detail in the 
next section. 

j-l*l Evaluation of the aerofoil effective incidence 

Since an aerofoil with a rounded trailing edge does not exhibit the 
well known 'Kutta' condition of conventional sharp trailing edged aero- 
foil sections, the determination of the effective incidence is complex 
and subject to error. Kind used a method of comparing pre-drawn 
potential flow pressure distributions with the experimental results. 
This enabled the downwash correction to bo calculated and the effective 
incidence deduced. This method has been adopted and improved in this 
study, by matching the experimental and theoretical pressure distri- 
butions over the leading half chord. Thu leading half chord was chosen 
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in nrdtr to avoid tht major diraot affaota of tha blowing jat which 
occur round tha trailing adga. 

Th. thMrttlcl pr..tu» dlitrlbution .round .n ,lUp,., ,i„„ pj, 
Jonaa for axampla, ia 



( 3 . 2 ) 


0 < 


0 


a - effectiva incidence, a 

eff 

< 27T such that J . ^j(i + CO10) 

o ■ 


THU., It e i. d.fln.d at th« .tatlc t.pplnp po.ltlon. th. th.or.tic.l 
preiaure coefficienta at theae pointa could be determined for a given 
\ff' The experimental leading edge half chord C waa 
calculated aa a fir at atep. The value of the eatimated waa input 
to the Scanivalve control program by the operator. The program then 
evaluated a aet of theoretical pointa and adjuated the value of the 
theoretical until the integrated areaa of the theoretical and 
axperimental leading half chord diatributiona were equal. 
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A co»nparl»on of th* maaaurad prtitur* dlitrlbution with th« pot.ntial 
flow .olutlon wai ihown on a CRT dlaplay in tha wind tunnal. Tha 
oparator than 'vlaually' itaratad to obtain tha baat fit to tha 
praaaura distribution ovar tha laadlng half chord by auccaialva 
astlfnataa of tha affactlva incldancat 

The three major problema with this technique werej 

I) the model geometry differed slightly from that of the true 

ellipse, due to the circular arc leading edge and the forward 
facing blowing slot. 

II) the effective incidence could only determined to ±> 3 ° at any 
single data point. 

lii) the technique could not be used when the leading edge was blown 
since the forward Jet flow caused too great a disparity between 
the theoretical and experimental distributions. 

Examples of the displays produced by this technique are shown in 
Figure 40. Further discnasion of the accuracy of this simple 
theoretical calculation is given in Section 5.1.1. 

3*l-2 Low pressure transducer calibration 

This system was for use where the required pressures were less than 
±30 in (0.76 m) of water, it consisted of a small pneumatic actuator 
which was coupled to the transducer and a water manometer. This 
simple system proved to be very reliable and accurate over its small 
range. Also it was simply adopted to enable both Scanivalve trans- 
ducers to be calibrated at the same time, thus avoiding any sensitivity 
errors between them. 
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3.1.3 High prg« «ur» traniduetr caUbratlon 

Thi. >y.t«m wai originally producad for uaa on tha auparaonic wind 
tunnel facility at Bath Univaraity. it conaiatad of a naadla valve 
controller which allowed aimultanaoua praaaura aupply to two 
'calibration atandard Wallace and Tiernan gaugea, giving differential 
and absolute preaeurea, and to the required prea.ure tranaducer. The 
supply was vented to atmosphere via a short section of fine glass 

capillary tubing which allowed the system pressure to stabilise 
quickly. 

The whole system was fed from the main blowing supply via the tip jet 
valve controller. This enabled the maximum calibration pressure to be 

set, avoiding any possible over pressures being applied to the trans- 
ducer face. 

The maximum range of this system was 80 p.s.i.g. 


All transducers were calibrated using either of these two systems and 
showed negligible hysterisis. An example of the calibration of the two 
Scanivalve transducers using the low pressure system is shown in 
Figure 41. 


It is important to note that the Scai.ivalve transducers were calibrated 
in situ as it was apparent that the operating temperature within the 
Scanivalve housing had a definite effect upon the calibration. The 
zero of these transducers was also found to be affected by the securing 
•enslon of the transducer retaining caps in the Scanivalvea. 


3.2 Hot Wire Anemometry 
3.2.1 Calibration rio 

All investigations ard u.iibrations of the anemometer equipment were 
performed on the nozzle rig shown in Figure 36. This is a similar 
arrangement to that described in Section 2.3.2. for calibratino the 
orifice plates. By use of different size nozzles, velocities from 
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w to 300 could b. produced, „c„l.. „„ .11 produed 
to th. ■«» 5.oin.trlc form .houn In figure 34, 

Tbec.llbr.tlon rl, u..d filtered .Ir, .. mentioned In Section 2 , 3,1 

t.°b le"Tb‘’'°‘* • -“tin/ 

teble Which eupported the v.rlou. probe conflgur.tlone, 

*b exp.rlment.1 Inve.tlg.tlon of the c.llbr.tlon Jet w.. performed 
using . disa du.1 ee„„r wire probe, the operation of which 1. 

scuseed In Section 3,2,4, The results obtains! showed the extent 
sod development of the constant velocity core and the turbulence 
level, within the Jet, These teste Indicated that a probe could be 
centrally positioned up to 1 inch ,25,4 mml from the nor.le face 
^Ithout css Of jet velocity. The results are sho. l„ 42, 43, 

e orifice Plate calibration dlscuased In Section 2,3,2 was used to 
determine the jet velocity at the nozzle. 

me ability to position a probe away from the face of the norzle 
lncre.sed the rotational limit of the rotary table by some 4o“ , This 

sTcI ir'^blT csUbratlon of the split film probe 

since It enabled the flow angle to become negative for the chosen 
probe geesnetry (see Section 3.3). 

gg,"sral comments on hot wire anemom.f.y 

tn general th. hot wire, and the associated 1 ns trmn. station were 
operate! In accordance with the relevant operation manual, 

dwldth at a typical expected flow velocity and «. upper limit 
freguency response of 40 K«s wa. attainable (... ri,„r. 44> , 

s“i:“ t.i.tio„- 


■ E ^ + Bu" 
0 o 


where B and n are the calibration conat«ita. 


(3.3) 
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The conetent n wa« to be measured ae it was felt thr.t the vai'je of 
0.5, usually assumed, was unreliable over such a wide velocity range 

It can be shown that the effects of ambient temperature drift are 
confined to a variation of B, the intercept of the log plot, rather 
than the slope n. Krietensen^^ proposed a solution to this problem, 
which has been extended by Starsmore^^ and as a final exact solution 
(see Appendix II for details) can be shown to be 



where 

O ■ overheci: ratio 

a «. ilfluid ~ '^cal^ 
"^fluid 

T measured in *K 


This relationship not only allowed correction to calibrations for 
ambient temperature drift but also allowed the calibration temperature 
to be corrected to a single value if the run temperature varied. The 
effect of varying the required calibration temperature at a fixed 
overheat ratio is shown in Figure 45. 

The standard equation can be reduced to:- 



and by plotting the log of both sides, the values of B/E ^ and n can 
be found. ° 


in. followln, corwction to Ui. r.ro flow volt.,. ... 

to t.k. .ccoimt of th. lnflu.nct of convKtiv. cooling. »lth th. air. 
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aligntd in the vertical plane i 


Eq “ 0,92 X E 

o measured 


( 3 . 6 ) 


An .xfnslve Inw.tlj.tlon w., p.„o».d Into th. .£,.cts of flow 
tompornture upon the overhe.t ratio c.llbr.tlon of wire probe. 
The ueuelly adopted probe eettln, procedure 1. outlined below, 


i) 


zero probe cable resistance on anemometer bridge. 


ii) measure probe resistance r 


iii) set 


operating resistance on anemometer decades. 


ro - (1 + a) r^ 


(3.7) 


iv) operate probe in flow field. 


Throughout this worJt, the di«?a 

, cne disa anemcaieters were unable to zero the 

cable resistances and consequently the resistance 
the operating resistance was assumed to be 


used to calculate 


“ ‘'probe + cable “ *^cable 


(3.8) 


The value of could be found from the resistance decades as that 

value necessary to give zero cable resistance. 

'ir "" "“>• v„i.tion. in 

^ient fluid tebperetur. while oper.tin,, then .llow«.o. ehould eleo 

BOvino*flTd^* '•■Iztance varying with the 

. ..lu. of „ ei,nific«,tly d.cr....d with incree.in, owerheet retio. 
Since the overheet retio i. . function of the dif,.r««e between the 
probe operetin, te,*.r.ture end the «a.ient fluid tebpereture, it 
-mma edvieeble to ev.lu.te the operetinp r.ei.tvic. fro. . probe 
r..i,t«.c. »..ured in the pr...„o. of th. «n,m, fluid. n,i. ... 
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thought to be the only way of maintaining constant overheat ratio anA 
hence constant calibration, over a period of time. 

Equation (3.8) now becomes 

’ Vobe + cable ‘ “^cable 

In the experimental arrangement already described, variations of fluid 
temperature from +20*^0 to -5°C were common, dependent upon the outside 
weather conditions and the air mass flow rate, whilst the aiiJsient 
tunnel air varied from +10°C to +35 depending on how long the wind 
tunnel had been running. The Importance of the revised technique is 
apparent and tests indicated an improvement in the repeatability of 
the calibration constants for a given wire from +5% to +1%. 

A typical wire resistance temperature coefficient, a, given by 

r - r^(l + 01 (T, - T^)) (3.10) 

1 o 1 o 

was found to be 

at = 3.695 X 10“^ n per “C 

Hence a 10 variation in ambient temperature would change the 
operating resistance by nearly 0.07 fJ at an overheat ratio of 0.8. 

This would have a significant effect upon the zero flow voltage and 
the probe calibration in general. 

The reason for the variation of calibration with overheat ratio is not 
clearly understood but may be a function of the change in thermal 
stress of the wire between the prongs caused by the change of operating 
temperature. The easiest solution was to adopt the above procedure 
and rvin at a constant overheat ratio of 0.8. This overheat ratio was 
chosen as a compromise between frequency response u,d probe life. 

A further series of experiments indicated little variation in the 
calibration constants over extended periods of time. The reason for 
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thli ii thought to be th« 'clean' air supplied by the 0,1 pm filter. 

The calibration equation (3,5) was initially used in these testae 
however, the small variations in n that have previously been discussed 
caused large variations in B because of the remoteness of the data 
points from the vertical scale. It was found more reliable to use the 
following equations 



Examples of the two techniques are given in Figure 46. The values of 
n obtained from both techniques were consistently close to 0.4 and 

the improved confidence in the intercept from the second method is 
obvious . 

The calibrations were all obtained by a linear regression method 
programmed on a PDPll/34. Typically, greater than 99.9% correlation 
was achieved and it was found sufficient to take only 5 or 6 
experimental points to determine the calibration constants. 

Boundary layer wire probes (5 pm diameter) and DISA D series bridges 
were used throughout the single sensor experiments. 

Two data reduction techniques are available for hot wire anemometry. 
The unlinearised signals can be analysed by either an analogue or 
digital technique. The analogue technique is particularly suited to 
single wire operation and will be discussed in Section 3.2.3, while 
the digital technique offers many advantages where more than one 
sensor is used (see Section 3.2.4). It was decided not to use 
^^^**^^***^* in this work since the controlling parameters ware 
susceptible to ambient temperature variations evident as changes in 
the wire calibration constants. 
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3.2.3 Operati on of ainglt ««nior wira prob«« 

All alngla aansor probai wera oparatad on tha almpla analogua ayatam 
shown in Pigura 47. Tha reduction of tha turbulence data and datar- 
mlnation of the mean velocity, including the effecta of ambient 
temperature drift, were performed on a pocket programmable calculator 
using the inverse calibration equation and the recommendations of 
Mojola^*. 



where a ” speed of sound ■ 



ms 


“1 



(3.12) 


(3.13) 


where Is the r.m.s. of the fluctuating portion of the bridge 
signal and the bar signifies a time averaged quantity. 


Extension of the simple analogue technique to dual and triple sensor 
probes is possible but requires the use of summing and multiplying 
circuits. Pull detailB of the necessary relationships are given by 
Mojola^^. 


An example of the effect of yaw angles in the plane of the prongs and 
normal to the wire is given in Pigure 48. It is clear from this 
figure that the contraction effect of the prongs whan the flow is 
normal to the prongs and wire produces approximately a 10% increase 
in the effective cooling velocity of the flow. Ifhereas when the flow 
Is angled in the plana of the prongs, the effective cooling velocity 
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!■ a function of the alna of the yaw angle, repreeenting the effective 
length of the wire in yaw, 

3,2,4 Operation of dual eeneor wire probee 

Whilet it was not expected to use these probes for investigations 
around the trailing edge of the aerofoil due to their lack of spatial 
resolution, their operation was investigated to gain experience in 
the use of dual sensor probes. Each sensor of an X~array wire probe 
was calibrated normal to the airflow in exactly the same manner as a 
single sensor probe . 
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The digital analysis technique used by Starsmore was adopted for 
this work since it was simpler and more adaptable than the equivalent 
analogue dual sensor technique. It consists of sampling the output 
of both wires at discrete time Intervals and calculating the 
instantaneous velocity vectors for each time step. It has been 
shown that the instantaneous velocity vectors for an X-array dual 
sensor probe in the XY plane {see Figure 49 for co-ordinate system) 
are: 


V - 


_ 1 _ 

l2 


\w * ■’21 ^ 1 

'i =1 

I’l * «2 1 

2 

1 2 21 
1^1 - **2 1 

2 

1 (1 t a^ / 

(1 + a^) 


(a^ - 1) 




2(a^ - DU 


(3.14) 


where “ effective cooling velocity wire #1 
* effective cooling velocity wire #2 
a ■ direction eeneitlvity coefficient 
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i 


A Jull derivation of thaia aquatlona ii given in Appendix ii. 


If a total of n aamplee ia taken at a frequency f, then the time 
averaged mean quantitiea are given by 


u - - 

n 


n 

E 


n*l 

n 


( 3 . 15 ) 


V - - 
n 


E 


n»l 


At ,.ch sample point, th, tluotuatln, valocltU. than glvan by, 


u' ■ U - u 
n n 


( 3 . 16 ) 


V* - V - V 
n n 
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Thu. th> tlmt avtragtd maan and fluctuating quantltUa art .imply 
d*ducd. Pull dttail. of th. program d.v.lop.d for thi. mthod ar. 
giv.n in App.ndix ll. Th. procdur. can^ r.p..t.d by reorientating 
th. prob. in th. XZ plan, to produe. w, w' ^ and u'w'. 

Th. valu. of a, where 

, , meaeured velocity 

actual velocity ' flow parallel to the wire 

ia generally accepted to be 0.2. This was ahown to be a reaonable 

aaaumption compared with the results measured in this work (see 
Figure 50) . 

The sampling frequency f, is desired to be as high as possible, or 
the sample size should be increased in order to 'capture' the 
maximum turbulence contribution. The sampling rate of the PDPll/34 
A/D converter has a nominal maximum value of 3.4 KHz using the supplied 
system library routines, in practice this figure can only be 
approached dependent on the type of storage file management employed. 

A maximum rate of 2.6 KHz was attained for this study. 

A total sample size of 2560 samples per wire was finally used as a 
compromise between accuracy and analysis time. This sample size 
showed negligible error in turbulence meas- ament compared to on. of 
32000, and gave an order of magnitude decrease in analysis time. 

Block diagrams of the dual sensor analysis system and the related 
computer program, are given in Figure. 51, 52 (and ... also Appendix 
11) . The filtered fluctuating components of each bridge signal are 
amplified to the maximum ±2.5 volts parmissabl. by th. PDPll/34 A/D 
converter. Th. amplified fluctuating input, war. monitored on an 
oacilloscope to avoid saturating th. A/D converter. Th. bas. signal 
r.m.a. and time averaged mean voltages, input to th. computer, were 
then used to recreate the actual voltage, and hence the velocity, 
time history, in this manner the maximvim fluctuating signal 
definition was obtained. 
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Thii data raduction tachnique waa uaad to Invaatlgata tha tunnal 
fraa straam turbulanca lavala and tha calibration Jat davalopmant 
(aaa Sactiona 2.2.2, 3.2.1). 

^•2-5 Sgurcaa of arror in tha hot wira anainoniatrv maaaur*m.n.., 

3 • 2 < 5 . 1 Tji_a tam parature qradiant aeroaa tha wall -ja*- 

Aa haa alraady bean atatad in Section 2.3.1, tha blowing air aupply 
tamperatura waa very dependant upon tha outalda weather conditiona. 
It waa uaual for the jet air temperature to be of the order of 10°C 
cooler than the surrounding tunnel air. Having indicated the 
dependence of a hot wire probe calibration on ambient temperature, 
the importance of the techniques previously described (see Section 
3.2.2) to account for these variations becomes apparent. However 
since the wall jet produces a strong entrainment effect between the 
jet and the external flow, the existence of a temperature difference 
between the flows produces a temperature gradient across the flow 
which varies with downstream posiitlon, the cold blowing jet air 
gradually being heated to the ambient conditiona. 


Since the temperature gradients were indeterminate for this study, 
the effects were allowed for by using the jet temperature at slot 
exit to adjust the wire calibration, in this way the error was 
reduced for the high velocity raeastirements, but some discrepancies 
were apparent in the wall jet edge velocities compared with the split 
film measurements (see Figure 98, 99). Had the tunnel air temperature 
been used, the edge velocities would be more realistic but an equal 
percentage error to the previous low velocity measurements would have 
produced increased velocity errors in the jet velocity measurements. 
Short of attaching a small temperature sensor to the wire probe and 
measuring the actual temperature at each measurement point, the 
adopted method seems to be an acceptable compromise. 


An attempt was mads to measure the temperature gradient across the 
flow by detecting variations of the resistance of an unheatcd wire 
probe. However, lacit of time for the necessary adjustment of the 
anemometer bridges prohibited precise measurenmnts. it was also 
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osnild.r.a um»ctM«ry .ina. th. bulk of th. tt.llln* .dg. .urv.y mi 
to bi pirforiMd uilng . .put film pcob. which ippirintly m. 1... 
•uiciptihli to inblint timpiricuri virlitloni (lu sietlon 3.3.1). 

Non-radial travwrtlno oath 

Du. to th. flnlt. .It. .„d po.ltlotilng of th. tr.wir.lng g..r prob. 
.upport. wd th. d..lc. to pl. 0 . th. hot wlr. prob. .. «,r to th. 
•lot .. po..lbl., It w.. n.c....ry to off..t th. wlr. from th. r.dl.l 
tr.«r.lhg plm. (... 5 .,,. th. prob. w.. tt.v.r..d io 

th. conflgur.tlOT .hown, «. ungul.r off.,t .t th. .wt.nt of tr.v .1 of 
.pproxli«.t.ly J typlcl. in th. pr..,n« of longltudln.1 
velocity gr.ai.nt., error. 1 „ th. m...„r.».nt of th. w.ll ).t outer 
layer Vvaiocitlea were apparent. 

The angular error was given by 

0 

* " — j ^*** figure 53 for notation) (3,18) 

0 was typically less than 10 ® which gives an error of less than 2 ° 
for R/y < 5 . 


It was possible to position the probes such that the offset angle was 
aerog however it was felt that the advantages in being able to approach 
the slot exit with the probes outweighed the errors produced. 

3 . 2 . 5 . 3 »^i and probe int erfe rence effects 

Th. met U.U.1 con.ld.r.tlon wh.n u.lng wlr. prob.. in do., proximity 

to . .olid W.11 1. h..t 10.. to th. .urf.e. It i. u.u.lly ,..um.d 

th.t th. h«t lo.. to th. .urf.c. de. not b.co». .l,nlflc«.t until 

th. prob. 1 . 1 ... th.n 0.004 in ( 0.1 .) from th. .urf.c How.v.r, 

in thl. .tudy, vdoclty m...ur«Mnt. .. clow .. phy.lo.lly .llow»3 

by th. prob. prong r.dlu. (0.003 In, 0.075 mm) w.r. .ppirontly 

un.ff.ct.d. Thl. wu thought to b. du. to th. unu.u.lly high .h..r 

gr.dJ.«it do., to th. w.11 nduclng th. nlitiv. Inorw.. m nM.ur.d 
valocity. 


The probe eupport geometry (eee Figure •5 3) wee deeigned to minimi ae 
Interference at the meaauring point. The Introduction of normal 
obetruction in the curved wall jet tendi to 'aplit' the flow and 
generate three-dimenaional inatabilitiea . It waa aaaumed that the 
■mail projected area of the prong tlpa would cauae minimal dietur-banco 
of the flow. The affect of probe interference upon thu measured lift 
coefficient is discussed in .Section 4.1.2, 5. 4. 1,2. 

The results obtained in Section 3.2.3 for the effects of yaw in the 
plane normal to the wire indicate that negligible error would arise 
from small changes in orientation relative to the curved surface. 

This conclusion also implies that the normal velocity measured on a 
single hot wire will also include the V contribution. 


where b is the direction sensitivity coefficient in the plane normal 
to the prongs, usually equal to 1,1. Hence in regions of high flow 
angle relative to the surface {due to entrainment and jet growth) , 
the hot wire will tend to overestimate the velocity. The wire probes 
were, in general, orientated slightly angled to the surface ('2 or 3°) 
to Improve the definition of the reference point for the position 
control system. 

To summarise, the maximum error in the wall jet velocity profiles 
occurred at the outer edge of the wall jet boundary layer and could 
be as high as +5 ms ^ dependent upon flow conditions and probe 
position. 


3.3 Split Fj]m Anemometry 

Split film probes are a relatively new type of probe marlteted by 

47 * 

Thermo Systems Inc. , and consist of two semicylindrlcal separate 
platinum films, sputtered onto a 0.006 in (0.15 mm) quartz rod. They 
offer a much improved spatial resolution (of the order of 8 times) 
compared with a conventional cross wire sensor and are approximately 
half the focussed volume of a laser anemometer system. A boundary 



layer type probe le aketchod in Figure 54. 

The probe operate! on tlie principle thati- 

i) The total heat tronefer on both filme gives a measure of the 
velocity vector perpendicular to the sensor. 

ii) The difference of heat transfer for the two films gives a 
measure of the velocity vector perpendicular to the plane of 
the splits on the sensor. 

Hence, from i) assuming a Kings Law type relationship, 

+ k^E_^ " A + B U " (3 

This assumption is valid, since if the splits are considered small 
compared with the total film surface area, then the total heat 
transfer would follow a law similar to that of a single wire. 

From ii) 

Ei^ - - f(Uj^)sin6 (3. 

where 0 is the angle between the plane of the splits and the normal 
velocity vector. 

Therefore 


U > U COS0 
n 


V - 0 line 


(3.22 
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The constant K is described by the manufacturers as a correction for 
non perfect matching of sensore. In the limit, K must be equal to 
the ratio of the zero flow voltages end T.S.I. s.iggssts that adjust- 
ment of the overheat ratios should be used to set K to ^;nity. In the 
light of the experiments with hot wires, it was felt to be more 
reasonable to adjust the probe cold resistances on the bridge 
resistance decades to obtain K ■ 1: this enables a finer control of 
the value of K. It should also be noted that in order to minimise 
the offset from the actual resistances, as one film is increased in 
resistance so the other should be reduced, or vice versa. This means 
that the total resistance is unchanged and the total heat loss due to 
normal velocity is also unchanged. In this manner, k is assumed 
equal to unity from here on. 

Originally, it was hoped to operate these probes from two 55D01 DISA 
anemometers: however, one of these anemometers was found to be 
insufficiently stable to run the films concurrently. The films were 
therefore driven by a Prosser 6100 twin channel anemometer system. 

This system has a high/low bandwidth switch which simplifies the 
running of these probes. The low setting limits the bandwidth to 
10 kHz while the high setting allows adjustment similar to that on a 
DISA anemometer. It was found that instability of these probes, due 
to thermal couplings across the quartz rod, was their major drawback. 
Thiree different modes of behaviour were experienced : - 

i) When using the DISA anemometers, since one channel was less 

stable, only 1 film could be run at a time. The large thermal 
couplings caused one film to drive the other and the situation 
was reversible but always of a 'hard over' nature. The 
phenomenon occurred even when the anemometers were run at the 
lowest gains and bandwidth with an overheat ratio of less than 
0.1, Independently of which film was started first and the 
duration of the probe cut-in time. 

11) Using the Prosser system with the low bandwidth selected, a 
stable situation occurred and the probe operated normally. 
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ill) With tha Proasar ayitam lat to high Isandwldth, a aituatlon 

ocourrad whara both films wara running and fllm^2 was found to 
ba slavad to fllm^l i.a. for an Incraasa in dua to a changa 
in flow angla^ Incraasad by tha saroa amount. 

The results of the above observations was that tha probe had to ba 
run at a reduced bandwidth. T)iis showed however, only a 2% drop in 
detected r.in.s. readings compared with the high setting at the chosen 
overheat ratio of 0.5 emd this was considered satisfactory. 


^.3.1 Calibration 

The calibration of these probes requires that the probes be set at a 
series of angles (in this case -5 to +30°) relative to the plane of 
the splits and the bridge voltages recorded for a series of speeds. 
This provides the data necessary for all the constants to be derived. 
From equation (3.20) 




+ BU 


n 


where 


E ■ E. measured x 0,92 

1,2 ^ 1,2 

to account for the effects of convective cooling upon the zero flow 
voltage . 

The constants B, n can be derived in a similar manner to those for a 
single sensor hot wire, usjng only the zero angle data. 




+ nlogU (3.23) 
n 
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No .llo«mo. for th. v.ri.tlon of tho anbUnt t«np.r.tgr. was 
inoludad In tha caUbration, If a similar oorraction to tha wira 
proba waa uaad, tha Indooad aff«:ta wara found to b. far mora 
algnlfloant on raaulta. Thla la thought to ba du. to running at a 
lowar ovarhaat ratio, o . 0.5, conparad with that uaad for tha wlra 
probaa and to tha large thermal oapaclty of th. probaa. Tima did 
not allow a full Invaatigatlon of this affeot and good rapaatablllty 
waa achlavad without th. Inclusion of a tamparatur. oorraction. Thla 
will be shown in Section .3.3.3. 


It ahould ba noted that, without further Invaatigatlon, It la not 

racommandad to run thaaa probaa up to tha aama wind apaada aa wlra 
probes for two reasons 


i) structural integrity 

11) tha frequency of vortex ahaddlng from the proba bacomaa 
significant relative to the effects of turbulence. 

Tha author la aware of axparlmanta at up to sonic speeds (Boeing Corp", 

private communication) but la unable to comment on th. reaulta 
achieved. 


The angular calibration constants are somewhat more difficult to 
obtain. From equation (3.21) 

- f(U^) sine 

0 

t - 0 , Ej ahould equal Ej, thla la arranged by netting K to unity 
and la also dependent upon th. actual apllt geometry relative to th. 
datum probe axis. The manufacturers suggest that, 

f (U„) - C U “ 
n n 


(3.24) 
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The syitero for derivation of thee», conetante C, m le auggeeted at 
followe . 


i) Plot and Ej veraua 6 for each speed. 

These are straight lines for a fixed and the crossover pointy 

0 la the deviation of the probe setting from the plane of 

error 

the splits. 0 should be a constant for a given probe 

unless the probe deflects under wind load. 

ii) Amend the yaw angles, using the deduced 

_ ^ 2 ^ versus sin(6 - • ^h^® results in a series 

of straight lines for each given speed. The lines should all 
pass through the origin if K has been correctly measured. They 
also show the equation (3.24) to be of the right form, since 
t>e gradient of each line, C , is a constant for any given 
speed . 

iii) Plot the log of the gradients in ii) against the log of the 
fixed speeds. The slope of this line represents m. 

iv) The value of C can be found to give the best fit to the 
experljnental data. It should be noted that extrapolation to 
the axis to give C is not possible as the response of the probe 
changes considerably at lower speeds, as will be indicated. 

Figure 55 - 59 show a typical set of calibration data for these 
probes. The normal and angular calibration for the probe used in the 
experimental work were found to be. 




(3.25) 


and 




9.B56 

n 


0.566 


sin6 


(3.26) 


It Is Interesting to note that the value of the normal index, n, was 
very similar to that typical of a single wire probe. This would 
seem reasonable since, if the splits are assumed small in area 
compared with the total sensor surface area, then the split film 
probe would lose heat to the air in a similar manner to that of a 
fine wire. 

It should also be stated that the validity of this type of angular 
calibration has not been demonstrated at high angles to the plane of 
the splits. Obviously, as 0 approaches 90° to the plane of the 
splits, the calibration loses its sensitivity to angular variation. 
Since this area was not of interest in the current work it was not 
pursued further. 

Also, at very low speeds, the split film sensor ceases to respond to 
angular variations altogether. This is due to the flow at low 
Reynolds numbers approaching the simple Inviscid flow case. It is 
suggested that the use of this type of probe in either, flows of 
high angular variation or, flows giving low probe Reynolds numbers 
should be regarded with caution. 

The calibration of the split film probe was not found to vary with 
time by any significant amount indicating a resistance to contam- 
ination and oxidation, problems associated with wire probes. It is 
however advised that the probes are 'burnt-in' for a short period 
(15 - 30 minutes) prior to establishing a calibration. This technique 
appears to stabilise the films and their quartz coatings to provide 
a more repeatable operation. 


3.3.2 Data reduction 

In terms of data reduction, split film probes can be treated exactly 
like a dual sensor wire probe. The discussion given in Sections 3.2.2, 
3.2.3, 3.2.4. regarding the use of analogue or digital analysis 
techniques applies equally well. 
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The disital technique la very ainiiUar to that described in Section 
3.2.4 except that the instantaneous velocities of each sample are 
derived directly from the calibration equations. Equation (3.25) 


can be used to give the value of and e can be obtained from (3.26) 
11^18 gives the magnitude and direction of the instantaneous velocity 
vector, allowing the normal components to be derived from (3.22). 

The time averaged of the normal components are found and used 

to deduce the normal fluctuations of each sample point. The normal 
and shear stresses are then simply obtained, similarly to the dual 
wire sensor system. Appendix II gives full details and a listing of 
the analysis program. The digital method was used throughout this 
work. The major advantage of the simplicity of the data reduction 
equations was that the computer analysis time was an order of 
magnitude less than the comparable cross wire system. This enabled 
* real time analysis to be performed. 


The analogue method is inherently more complicated electronically 
since it involves the use of a number of summing and multiplication 
circuits. Since little has been published regarding the use of split 
film probes, a full derivation of the analogue turbulence equations, 
and the assumptions made, is given here for the benefit of other 
researchers^^ These equations have previously been published in 
T.S.I. TB20 but unfortunately the equation for ^ was incorrect. 
The equation was given as 


u'v' 








V2> 


(3.27) 


where e^ and are the fluctuating portions of and E . This is 
obviously wrong, since ^ 




0 


1 


f It 

1 ‘ 


65 


The correct derivation for u*^, and u'v* is as follows j- 


From equations (3.20), assuming K ■ 1 and 


® 1,2 " h ,2 * 1,2 


U - U + u' 


V - V + v' 


then 


"" 2 *“ 2 n 

(E, + ej + (E_ + e.) = A + BU 

11 2 2 n 


n 


u' , aV 


A + BU 1 1 + -£- + _ 

U 2U 


Multiplying out the L.H.S. and expanding the R.H.S. gives, ignoring 
third order terms 


(e/ t ♦ 2(E^e^ + E^.-^) t * e/) 


A+Bu"*Bh" 


lU 2U^ 


2U^ 


now, if the mean flow is aligned with the plane of the splits 


Z 2 ^ - 2 ^ ^ --n 


then 


-n / u' ^ a^v*^ ^ (n-l)u'^ 


2(Ejej^ + E 2 « 2 ) + («i + e, ) ■ nBU (— + — + 


U .3U 


-2 


2U 
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squaring both aides and .ignoring triple products 


♦ V2'] 



therefore, time averaging and taking the square root 



2 

nBU " 






V2> 


(3.281 


To obtain the shear stress equation, consider (3.24). If the mean 
flow is aligned in the plane of the splits, such that 

- 0 


and 

then 


sin 0 * ^ 


- W " <»i' - '2 ’ ■ 

Therefore, by expanding as before 


ra-1 


R.H.S. - Cv’u"’”^| 1 + (m-1) 


aVI 

2U ^ 


2(E^ej^ - E^e^) + - e^^) » CU™ ~ + (m-l)u'v' 


(3.29) 


time averaging 


2 2 , ....rrJn u'v' 

e. - e • (m-l)CU -3^' 

X ^ U ^ 




(m-l)CU 


, 2 2 , 
^ ‘*1 - *2 > 


(3.30) 
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This 1« an int.raatln, comparl.on with tha usual analoflua cross wire 

aquation which is alio a function of tha diffaranca of tha two wira 
rtin.s. valuas. 


If aquation {3.29) is squarad, chan, Ignoring triple products 


Time averaging and taking the square root yields 




cu" 


<Vl - «2'2>‘ 


(3.31) 


These equations are not particularly difficult to solve but they do 
require that the mean flow vector Is aligned in the plane of the 
splits. Any deviation from this is significant since (E ^ - e ^ 
h.. ta.„ ignored „d thi, could be orders of eegnltode gLter'then 
the other terms left in the equation. The digital method is far 
simpler in assumptions and does not require the mean flow to be 

aligned in the plane of the splits, providing obvious advantages in 
a highly curved flow. 


3.3.3 Operation 

■me .put ,11b probe, were povered from . Pro...r 6100 twin ch.nnel 
.n«ncb.t«r .. dl.cu..«l m section 3.3. me probe, were mounted with 
the plene of the epllt. t.ng.ntl .1 to the loci .urf.oe .nd tr.ver.ed 
radlelly u.ln, the tr.cr.lng gear pr.vlou.ly d..orlbed In Section 
2.1.4. me po.ltlon control ey.tem (.ee Section -2. 1.5) wa. und to 
.et the datum laro poaltlon of the probe, rel.tlv. to th. local 
.urfac and the eenaor .* 1 , could be brought to within 0.009 Inch.. 
10.23 mm) of th. .urf.oe. mi. wa. th. minimum dl.unc pcalbl. 
dlctatw) by th. cranio aupport tube, at th. nd. o, tha cnaor. It 
wa. a..um«l that the m.a.ur.m«.t poaltlon coincided with th. ..n.or 
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axil. Tr«v«riai of th« w»U j«t could th«n b« yarformcd ai for a 
wlra proba daacrlbad In Section 2. 1.5.1. 

excellent repeatability wax obtained for velocity profiles around 
the trailing edge of the aerofoil (see Figure 60) Justifying the 
assumptions regarding ambient temperature variations. The data from 
the probes was analysed at each measurementj^int and took approx~ 
imately 10 seconds to give values of U, V, u'^, v 


^ 3 ^ Spurces of errors in the split film anemometry measurements 

Unlike the single wire probes, the split film probes did not suffer 
from being offset from the traversing gear radial axis nor did the 
flow temperature gradient across the flow appear to affect the 
results. However a number of other possible sources of error must 
be mentioned. 

3 . 3 . 4 . 1 Wall effects 

Because of the large thermal capacity of the split film probe con^jared 
with a wire probe (the diameter of a split film probe is approximately 
30 times that of a wire probe) and the orientation of the plane of 
the splits, the split film probe is more subject to errors caused by 
heat loss to the surface. Since it was necessary to bring the probe 
as close to the wall as physically possible, a technique of equalising 
the bridge voltages, in the presence of the airflow, at some small 
distance from the wall was adopted. The technique was similar to 
that for reducing K to unity, (see Section 3.3) and involved small 
changes in the film cold resistances on the anemometer resistance 
decades. These adjustments seemed valid since at some small distance 
from the surface, the local flow angle should still be nearly 
tangential u the local surface. Also, it was attempted to m.^e 
these adjustments at or near the velocity maximum to avoid problems 
of shear flow on the finite probe diameter (see Section 3. 3.4.6). 

This technique obviously is open to criticism but appeared to provide 
good representative results for the time average mean velocities 
down to the minimum distance from the wall. 
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How«v«r a distinct «ff«ct upon th« turbuUnc. maasurtmtnti cloi. to 
th. wall waa alao apparant. Tha haat loa. to tha aurfaca appaarad to 
Incraaaa tha maaaurad turbulanca within a givan diatanca from tha wall 
(aaa Saction 5.4. 3.1). Thia would aaam raaaonabla ainc<* tha praaenca 
of tha wrll would cauaa an aaymmetry in tha raaponao of tha filma and 
hanca an incraaaa in tha maaaurad turbulanca. Prom tha data racordad 
it waa dacidad that all split film turbulanca data within 0.036 in 
(0.9 mm) should be regarded as unreliable. 

3. 3.4. 2 Probe geometry 

From Plgur. 52 It le clo.r that th« plane of the .puts 1. offset from 
the probe support axis by approximately 1°. TOe offset whilst appearing 
positive from Figure 52 would be of the sense to produce a negative 
offset in measured flow angle. 

3.3.4. 3 Plow curvature 

Due to the finite site of the sensor cylinder, the high curvature of 
th. flow produces a movement of th. front stagnation point towards the 
surface. The effect can be shown to be small compared with the 

gecmetrlcal considerations of th. previous sections, being of th. order 
of only 10 seconds of arc. 


^•3.4. 4 Frequency response 

As has already been stated, (see Section 3.3), the split film probes 
had to be operated at a reduced bandwidth in order to maintain probe 
stability. This caused a severe reduction in the detected turbulence 
levels, in some cases by a factor of 4 compared with hot wire results 
Also it is reasonable to suggest that the large relative sise of the 
-put film prube will cause it to be unaffected by eddies of smaller 
else than the probe diameter. This would also appear as a reduction 
of the detected turbulence levels. Example, of these effects will be 
shown and discussed in more detail in Section 5.4.3, 
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Prob« intarf»r«nc 3 « 

flow .ttoet. producd by iotroduoln, . .put m» prob, into . 
hlphly ourvd W.11 j.t with ..vwr. pc.ltlv. „d „.,.tlv. pr...„„ 

With th/flT* d«t«nnln«. in gtn.ral any infc.rfaranca 

with th. now ...... th. J.t to .h..t nnd „u.„p ^.hlnd th. 

ob.ttuotlon. Th. .wttnl to whloh th. fll„ prob. .ff.ofd th. flow 

... hot lnv..tl,.t.d but .11 po..lbl. pr.o.utlon. w.r. t.h.„ 1 „ Z. 

p^. .upport ,.oo.try to .mini., th. .ff.ct., » .i„n.r p„t 

.Uohtly .o.ll,r reduction In lift =o.fflcl.„t to th.t noted In 

film prob*^'!/ observed auggeating that the 

.nd thir* “““ 

and their aupports. 

3.3.4. 6 Shear flow effecta 

we. effectively error .p.tlei reeolutlon. if the probe le 
the pretence of po.ltlv. .hear en effect oppc.lte to he.t lot. to 

we“.“ Ue'drtV'"'" 

«. applied to the reeulta to account for the effeot. of .hear 

gradient! however it is suggested that th« 

uggested that the correction would be of 

ollar form to th.t for a pitot tub. l„ .h..r flow ,... s.«inn 3 . i : 
Uttl. notlc.bl. error wa, produced between the wean velocity pronue 

10™ ctiontl.r •• - 


teJl.l static Preteur. 

in the eajorlty of bouod«ry layer reaearch th. d.teroln.tlon of nornel 

'p L .t:*"” •”- --btion or 

a aJtl "““.ver. In th. c... of a hljhly curvrtflow 

c pressure gradient must exist due to the streamline curvature 
P^tlcular, Where J.t flow .xl.t. and th. velocltle. end eh«r ‘ 

h“* “• “• 

.Ration. b.C0« dominant. I„ p..t r....rch.r. have ..d. 

two basic assumptions regarding the flow.- 
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i) that th« «tra<uttiinai art all locally parallal and concentric 
about the centra of curvature of the local aurface, thuai- 


R ■ R 

o 


(1 


X 

R 


(3.32) 


Thie aeaumptlon will be shown to be invalid in Section 5,4,2, 


ii) that the radial static pressure distribution can be given by a 
simple force balance 


3p PU^ 

m !■ 

3y R 


(3.33) 


Provided that U could be defined as a function of y and that the strear.- 
lines were concentric about some known centre surn that R could ba 
found also as a function of y, this simt.io bale, ice could be of use. 
However neither U nor R is a simple function of y. Many r^searchcLs 
have recorded the overall static pressure diffarenct across curved wall 
jets and have shown somn agreement between exj srim«nt and the liimple 
force balance if a mean Vv'\lue for 0 is absui^iod to apply a :ros!i the 
flow. However, a detailed knowledge of the radial static pressure 
distribution is of great importance if a satisfactory clo'^ure of the 
governing equations of highly curved wall jets is co be achieved. 

28 

Englar attempted to measure tl*e distribution within the jet by 
using a vertical sharp edged plate mounted on the trailing edge, with 
a series of static tappings normal to the trailing edge cylinder. 

These results irdlcatad a departure from the expected pressure 
distribution, Figure 61) and will be further discussed in 

Section 5.3.1, Englar' x technique was thouglit unsstis factory due to 
the interference effects at the junction of the plate and th.i trailing 
edge cylinder. The effect of having a veitlcsl plate within the wall 
jet flow was also suspectsd to be causing large three dimensional 
disturbances. 
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Conwntioniil itatic praiiur* probas cannot ba uaad baoauaa of tha 
highly ourvad natura of tha wail jat flow, tha larga Induoad flow 
anglaa and tha lall icala of tha flow to ba axaninad. 

48 

Pacha duacrlbad a mathod uaing a hot wire proba poaitionad alongside 
a total haad tuba, the atatlc preaaura being deduced from Bernoulli 'a 
Incocnpreaeible equation. Tha method was shown to work for flat plate 
boundary layers and was adopted in this work for a survey of the 
trailing^edge wall Jet. This method was also used by Dvorak and 
Woodward to determine the radial static pressure distribution over 
a slotted flap arrangement. The major problem with this method was 
the relative magnitudes of the terms in Bernoulli's equation. 

■ ^bno.' ■ - P.tno.’ - t ♦ w'^j (3.34) 


Usually, the working section static pressure is set to atmospheric 
pressure by venti g the high speed working section. This implies that, 
Pq (relative to atmosphere) and the dynamic head are of equal 
magnitude. Hence the dif f'<*rencing of two equal large quantities to 
obtain a small quantity is subject to large errors. As mentioned in 
Section 2.2.1. the tunnel used for this research was run with the 
high speed working section vents closed to avoid air inflow problems. 
Consequently, the total head of the airflow was approximately atmos- 
pheric. This Implies a much greater confidence in the values of the 
static pressure deduced from (3.34) since (p - p^^^) was of the same 
magnitude as the dynamic head and the total pressure relative to 
atmosphere was comparatively small. 

A sketch of the system used is shown in Figure 62. The total head 
tube was a flattened stainless steel tube (0.019 in, 0.48 mn thick) 
and the pressure was measured on an alcohol manometer. This type of 
probe is far less seneitive to flow angle compared with a conventional 

static pressure probe. A correction, suggested by Young and Maaa^^ 

aeio 



— “ 0.13 + 0,08 “ 
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(3.35) 
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was applied to taka account of tha affect of the aheared flow field on 
the pitot tube. 



Itie pitot tube and the hot wire were mounted approximately 0.25 in 
(6.35 mm) apart to avoid lateral interference effects and were 
positioned an a common axis using a travelling microscope. The 
position control system described in Section 2.1.5. was used to 
control the datum zero position of the probe combination. The hot 
wire was used as the position detector and the rigidity of the system 
could be estimated since the reference point cut off could be set to 
allow tne pitot tube to just contact tha surface. As the pitot probe 
contacted the surface, so it became energised with the 5 kHz signal. 
This gave a step in the field intensity around the wire sensor and 
increased its positional sensitivity. Relative movement of the two 
probes could then simply be judged by the difference in the datum 
positions as shown on the step counter. The deflection using this 
probe arrangement was shown to Ise negligible. 


3.4.1 Sources of error in the radial static pressure measusrement 


A summary of the affects of possible errors in tha measurements on a 
typical wall jet static pressure distribution is shown in Figures 63, 


3. 4. 1.1 Effect of turbulence 


Equation (3.34) contains an allowance for the turbulent energy in tha 
flow. If it is assumed thati- 


» • 'o - 


and that 


U ■ U + u' ♦ v' ♦ w* 
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then 

• U ^ I (u* + V' + w')^ 

' U U ^ 

Time evereglng glvee 



asBuming the shear terms to be small. 

Then since the measured pressures are automatically time average^', 
quantities, equation (3.34) can be deduced. 

If the longitudinal turbulence intensity ^u' ^ /U w as of the order of 

-1 2 

10% at a velocity of 50 ms , then the correction u* was approximately 
1% of the dynamic pre ssure . Variations of this order can be ignored. 
Similarly even if v' and w* were assumed of the same magnitude, the 
overall effect could still be ignored. 

Turbulence levels of 10% were assumed not to affect the total head 
tube readings. 

3 . 4 . 1 . 2 Eff ect 

The effects of the highly sheared flow are significant only for the 
pitot tub assuming the hot wire to be negligibly small. The existing 
data on the effects of sheared flow on pitot tube measurements are 
very scant and most researchers use the standard corrections proposed 
by Young and Maas^° (see Section 3.4). The effect of this correction 
was to provide a constant effective tube centreline offset in the 
direction opposite to the shear gradient and takes no account of the 
actual magnitude of the shear. This is not considered satisfactory 
but must suffice in the absence of an alternative. For the pitot tube 
used the offset was only approximately 0.002 in, (0.05 mm). Results 
for two different sised pitot tubes showed little difference in the 
measured total pressure. It was therefore asswed that while a 
eorraotiMi should be applied, it had no pronounced effect t^on the 
trends of the results. 
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A third •vn imAlltr prob* war also tastad, thia balng round, aquara 
facad and 0.011 in (0.27 nm) outaida diamatar, but two problama 
pravantad ita uaai 

i) tha axcaaaiva lag tima raquirad to achiava a ataady praaaura 
raading. 

ii) tha proba was inaufficiantly rigid and tha highly curvad flow 
accantuatad thia problam. 

3. 4. 1.3 Errors in hot wire ananometry 

All tha ralavant precautions discussed in Section 3.2.5 were used and 
the anemometry results were considered satisfactory. 

3.4. 1.4 Probe geometry 

The two probes were mounted approximately 0.25 in (6.35 mm) apart 
laterally and the cross interference was assumed to be negligible. 

The probes were assumed not to have deflected under aerodynamic load. 
This seemed reasonable considering the accuracy of the position 
control system for the hot wire and the rigidity of the pitot probe 
support. No deflection was apparent when checked after each run. 

The flow angle data provided by the split film probes showed angles 
not greater than ±10 . Hence, if it could be assumed that both the 
hot wire and total head tubes would provide acceptable data for flow 
angles up to ±10° from the probe axes, then the angular effects on 
the radial static pressure distribution should be negligible. 

3. 4. 1.5 Air density 

The problem was not simply of compreseibility, although density 
gradients must exist within a high velocity wall jet, but of defining 
density value in the flow which had a temperature gradient. 
As a compromise, the density of the air at the Jet exit was used, as 
oaloulatsd by the Cm calculation program (see Section 2.3.4). This 
vsluB was based on the flow rate and ambient conditions. 
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4. PMBBWTATIOW OF RliOLTS 

Thli ••otlon lummariMi th« r««ult« obt«in«d| d«t«ll«d dlacuasien of 
tha raaulta la raaarvad for Chaptar 5, 

4.1 Ovarall Parfomanca 

4.1.1 No blowing parformanca 

Figure 65 ahowa tha lifting parformanca of tha unblown aerofoil 

maaaurad uaing the preasure data acquialtlon ayatam (aaa Section 3.1) 

at a test Reynolds number of 1.3 x 10® based on the aerofoil chord. 

Tha free stream velocity was 31,9 ma~^ including a solid blockage 

correction (see Section 2.2.4). An estimation of tha performance 

38 

obtained from the method of Dvorak is also shown and ia In good 
agreement. 

The measured three-dimensional lift curve slope was found to be 3,25 
per radian. The value for a two-dimensional ellipse at a similar 
Reynolds number frcm Hoerner®^ was found to be 4.0 per radian. Tha 
experimental results contain no allowance for the low aspect ratio of 
the model or wall/end plate interference effects. The tip jet blowing 
system was not used for these tests. 

Examples of the measured pressure distributions for the unblown ellipse 
are given in Figures 66, 67. Theoretical preasure distributions 
obtained from equation (3.2) assuming the measured lift coefficients 
and incidences are shown together with distributions obtained from the 
viacous/potential flow calculation of Dvorak and Kind^®. 

4.1.2 Trailing edge blowing only 

Figures 68, 69 show the overall performance of the aerofoil ae a 

function of trailing edge blowing nosMntum eoeffleient for a range of 

geometric incidences . The affects of slot height, hot wire probe 

interfhrenoe and tip jet blowing are indicated for a range of blowing 
0 

rates at ■ o . The increased slot height shows an improved 
performance, whilst it is obvious that the probe interference nay 
aeeount for a reduction of ae much as 20% in the lift eeeffioient by 
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dliri^ting th« trailing adga auction peak. Tha hot wira proba waa 
poaitionad approximataly 1 m (25 nm) to ona aida of tha atatio 
tappinga , 

Figuraa 70, 71 ahow axamplaa of tha CRT diaplay from tha effactiva 
incidanca calculation (aaa Saction 3.1.1) and alao aoma axamplaa of 
tha maaaur^d praaaure diatributiona in tha praaanca of trailing adga 
blowing. The effactiva incidanca calculation waa uaed to produce tha 
downwaah corrections for a range of and Cu ahown in Figure 72. 

The corrections deduced by Kind^^ are indicated and ahow good agree- 
ment with the positive results. Using th-se corrections, curves of 
against Cp for constant effective incidence are obtained (see 
Figure 73) . Results from Kind^^ and the program CIRCON^^ are also 
indicated in Figure 73, the agreement with the experimental work 
being satisfactory. Some further typical results for different 
ellipses and aerofoil section are shown in Figure 74 and indicate the 
high lift augmentation capability of circulation controlled aerofoils 
compared with jet flap arrangements. 

The measured pressure drag coefficients are shown in Figure 75 for a 
range of and blowing momentum coefficients. The increase in the 
pressure drag due to the high suction around the trailing edge is 
clearly evident. The effect of including the jet reaction component, 
equation (3.1) is indicated. 

4.1.3 Leading edge blowing, only 

Figure 76 shows the effect of leading edge blowing alone i 4 >on the 
aerofoil performance at - 5°. Initially the reduces quickly, 
then a more linear reduction at eWa greater than 0.003 - 0.004 occurs, 
The reason for this change in slope is not fully understood since the 
changeover in flow pattern (see Section 1.2) occurs at a much higher 
Cy. Figure 77 shows examples of the measured pressure distributions 
tax a range of leading edge blowing momentum coefficients. 


78 


^•1*4 LMding and trailing tdo« blowing 

Tht affect of leading edge blowing at aero geometriu incidenoe for a 
range of trailing edge blowing monentum coefficiante ia ahown in 
Figure 70. The expected reduction in performance ie apparent. 

Figure 79 showa axanplea of the meaaurad preaaure diatributione 
around the aerofoil for various araounta of blowing and the two poaaible 
flow caaee are clearly indicated. The firet flow eaae, at low 
and high incidence indicatea that the flow 'folda-back* and flows 
along the upper surface causing a more positive value of Cp at the 
leading edge. The second case, at high Cy-_ and low incidence 
indicates the continuation of the leading edge jet flow onto the 
lower surface of the aerofoil as shown by an increased leading edge 
suction. 

The inability of the effective incidence calculation to match the 

measured pressure distributions in the presence of leading edge 

blowing is shown on the CRT displays. Figure 80 shows the effect of 

dual blowing for a variety of geometric incidences and trailing edge 

blowing rates at a fixed Cy._ of 0.0125. 

L£ 

4.2 Trailing Edge Investigation 

The numerical results of the trailing edge investigation are given in 
Tables 1-17. Tip Jev owing was used at its estimated optlmvm 
throughout these experiments (see Section 2. 2. 3. 7). 

4.2.1 Upstream boundary layer 

The boundary layer approximately 4 slot widths (0.084 in, 2,1 mm) 
upstream of the blowing jet was measured using a boundary layer single 
sensor hot wire, in the presence of the leading and trailing edge jet 
flows. The no-blowing case was also investigated. Values for the 
shape factor of the profiles were determined but are subject to some 
error due to the lac)c of definition of the boundary layer edge . The 
longitudinal turbulence Intensity distribution across the layer was 
also measured. A comparison between the results of Klnd^*, the 
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praaant itudy and the uiual powar law aaawptlon 



(4.1) 


ia alao .hown In Flgura 83. Th. ra.ult. Indicata « 
deviation tram th. power law in the pre.enc. of th. blowing jet 
T7»e leading t .ga flow i. aUo .hown to incr.aa. th. turbulence 
intensity within the boundary layer, as night be expected. 


^•2.2 Wall jet velocity profile. 

Pi,«r.. e. - 91 .how. for o“ ,.c.trlo i„ola.„«, th. vloolty 
(U) profll.. throuoh th. tr.llin9i .19. w.ll J.t for th. four 
t.bul.t.d blowlnp C.W. („. T.bl. 191 and «...ur.d by . .pUt fu„ 
..h«.r, n,. r..ulto .r. ,lv.„ for v.lu.. of 9 .t 10° l„t.rv.l. up 

to Mp.r.tlon, .tMtln, .t 15° fro. th, .lot. In o«.. iuttr- 

mediate traverse, have alao been performed. 


Th. no.0,,1 „ioclty 1. .hown .. th. loci flow c,l, 5/u i„ 

ri9ur.. 92 - 95. Inflow ™,l., .. high .. lo° toward, th. aurfac. 
ar. indlcatad n.ar th. .lot. a. .„tr.i™»nt proo... and J.t growth 
rat. ar. al.o clcrly vl.lbl., Plgur. 96 .how. th. aotual variation 
of V aoro.. th. wall J.t for two blowing rat... Plgur, 97 1. «. 

indication of th. ooopatabUlty of th. r..ult. with th. continuity 
equation 


1 3u . 3v 

R ly ■ ° (4.2) 


Plgare.^98, 99 show comparison of U/U^ .gainst y/y with the wor)c 
of Kind and also direct compari^jn. of the various methods ueed in 
the c-irrent .tudy. hot wire, split filn «nd the velocity data used 
to derive the radial static pressure. The hot wire results are subject 
to the «rror. discussed in Section 3.2.5., »ost obvious of which i. 
the tempsrature gradient across the wall Jet. 
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riguMi 100 - lOS show tht lirtAmlln* pattarnt, Unas of oontMnt 
(y.U/U^), daducad from tha apllt film raaulta for tha four blowing 
rataa. Tha Inflow and Jat growth quantltlaa ahow raaaonabla agraa- 
mcnt with thoaa Indlcatad by tha local flow angla maaauramanta • Tha 
curvatura of the ataamllnaai daducad by a numarical curva fitting 
tachnique and the uie of tha standard aquation 


Radius of curvatura 



(4.3) 


ara shovm for three blowing rates and compared with 



usually assumed in theories and computations. The results clearly 
show that the streamlines are not concentric with the local surface 
centre of curvature. It is interesting to note that at higher 
blowing rates, the distributions tend to be parallel to tha results 
from (4.4). At positions closer to the slot, the outer layer 
curvature is less than expected indicating entrainment affacti.. 

Tha affects of different slot heights and slot lip thic)cnassas are 
shown in Figures 106, 107. The results were obtained using a single 
sensor hot wire probe. 

Tha increased slot height allowed a constant velocity core to exist 
in tha Jat and enabled estimation of tha boundary layer thioknass 
within tha slot contraction. 

Tha increased slot lip thicitnassas ware obtained by sucoaaalva 
addition of adhesive tape (0,007 in, 0.18 mn thio)i) to tha t«>por 
•urfaca of tha slot lip. Intarastingly tha wake of the lip appears 
to altar little in downstream extant. 
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A oomparlfion of Angular position of j«p«»tlon moAsurtd from th« slot 

12 

for ths prsiont study oompsrsd with ths rssults of )(ind ind 

90 

Bnglsr Is shown in rigurs 106. Xt would appssr that tha currant 

26 32 

rasults show good agraaroant with Englar whlla Kind was not 

oparatlng at tha optimum blowing rata for tha tast conditions. Two 

26 

sats of rasults ara shown from Englar « tha first baing tha publlahad 
rasults daducad from a surfaca mountad shaar strass probai tha sacond 
baing takan from tha prassura distributions putlishad in tha sana 
raference. Tha separation point was assumad to ba tha point whara tha 
static prassura reached a constant value , as in tha separation bubble. 

A marked difference exists between the two sets of datb> this will be 
firrther discussed in Section 5.1.2. 

Figures 109 - 112 show the variation of the jet parameter y„/2' ^a' 

U f U . against angular position from the slot for a varitity of 
m min 

trailing edge blowing rates. These results are taken mainl.' from tha 
split film velocity profiles. 

4.2.3 Turbulence results 

All the wall jet turbulence meaeurements were taken at saro geonetric 

incidence. Figures 113 - 116 show the longitudinal turbulence intensity, 

/u^/U , distribution across the wall jet for four blowing cases. The 

measurements were taken at 10° intervals of 6 starting at 15° from the 

slot exit. Radial positioning closer than 15° from the slot was not 

possible with the split film probes and although the hot wires could 

ba positioned almost anywhere, tha error due to radial misalignment 

(discussed in 3. 2. 5. 2) became increasingly noticeable. Tha results 

shown ara from both tha split film probe and a single sensor hot wirai 

they clearly indicate the affect of probe sisa and operational band- 

32 

width. Tha results of Kind ara compared in Figure 117 with tha split 
film rasults at the one data point at which the majority of tha flow 
paramotarn (a.g. angle from slot, Cu) are similar. Figure 118 shows 
tha affeota of varying tha slot height and slot lip thickness upon 
the longitudinal turbulanca intensity an recorded by a hot wire probe. 
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Flgurai 119, 120 ihow tha ralattva pealtiona of U 


„(P 


and Wu' fo>r two blowing oaaaa. Tha OMKlnum turbulanoa 

Intanalty ia ahown to ramaln juat outboard of tha cantra of tha fraa 
a^iaar layar whila tha minimum intanaity cloaaly followa tha valoeity 
maximum in tha jat flow. Tha origin of tha turbulanoa intanaity 
maximum ia ahown to ba tha a lot Up. 

Piguraa 121 - 124 show tha normal atraaa for tha four blowing 

caaaa aa racordad by tha aplit film invaatigation. Tha raaulta again 
ara ahown for 10° intarvala of « through tha wall jat, atarting at 
15° from tha riot axit . 


Pigura 125 ahowa a direct comparlaon between the longitudinal and 
normal r.m.a. fluctuation. The normal r.m.a. ia ahown to be roughly 
60% cf the longitudinal intanaity. Ihe ralativa radial poaitlona of 
tha longitudinal and normal turbulence intancitlaa maxima for a 
variety of Cjj'a ia shown in Figure 126. 


Wje results obtained for the Reynolds's shear stress, u'v* , fr^ the 
split film prcdse are ahown in Figures 127 - 130. An exanple of the 
correlation coefficient, given by 


correlation ooerfieiant ■ 




(4.5) 


ia shown in Figure 131. The usually accepted value is 0.4 although 
this ia subject to error at points near areas of high 8^0/9y* (i.a. 
the velocity maxima and minima) . 

It is suggested that tha correlation coefficient as measured by a 
split film probe can be used to produce results representative of a 
higher banAridth device. If it is assumed thAtha relative 
magnitudes of the turbulence intensities are constant compared with 
each sensor bandwidth, then let 


( 4 . 8 ) 


film 
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Th«n 



■ 1 X 



ALio AMunlng • oonalatant oorralation coaff Iolanta thani- 


u'v'^ir^ ■ /Fx X film oorralation coaffloiant (4.7) 

In t hia mannar it ia poaalbla to produce mora raallatlc valuaa for 
2 

u'v* and v' and to indicate the bandwidth limitationa. Hia raaulta 
of thaaa adjuatmanta are ahown in Flguraa 132, 133 for one blowing 
caaa. Thia technique ia only approximate aince the relative aiaaa 
of the eanaora will affect the ability of a given eddy aiae to cauae 
a fluctuation of the bridge voltage. 

The ahear etreaa raaulta are compared with the aaauned diatribution 
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in the program of Dvorak at al and with the meaaurementa of Wilaon 

21 33 ' 

and Qoldatain and Jonea in Pigurea 134, 135. 

4.2.4 Radial atatic preaaura dietributiona 

4fi 

Aa a check on the validity of the tech.'lque of Pache (aee Section 
3.4), the atatic preaaura diatribution actr a the boundary layer, 
juat upatream of the trailing edge blowing alot waa determined. The 
raaulta, with and without trailing edge bloving, are given in Pigurea 
136, 137. The raaulta ahow excellent agreement with the uaual 
aaauaq^tion within the boundary layer that, for low curvature 



The indicated preaaura coefficient for the unblown oaae ia alao in 
good agreement with a value obtained from the aurfaca preaaura 
meaauring equipment. 

Figure 136 ahowa the meaaured wall atatic preaaurea around the 
trailing edge for three blcwing oaaea. nieee raaulta were used 
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to provld* the w«ll datum prtaaura for tho radial atatie praaaura 
diatrlbutlona givan In Figuraa 139 •• 142. Tho oonvorganeo of tho 
flow quantltloa with tho wall datum polnta la axeallont. riguroa 6J, 
64 ahow axanqploa of tho raw moaauromonta uaod to doduoo p« tho atatie 
proaauro. Tho offocta of total hoad tuba aiao and ahapo* ahoar flow 
corroction and orrora in volocity moaauramanta aro indicatod, and 
diaouaaod in Saotion 3.4.1. Tho froo atroan total proaauro for thooo 
oxporimonta waa approximataly 0.5 cm of Alcohol bolow atmoaphoric. 

49 

Tho roaulta of Dvorak and Woodward uaing a ainilar twin probo 
tORhniquo to avaluato the radial atatie proaauro diotribution ovor a 
alotcod flap are ahown in Figure 143. Iho aimilaritv in ahapo of tho 
diatributioua ia romarkablOf and parhapo tho flow fiolda could bo 
taken to ba aimilar in nature. 
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5. DIbT‘ bt ION or RBBULTS 

5.1 Ovarall P«rfonMnc« 

5.1.1 Unblown Mrofoll 

The three dimensional lifting performance of the unblown aerofoil, 
Figure 65, was shown to be Iwss than the two dimensional value 
predicted by Hoerner^^ for a similar Reynolds number, as might be 
expected due to the low effective aspect ratio. 



exp 


3.25 per radian 



Hoerner 


4.0 per radian 


The low, indeterminate, aspect ratio of the model clearly affected the 
unblown performance. Due to the end plates the downwash at the trailing 
edge was increased compared to the two-dimensional case, affecting the 
estimation of the lift curve slope. It is interesting to note that the 
effects of the end plates are usually confined to reducing the effective 
incidence of a sharp trailing edged aerofoil where the rear stagnation 
point is fixed. This has the result of reducing the measured lift 
curve slope. The ro'..r ded trailing edge of the current model however 
allows the rear stagnation point to move to produce an opposite, but 
apparently, smaller effect upon the lift curve slope. In general, the 
measured iinblown performance and pressure distributions were 
satisfactory. 

Examination of Figure 66 indicates one farther contribution to 
the lift curve slope discrepancy. The theoretical pressure distribution 
shown was calculated for a true ellipse of identical thickness i chord 
ratio to the model. However the model is not a true ellipsei- 

i) the leading and trailing edges are circular cylinders 


km 1 .: 


oo 



ii) blowing ilott oaus* •urfae* irragularltlai at tha loading and 
trailing odgoa. 

It ahiould bo notod that tha calculation dooa not taka account of 
viacous off acta, l.a. boundary layar growth over tha aurfaca, and 
tharofora cannot be totally rapraaentatlva of the flow over tha latter 
half of the model. However It doea indicate that the Increaae in tha 
leading edge auction peeJc obaerved experimentally waa produced by tha 
circular leading edge compared with a true ellipae. The effect of the 
leading edge slot was not obvious from the measured pressure distri~ 
bution although oil film flow visualisation indicated a small (< 0.1 in 
2.5 mm) separation bubble followed by a turbulent reattachment . The 
lack of asymmetry in the lift generation at positive or negative 
incidences suggests the effect of the sloi lip on C to be small. 

£j 

Figure 67 shows comparisons of the experimental three-dimensional no 
blowing pressure distributions and those obtained from the two- 
dimensional viacoua/potential flow analysis of Dvorak et al^^. These 
results clearly show the presence of the non-elliptic leading edge. 

On the suction surface a high suction peak is indicated in a region 
short of pressure tappings, and on the high pressure surface a similar 
effect is produced. Ibe discrepancies in the pressure distributions 
over the trailing half chord are not fully understood although this 
may be a surface roughness effect. 

No boundary layar trips were included on the model due to the presence 
of the leading edge slot lip and the relatively large teat Reynolds 
number of 1.3 x 10^. If the model surface had loeen ssnoth, natural 
transition would be expected at approximately 35% chord. 

No signs of approaching stall were indicated over the limited incidence 
range (±7.5°) examined. 

measured pressure coefficient at the leading edge stagnation point, 
(see Figure 67) was shown to be within 1% of unity, indicating a 
reasonable evaluation of the bl<’ckage oorreotien. 
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5.1,2 Trallino •do« blowing only 

Th« p*r£or«nanc« o£ th« ••ro£oil ihown in Figur«i 60 - 69 wai vary touch 
ai oxpaetod. Tha a££act o£ tha tip Jata waa aignificant and la 
Indlcatad in Pigura 60. 

Tha advantagaa of operating at lowar valuaa of Cu ara obvioua and lift 
augmentationa OCj^/3Cy) of the order of 60 were attained. Figure 69 
indicatea the effecta of increaaing the alot height: chord ratio and 
the effect of tha presence of a hot wire probe in the wall jet. The 
Inproved performance obtained at the higher slot height was expected 
from the information given in reference 40, which suggests that an 
optimum value for the slot height:chord ratio exists (typically 0.002) 
emd would appear to be a function of the slot design, exit angle 
relative to the local surface and of the relative mass flows between 
the jet and upstream boundary layer flow. 

The presence of the hot wire probe causes a premature local separation 
of tha flow and since the probe was located approximately 1 in (25 mm) 
to the aide of the static pressure tappings, the full effect was not 
measured. In some cases the measured lift was reduced by 20% with the 
probe in place at the trailing edge. It could also be shown that the 
proximity of the probe to the surface affected the measured performance, 
the maximum reduction of being when the probe was closest to the 
surface (see Figure 144) . All the overall performance measurements 
(static pressures) ware ta)cen with the probe and probe support arrange-* 
ment well removed from the static pressure tappings. The question of 
probe Interference with regard to the anemometer results is not so 
serious since the hot wire is a relatively long diet* nee in front of 
the probe holder, and the split film probe, whilst buDciar compared 
with a hot wire, does not require the same physical support (see 
Figure 53 and Sections 3. 2. 5. 3, 3>3.4i5.)> 

The pressura distributions recorded (see Figures 70, 71) exhibit the 
typical • saddle-back' form assoolatsd with circulation control 
aerofoils. 
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rigur* 72 shows ths sstlmatsd downwash corraotlons from tho sffsotivs 

Incldancs calculation (Saction 3.1.1). At posit! vs gaomotric 

32 

inoidanoa tha agraamant with tha rasults of Kind was axoallant. 

Tha changa in slopa for nagativa gaoinatric incidanca is suggastad to 
be due to the reduced affect of tha \ all jet wake upon tha upstream 
flow. Tha more negative the incidanca '-^a closer the wake trajectory 
approaches the tunnel centreline, reducing the under surface wake 
blockage effect. 

Figures 72 and 68 were used to produce Figure 73, the lifting perform- 
ance of the aerofoil for constant effective incidence. The results 

32 38 

are shown compared with Kind and those predicted by CIHCON for the 

32 

current aerofoil. Vftiilst the agreement with Kind is excellent, some 

doubts do exist since the current results were obtained at a higher 

Reynolds number. It is also apparent from Figure 106 that there is a 

large discrepancy in the angle to separation for the various comparable 

flows. In general, the current modal exhibited much lower euigles to 
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separation for similar cu 's than those of Kind or CIRCQM while 

still producing similar lift coefficients. However, the current 

28 

results appeared to fit well with those of Englar , certainly for 
angle to separation. There are several factors which may account for 
these inconsistencies: - 

i) the theoretical calculation procedure was incorrect. 

li) variations of slot design, slot height: chord ratio between 
experimental aerofoils. 

ill) teat Reynolds number variationa. 

iv) dafiniticn of separation point from axpsrlmantal rasults 
a.g. onset of constant surface static pressure, saro shear 
stress or 3U/9y - 0 at wall from velocity profilaa. 

v) surface roughness . 
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Zt appaari that of thaaO) tha two main roaaona for than 
Inconalatanelaa arai* 

1) Firstly tha axparlmantal rasults wara obtalnad for a ralatlvaly 
low slot haightichord ratio (0.000897) conparad with Klnd^^ 
(0.0012) and aa shown in Plgura 69, Incrsasad slot halght 
producad a significant improvement in the lifting performance 
of the aerofoil. 

ii) Secondly, the trailing edge cylinder was not quite smooth and 
this reduced the angle to separation of the Coanda jet and 
hence the circulation for a given Cp. Ihese two effects have 
opposite results on the generation of and hence the dominant 
factor is difficult to determine. Further wor)c to investigate 
the effects of the surface finish of the trailing edge cylinder 
on the current aerofoil may be performed at Bath U:ilversity. 

Figure 74 shows a family of typical performance curves for various 
aerofoil shapes euid indicates the advantages of circulation control 
over the simple jet flap. 

Figure 75 shows the expected increase in with Cp caused by the 
reducticn of surface pressure around the trailing edge and the expected 
Influence of the jet thrust. Since the drag of the aerofoil was not 
the main area of interest in the current study, it wts not fully 
investigated. 

5.1.3 Leading and trailing edge blcwinq 

The effects of leading edge blowing indicated in Figures 76 - 80 show 
a smaller lift decrement than might be expected. The two distinct 
leading edge flow fields are indicated by the pressure distributions 
with and without trailing edge blowing. The boundary between the two 
flow oases was not easily defined since many parameters are involved. 
No hysterisis effect was apparent. It is suggested that the 
parameters involved are those which detesraine the position of the 
leading edge stagnation point relative to the leading edge blowing 
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■lot, naiMlyi- 

I) CM._. This obviously Is ths major paramstsr Involvsd sliirs it 

LE 

dictates ths ability of ths jst to rsmsin attached to ths 
curved surface. 

II) Cw_. Since the trailing edge blowing momentum coefficient 

TE 

controls the trailing edge separation point and hence the 
circulation, it must also have great influence upon the position 
of the leading edge stagnation point. Thus, as was found 
experimentally, the higher the value of the longer the 

flow 'folded-back' for increasing CP^^^ at constant effective 
incidence. 

ill) Effective incidence. The evaluation of the boundary condition 
for the two flow fields was not simple when attempted in a 
closed tunnel since the effective Incidence was a function of 
Cp_. Hence, controlled experiments at constant effective 

XE 

incidence were impractical. However, the variation of 
incidence must have a direct effect since it requires a move- 
ment of the leading edge stagnation point. This in^lies that 
at more positive incidence the leading edge flow would 'fold- 
back' until a higher Cp^^ was reached. 

The technique used to determine the effective incidence of the aero- 
foil (see Section 3.1.1) was unsuitable when leading edge blowing was 
used. The theoretical pressure distribution used in the technique 
was not able to account for the jet flow around the leading edge and 
hence no practical comparison with the experimental results was 
possible. There is as yet no theoretical method capable of solving 
the leading edge jet flow cases. 

The lack of a large number of static tappings in the leading edge 
cylinder and in the very thin slot lip precluded any detailed 
estimations of either the extent or natvure of the separation bubble 
under the fold beck flow or whether the stagnation point wee detached 
from the surface when the jet continued on to the lower surface. 
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5 . 2 V«looity and Turbultno^ Profil«» in th« Uoitraam Boundary 

PlguTM 61, 62, 63 ihow the velocity and longitudinal turbulanca 
prof ilea masurad upatraam of tha blowing a lot. Tha trailing adga 
blowing jat waa axpactad to influanca tha velocity profllaa in two 
waya. 

i) Tha incraaaad circulation and aaaociatad atatic praaaure 

diatribution waa expected to change tha boundary layer growth 
along the upper surface. Thia la clearly r own by the profiles, 
the boundary layer thickness being reduced in the presence of 
the blowing jet. 

il) Close to the slot, the strong entrainment of the boundary layer 
by the blowing jet was expected to change the form of the 
boundary layer velocity profile from the usually assumed 
power law 



( 5 . 1 ) 


This is shown in Figure 83 and the effect of blowing is to increase 
tha velocity gradient close to the wall. The difficulty in 
accurately defining the edge of the boundary layer makes qualitative 
comparison of the velocity profiles unreliable. 

'n\e effect of the leading edge blowing (Figures 61, 62) upon tha 
velocity profiles waa to produce a nett deficit in the circulation of 
the aerofoil and hence a nett decrease in the velocity at the edge of 
the bovndary layer. When the leeding edge jet flow folded back on 
the upper surface, no wake of the jet was observed in the velocity 
profile just upstream of the blowing slot. 

The shape factors shown in Figure 61 are in good agreement considering 
the indeterminate nature of the boundary layer edge and compare well 
with those of Kind^^. 
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Th« longitudinal turbulanca intanaltlaa, (rigura 82) aoroaa 

tha boundary layor wara vary much aa axpaotad, ahowing an inoraaaa 
aa tha wall waa approaohad dua to tha auppraaaion of tha normal 
turbulanca and tha valocity tandlng to zarn. At tha adga of tha 
boundary layar tha valuaa approachad thoaa maaaurad for tha wind 
tunnal frae atream turbulence lavela (aea Section 2.2.2). 

T^a Introduction of the trailing edge blowing jet greatly reduced 
the meaaured turbulence intenaititeSf primarily due to the removal 
of tha onset of separation, the change in the boundary layer growth 
and pressure distribution over the upper surface and the strong 
entrainment of the blowing jet. The introduction of leading edge 
blowing Increased the turbulence intensities compared with the 
trailing edge blowing only case. At the flow conditions Indicated 
in Figure 82, the leading edge jet flow was found to be folded bac)c 
along the upper surface and the increase in the turbulence is ta)cen 
to be a residue of the dlstxurbed upstream flow. 

5 . 3 Trailing Edge Wall Jet 

Of all the detailed trailing edge wall jet measurements, the radial 
static pressure results indicated the most unexpected trends. Prior 
to this work, the radial static pressure distribution had always been 
represented by a simple radial force balance based on a streamline 
curvature, assumed concentric with the local s\urface. This was shown 
to be invalid and was then the area in which the interpretation of the 
results was concentrated. 

5.3.1 Radial static pressure measurements 

The results, shown in Figures 139 - 142, were at first viewed with 
some scepticism, and a thorough examination of the possible sourcea of 
error (Section 3.4.1) was undertaken. However it was shown that tha 
measured distributions indicated valid trends and In scbm oases tha 
removal of possible errors merely accentuated those trends. 
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Figur* 136 ihot't th« radial atatlc praaaura dlatrlbutlon juat up>traam 
of tha blowing >\aot with and without tha blowing jat. Tha agraamant 
i#lth tha uaually aaaunad aaro praoaura gradlant for tha no blowing 
caaa providaa an Mxcallant chaok on tha validity of tha axparinantal 
tadiniqua. Tha maiaurad static praaaura valuas ara in good agraanant 
with tha maasurad airfact praaaura coafficiant valuaa from Flgura 67. 

5. 3.1.1 Physical aicnificanca of raaulta 

Whilst tha maasurad radial static prassura distributions wara not as 
axpactad, thay wara not without some physical justification. 

i) Close to the wall, tha pressure gradients became negative as is 
predicted if the radial momentum equation is allowed to tend to 
y “ 0, giving 

1 3£ . _ 8(v'^) 

P 9y 3y (5.2) 

since v'2 must tend to aero at the wall due to the suppression 
of V' fluctuations by the wall, then 3v'2/3y i® positive. 

Hence, 9p/8y is negative. 

ii) the large positive gradient, in excess of that required for a 
radial force balance, suggeste a reason for strong flow attach- 
ment on a highly curved surface - the Coanda effect. 

iii) the two inner pressure gradients tend to equalise as separation 
is approached. Ibis could )9e suggested aa a separation 
criterion, since a large inner negative gradient would tend to 
reverse the flow curvature. 

iv) the outer negative pressure gradient, which extends approximately 

^ ^min 9 *n«ral less sleep than tha positive 

inner gradient. Biis pressure gradient was thought to influence 
the growth of the shear layer and the mixing rate of the two 

atreaaw . 
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V) At th« adg* of th« flow, th« prottur* gradiont «g«in boooMi 
poAltlvo with a alopa oonaiatant with an ovarall radial feroa 
balanoa. 

For oonpariaon, praaaura diatributiona calmiatad uaing a aimpla 
forca balanca appliad ovar two layara, tha ^at and tha boundary layar, 
ara ahown on Figuraa 139 - 142. nia alopaa wara calculatal aaatming 
a conatant avaraga velocity acroaa each layer from tha maaaurad 
velocity profllaa (Figuraa 64 - 91) and a mean radiua of curvature 
for each layer from tha deduced radii of curvature (Figuraa 101, 103, 
105) . Good agreement between the gradienta in the outer layer ia 
indicated. 

Static preaaure diatributiona of thia form have prevloualy been ahown 

by other reaearchera, but little diacuaaed in their reporta. Aa 

49 

mentioned in Section 4.2.4, Figure 143, Dvorak and Woodward uaed a 
aimilar experimental technique to Inveatigate tha atatic preaaure 
diatribution over a trailing edge Blotted flap, nie reaulta are very 
aimilar to the preaent atudy but are not diacuaaed in their report in 
any detail. The flow through a flap a lot could be ainply repreaented 
aa a low apeed wall jet and conaequently the similarity in thw 
reaulta should not be unexpected. The large physical scale of the 
flap results would also indicate fewer measurement errors due to flow 
angle, wall interference and shear flow, providing further justificat- 
ion for the present results. 

Kruka and Bskinasi^^ also indicated radial static pressure measure- 
ments of similar form for a plane wall jet, but failed to comment on 
their reaulta. 

52 

Millar and Cominga showed aimilar pressure distributions in a 
round free jet flow. Their work uaed a flat plate probe which waa 
optically aligned in tha mean flow direction. The reaulta indicated 
a positive /negative pressure variation in the vicinity of the shear 
layer very similar to that of the present study. 
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26 

EngLar uiad a radially mountad trailing adge flanga plata with 

radial atatlc tappings and gava tha first raaultt for a highly 

curvad trailing cdga wall )at. Intaraatinglyi tha raaulta ahowad a 

nora anvara potitiva praaaura gradlant much cloaar to tha wall than 

would normally ba axpactad. At aoma radial poaltiona and blowing 

conditiona thara waa a tandancy for a nagativa gradlant to appaar 

(aaa Figura 61). Howavar thaae raaulta must ba quaationad and 

perhaps considered unreliable due to plata Interference and corner 

flow interaction in the region of high positive and negative 

28 

longitudinal pressure gradients. Unfortunately Englar published 
no velocity profiles for this study. 


Tha present experimental results for the normal stress (p pv'^) 
were compared with those predicted by Kind^ . Kind used a simple 
numerical method to predict the normal stress based on satisfying a 
radial force balance with two equations to represent tha velocity 
profiles. The normal stress was given as 



suffix a denotes surface 
suffix 6 denotes adge 


where 


P + Pv* 


-P 

2 ^ \af 


(5.3) 


The results are shown as Figure 145 and indicate reasonable agreement 
until tha experimental results change slope in the outer part of tha 
shear layer. This is associated with tha negative static pressure 
gradients obtained in this region. 


5. 3. 1.2 comparison of experimental results with the radial 
momentum equation 

Mhile the measured pressure changes ware small > the small distances 
over which they occurred gave rise to very large pressure gradients. 
Xt was hoped to determine values of ip/iy by direct substitution of 
the measured velocities » radii of curvature and fluctuating 
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compontnti into th« radial momantum aquation* 


U 8v 

R ^ 


+ 


v3v (U^+U'^ -v'^) 

^ - j. 





(u'v' J 


(5.4) 


All of tha tarroa of thla aquation could ba diraotly obtalnad from 
tha axparlaantal raaulta and uaad to provida valuaa for 3p/9y at 
varloua radial poBitiona. Zt waa aaaumad that 


(u*2 _ v'^) << 

and that valuaa for v'^ and u'v' ahould ba corractad aa diacuaaad in 

Saction 4.2.3, Pigurea 132, 133 to taka account of tha low oparational 

bandwidth of tha aplit film probaa. The raaulta obtalnad ara givan 

in Table 18 and Figura 146 and indicate that the equation under- 

aatimatea the meaaured gradienta in aome caaaa by more than an order 

of magnitude. The errors appear to be largest at the limits of the 

y values examined i.e. the edges of the free shear layer. It is 

of interest also to note the magnitude of the normal stress 
■■ '■ ! 

gradient 3v*^/3y in relation to the other terms in the radial equation. 
Clearly if a solution to the wall jet problem uses this equation then 
extreme care must be taken in eliminating or modelling these 
fluctuating terms. 

Since great care had been taken in the evaluation of the experimental 
results, the inability of the radial momentum equation to produce the 
meaaured pressure distributions posed a serious problem, it was 
suspected that tha actual flow field was not adequately described by 
the time averaged results. 
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5. 3,1, 3 An obiTvation ot th« •n>mowf r tlgnaU 

During th« hot wlr« anomonator lurvay of tha trailing adga flow fiald. 

It waa obsarvad that tha high paaa filtarad aignal conaiatantly 

axhibltad a 'blaa* dapandant upon tha poaition of tha proba in tha 

flow, Figura 147 ia an llluatration of tha obaarvad aiqnala, Naar 

U , tha aignal showad a diatinct poaitiva bias which ravaraad naar 
Q 

U Thla affact waa mainly conflnad to maaaui ,r"anta within tha 

min 

firat 30° of tha flow from tha jat. Away from thaaa polnta the 
aignal waa aeen to fluctuate around tha zero voltage point aa would 
normally be expected. Thia observation coupled with the unaatiafactory 
solution of the radial momentum aquation led to a auggeetlon that the 
flow field was not of the simple nature usually assumed. 

5.3.2 Proposed flow field 

It is now proposed that tha actual Coanda jat flow is of tha general 
nature shown in Figure 148. The flow consists of a stream of coherent 
vortices emanating from the slot lip region. The vortices are formed 
from tha excess vorticity produced by the difference in velocities of 
the two flows and the instability of the ensuing vortex sheet. An 
instantaneous flow pattern between the vortices is suggested and 
indicates the mechanism for the strong entrainment associated with 
Coanda flow. The attachment of the wall jet to the surface, Coanda 
effect, may be explained in two ways. 

i) Smith et al^* have shown that in inviscid potential flow, tha 

flow remains attached and spreads due to tha mutually induced 

valocJtiaa of tha vortices. This phenomenon also produces 

aaparation of tha flow at soma distance from tha slot, dapandant 

upon tha initial vortex strength assumed. It is interesting to 

also note tha similarity between tha proposed flow field and tha 

39 

calculation method of Smith at ai 

ii) tha concentration of tha streamlines between each vortex and tha 
surface must produce a local reduction of static prassura hanca, 
each vortex attaches itself. The antzainmant produced by such a 
vortex stream ia sufficiently strong, dapandant t^n tha jat 
strength, to delay separation of the axtamal flow. 
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It li lugg^iUd that it la tha growth of tho diaorata vortloaa dua to 
antrainnantf tha loaa of angular momantum dua to aurfaoa vlaooua 
affaota and tho loaa of total angular momontun to nagotiata tha ourvad 
aurfaoa which dlotato tha axtant of tha attaohmont of tho wall Jot. 

Tha froquancy of tho vortax atraam* tha dacay of tha vortax atrangthr 
tha aitimation of tha initial vortax atrangth and tha affaota of 
gaomatry, axtarnal flow and blowing rata will ba further dlacuaaad in 
relation to the production of an improved theoretical wall jet modal 
in Chapter 6. 

It will be ehown in the following eactione how the time averaged 
reeulta can be deduced from the euggeated inatantaneoue flow flold, 
and ita relation to the meaaured turbulence and preaaure diatributiona 
will alao be ehown. The paaeing frequency of the propoaed flow field 
vorticea waa indeterminate but muat have been above tha upper limit of 
the hot wire anemometer to explain the obaervationa of Section 5. 3. 1.3. 
If it ia aaaumed that the vorticea are convected at approximately tho 
half velocity, \/2 and are initially a alot width in diameter, then a 
typical frequency of 100 KHz ia eaaily attainable. It ia auggeated 
that the vortex centree follow the locua of the half velocity point. 

5. 3. 2.1 Supporting evidence for the piopoaed flow field 

Evidence of the auggeated flow etructure haa been noted by a variety of 
reaearchera although the influence of theae vorticea upon a Coanda flow 
field haa bean neglected. 

54 

Horne and Karamcheti obtained aoma detailed achlieren lAtotographa of 
thia type of flow in a plane laminar jet in the preaence of a variable 
length wall. The photographa were coupled with hot wire probe outputs 
which clearly indicated the periodic nature of the flow. The vorticea 
were shown to exist even in the preaence of a curved wall, although 
the growth rate of the vorticea waa far greater than that for a plane 
wall. 
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55 

Daviaa and Baxtar dlicuai tha davalopmant of ring vortlcaa in a 
fraa air jat from tha Kalvin-Halnholts inatabiJlty of tha vort«c 
shaat and tha Importanca of the vertical upon tha traniltion of tha 
fras ahaar flow. Hiair raaulta indioatad a strong Raynoldi numbar 
dependanca and tha flow appaarad to ba wall ordarad balow a Raynolda 
numbar of ISOOO basad on tha jat diamatar. At Raynolda numbara abova 
20000, the onaat of chaotic motion is incraaaad. A hydrogen bubble 
technique was used to visualise the flow field. 

Damms and Kuchemann^® proposed a similar flow field to predict the 
mixing region between two parallel streams of different velocities in 
the wake of a splitter plate. Their vortices were allowed to have 
cores and the effects of either equidistant or exponential spacings of 
the cores were Investigated. The growth rates of the cores were 
investigated, and time averaged velocities across a shear layer 
produced. The work was based upon the experimental results of Brown 
and Roshko . They used a shadowgraph technique to trace the vortex 
cores formed behind a splitter plate between two streams of different 
density (nitrogen and air) in a pressurised working section 6 atmos) . 
A vortex streaming frequency of approximately 10 KHz was found at a 
streaming speed of approximately 11 ms“^. The vortex cores were 
apparently formed as the vortex sheet between the fluids rolled up at 
regular intervals. The experiments were run at relatively high 
Reynolds numbers. Other experimental evidence is given in reference 
56, in particular, Clark and Markland^® found evidence of a core type 
structure in the free shear layer of a plane wall jet. 

Further evidence and discussion of the cause, structure and persistence 

of large vortical structures is given in references 59-62. in 

59 

particular, Wygnanski et al show photographic evidence of the 

62 

structure and Chandrsuda et al show the process of helical pairing 
which will be further discussed in Section 6.2.3. 

The proposed flow field is therefore well supported by previous 
research . 


Th« ipatlal corr«l«tion of th« propoicd vortox itraam and tha 
axparimantal valocity profilaa, turbulanca quantitiaa and radial 
atatie praasura diatrlbutlona for a alngla angular poaition and 
blowing rata ia given in Figure 149. A full interpretation of each 
of tha experimental parametera with reapect to the vortex atream ia 
given in Sections 5.4.2r 5.4. 3^ 5.4.4. It is immediately obvious 
from Figure 149 how the vortex stream relates to the time averaged 
velocity profiles and also in particular how the stream allows the 
radial static pressure to depart from the normally smooth distri- 
butions within the bounds of the free shear layer. 

5.4.1 Overall performance 

The proposed flow field now provides a simple physica] explanation of 
the Coanda attachment* its dependence upon the jet blowing rate and a 
number of the variations of overall performance previously Indicated 
(Section 5.1). 

5 . 4 . 1 . 1 Effect of increasing slot height 

For identical Cy values an increased slot height produces a lower jet 
velocity; this would tend to suggest a weaker vortex strength at the 
initiation of the vortex stream. Also* the streaming spaed will be 
reduced. In general* it would be expected that an increase in slot 
height wonld reduce the lifting performance of a given aerofoil. This 
is shown to be true in reference 40 for slot height* chord ratios of 
9Jf**ter than 0.0166. il>e current tests ware conducted at ratios lower 
than this figure (generally 0.000897) and hence it would appear that 
an optimum ratio would exist dependent upon the relative width of the 
J*t* the boundary layer within the jet and the \ipstream boundary layer* 
below which the trend was reversed. Section 5. 4. 2.1 contains a further 
discussion on this point. 


5. 4. 1.2 Effect of probe interference 

Nhilst it ia not suggested that the individual vortices are of full 
spanwiae extent* the shedding frequency should be high enough to avoid 
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dltcontlnultlcs In attachmant# nia affact howavar of introducing an 
obstruction into tha Coanda flow (an anamonatar proba for axjnpla) 
would ba to ruptura tha vortax ihaat and causa a 'roll-up' of tha 
vorticas eithar sida of tha obstruction. Tha rasultlng longitudinal 
vorticas ara a wall known phanomanoro, as is the tendency of the jat 
sheet to split in the wake of an obstruction. 

5.4.2 Velocity profiles 

The time averaged velocity profiles U, as measured by a split film 
probe, (Figures 84 - 91) were much as expected. It is obvious trom 
Figure 148, how the proposed flow field could be used to produce 
these time averaged results. Further discussion of this point will 
be given in Section 6.2.1 in relation to determining the true vortex 
strength from a given time averaged velocity profile. 

The U results show excellent agreement with the one comparable profile 
32 

of Kind , (Figure 98) and the other two hot wire techniques used in 
this study, (Figure 99) . A discussion of the profile errors of the 
anemometer equipment is given in Section 3.2.5. 

5 . 4 . 2 . 1 Effect of slot height variations 

The effect of Increasing the slot height upon the measured velocity 
profiles (Figure 106) is seen to be primarily a change in jet velocity 
profile. The existence of an approximately constant velocity core is 
apparent although outside the velocity minimum, little difference 
exists in the profile. Unfortunately time did not permit further 
downstream investigation of the wall jet at the increased slot hel^t. 
With reference to the effects of increasing slot height upon overall 
performance discussed in Section 5. 4. 1.1, it is suggested that the 
optimvun performance is Influenced by the bovndary layers in the jet. 

Ihe presence of the boundary layers in the jet produces a concentration 
of vorticity within those boundary layers proportional to the velocity 
gradient. The presence of a constant velocity core provides the 
optimum conditions for this vorticity concentration, i.e. maximta 
au/ay. Therefore since the magnitude of the vorticity dafines the 
distance to separation of the wall jet, it ia proposed that it is the 
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■lz« and gradiant of tha jat boundary layar (dapandant upon tha 
axiatanca of a cora, or not), couplad with tha af facts of tha 
proximity of tha vortax atraam to tha wall which produces tha 
variation of parformanca of circulation control aerofoils . 

5.4. 2. 2 Effects of slot liP thickness 

Tha effect of increasing the thickness of the slot lip, (Figure 107) 
is very interesting. It was expected that the effect would be to 
lengthen the starting region of the flow in the wake of the lip and 
reduce the lift coefficient. However, this was not apparent from the 
results, which implied that the entrainment process was strong enough 
to produce inflow angles of up to 25°. This also suggests that the 
formation of the proposed vortex stream does not require a starting 
length dependent upon geometry, but rather on velocity ratio. 

5.4. 2. 3 Variations of the parameters which describe the velocity 
profiles — 

Figure 109 shows the variation of the position of the half velocity 

point, y with 0 for a variety of blowing rates. The half velocity 
^ ^a/2 

point may be defined two ways; 




U + 0 , 
m min 


where 

- 

2 



2 


(5.5) 



0 + u , 
m 0 

2 


or 

\ - 



2 




Figure 109 gives both variations. However it is shown (Figure 109) 
that at higher blowing rates where a fully developed wall Jet profile 
is attained, the two half velocity points coincide 0^) . 

Figure 110 shows the variation of the wall Jet profile edge velocity 
with 9 and Cp . In sons cases, detemination of G was difficult due 
to the variety of wall Jet velocity profiles obtained. 
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Pigur* 111 ihowc th« variation of with 6 and Cu. 0 ii 
aaaumad to tand to aero at tha alot. Pigura 112 .howi thi"variation 
of with e and Cu. Thaaa two figuraa giva soma indication of tha 
variation of vortex strength in tha wall jet. it could be aaaumad, 
allowing for the fact that tha velocity profilaa are time averaged, 
that the vortex strength ia given by 

K‘ W- '>'m- W] 


Where G is a factor to allow for the effects of time averaging upon 
the individual vortex strengths, (y^ - y^^^) represents a typical 
length scale for the vortices, assumed to be their diameter, and is 
equal to the width of the free shear layer. The velocity difference 
across the shear layer (U„ - , coupled with G, Indicates the 

velocity increment across the vortices. Figure 150 shows the 
variation of with 0 and Cp which can be seen to be an 

exponential type decay, ihe method of Smith et al^® uses a similar 
type decay to model the vortex strength distribution along the wall 

jet. Some comments upon the relative size of G will be made in 
Chapter 6. 

Streamline patterns and their radii of curvature 

The time averaged results, U, shown in Pigure. 84 - 91 were used to 
produce the streamline patterns and their radii of curvature 
(Pigures lOO - 105) . the effect, of the strong entrainment close to 
the .lot and the growth rate of the jet are obvious from the stream- 
line patterns. The radii of curvature were produced, a. described in 
Section 4.2.2 and show a distinct departure from the usually assumed 
concentricity with the local centre of surface curvature « 


Equation (5.7) is shown on the figures for comparison and it can be 
aeen that at higher blowing rates and angles from the slot, tha 
variation of the deduced radius of curvature becomes parallel to 


•quatlon (5.7). This would appear to indicate the occurence of a 
fully developed wall jet profile once clear of the slot and the 
aeeociated strong entrainment. The deduced radiue of curvature can 
be aa much as twice that usually assumed by equation (5.7). This 
would produce significant errors in many theories particularly where 
the radial static pressure gradient is given by a simple force 
balance (equation 3.33). 

The measured local flow angles V/U, (Figures 92 - 95) seemed to 
indicate the expected trends of entrainment and jet growth. However 
the results when compared with the requirements of continuity, 
(Figure 97) 


^ 3V 

3y " (5.8) 

were shown to be too negative. It was suspected that the errors 
discussed in Section 3.3.4. produced this apparent negative offset 
and allowing for this, the results provide a reasonable representation 
of the streamline patterns shown in Figures 100, 102, 104. 

The derivation of the normal velocity profiles from equation (5.8) is 
suggested as realistic since the vortices in the proposed flow field, 
although markedly different from the time averaged streamline patterns, 
would produce no nett effect upon V; as each vortex passes, then 
equal and opposite V influences are experienced producing a nett aero 
effect for a constant radial position above the surface. 

5.4.3 Turbulence results 

The first consideration regarding the measured turbulence parameters 
in the light of the proposed flow field is, exactly what has the 
equipment measured? Ihe instantaneous velocity vector instead of 
being simply represented thus 
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could now be 


U - U + Q + u' 


(5.9) 
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consiitant with th« suppresilon of th« normal turbulanca at tha wall^ 
V' ^ * 0 aa y 0 

The effects of heat loss to the surface from the hot wire probe did 

not seem to be a problem, although the closest points were within 

the region (y < 0.004 in, 0.1 mm) generally regarded with suspicion 

for wire probes of 5 gm diameter. The effects of the wall upon the 

split film probe results are however more noticeable. Whilst heat 

loss to the surface may have affected the mean velocity measurements, 

the effects of probe interference are more evident in the turbulence 

measurements. This is because the turbulence parameters are 

calculated as the fluctuations upon the mean, which is calculated at 

each step) hence the turbulence levels are relatively insensitive to 

changes in mean flow angle or heat loss to the surface. All split 

film probe turbulence measurements within 6 probe diameters 

(~ 0.035 in, 0.88 mm) of the surface have been ignored. It should be 

remembered that the film probes were run at a much lower overheat 

ratio (O.S) compared with the wire probes (0.8). The effect of the 

wall interference on the split film probe is clearly shown in 

Figure 117 as an increase in Iseyond that expected as the wall 

is approached. Figure 117 also shows a comparison Isetween the 
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present split film results and the hot wire results of Kind . The 

agreement is remaricable and not only provides justification for the 

assumptions regarding the extent of probe interference but also 

32 

indicative that the results of Kind were ta)(en at a similarly low 
operational bandwidth. 

In the region of the shear layer, Figures 113 - 116, the turbulence 
results must now be reinterpreted for the presence of the vortex 
stream. The time averaged longitudinal intensity (/u^^/5) indicates 
a peak value in the outer half of the shear layer. This is explained 
simply as the region In which the longitudinal velocity vector 
experiences the greatest variations. Dependent upon the relative 
position of a vortex, the velocity vector can be rotated through a 
full 180^ as Indicated In Figure 151. 
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Plgur** 119, 120 mor* cltarly illuitr«t« th« r«latlv« poaitlon* of 
tha naxinun and mlnlinuiD longitudinal tufbulanca intanaitioa to tha 
valocity maximun and minimum. 

The affacta of tha excaaaiva surface roughneaa raantioned in Section 
5.1.2, upon the turbulence raaulte is not obvious since no suitable 
data exists to provide a comparison. 

The longitudinal turbulence intensities, (Figures 113 - 116) exhibit 
a minimum time averaged profile a short distance downstream of the 
slot, in common with the other turbulence parameters. This is 
thought to be an indication of the change in sign of the longitudinal 
pressure gradient. The negative pressure gradient suppresses the 
turbulence while the positive gradient has the opposite effect. The 
e position of this minimum in the turbulence profiles appears to be 
nearly constant for varying Cy emd was consistent with the position 
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of the measured suction peak. A similar effect was noticed by Englar 
on measurements of surface shear stress, an example of which is given 
as Figure 151. 

As separation is approached, Figures 113 - 116 Indicate that the 
intensity of the turbulence increases rapidly with distance from the 
slot, as would be expected. 

Figure 118 shows the longitudinal turbulence intensity as measured by 
a hot wire probe at similar Cy values but for variations of slot 
height and lip thickness. Increasing the slot height reduces the 
minimus intensity due to the presence of a near constant velocity 
core. However, the maximum Intensity near the slot lip has been 
nearly doubled. The doubled slot height effectively reduces the jet 
velocity for a given Cy by nearly 30%, thus reducing the initial 
mixing procees and the formation of the vortex stream, this results 
in a higher nett turbulence intensity in the wake of the lip. 

Ihickening the lip shows an increase in the value of the peak intensity 
in the shear layer and a corresponding expansion of the wake of the 
lip. The effects upon the velocity profiles, (Section 5. 4.2. 2, 
figure 107) were shown to be small and this indicates that the 
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■tarting langth was littla changad by tha lip thicknaaa. Tha 
affaota of thickaning tha lip upon tha ovarall lift coafflclanta 
wara within tha axparlmantal accuracy of tha maaauramant ayatam and 
hanca could not ba uaad to Indlcata any tranda. 

In both doubling the slot heights and thickening the slot lip, tha 
typical eddy size of the flow In the wake of the lip has been 
Increased by reducing the velocity gradients and increasing the 
typical lengths. This could also be a factor in the apparent increase 
in the measured turbulence in the wake of the lip. 

5.4. 3.2. The normal turbulence intensity 

The normcl turbulence distributions are shown in Figures 121 - 124. 

The meaaurvnents for this parameter ware only obtained using the 
split film probe and as such are subject to the operational limitations 
discussed in Sections 3.3.3, 3.3.4. The values close to the wall 
have had to be ignored due to probe interference effects as discussed 
in Section 5. 4. 3.1. It is interesting to note that as the wall was 
approached, the lower film would cause a small contraction between 
the probe and the surface and therefore reduce the fluctuations 
measured by the lower film. Now, since the normal fluctuations were 
a function of the difference of the two film signals, this asymnetry 
would cause an apparent increase in the measured fluctuations. This 
can be seen to be the case as shown at one blowing rate in Figure 124. 
Zn general the values for the nomal turbulence were approximately 
60% of the longitudinal values as shown in Figure 125. 

Ihe normal turbulence intensity maximum appeared to be located at the 
centre of the shear layer, this being the proposed centreline of the 
vortex stream. This result seams reasonable since, at tha vortex 
centreline, as each vortex passes, the velocity vector rotates 
through 160° to give the maximum variation of V, (see Figure 151). 

Figure 126 shows that, at lower blowing rates, the positions of the 
maxima of the normal and longitudinal turbulenoe Intensities do not 
eolnelde, the normal turbulenoe being closer to the surfaoe. This 
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relate* well to the propoeed vortex atream. However » at higher 
blowing rat*a« the longitudinal turbulence naxlmum move* oloaer to 
the eurfao* to colnold* with the nonnal maximum. Thla ph*n<«n*non la 
not fully underatood but may be due to the vortlcea not being 
circular. A* the jet velocity la Increaaed, ao the vortlcea are 
deformed allowing the turbulence Intenalty maxima to converge. 

5 . 4 . 3 . 3 The Reynold* ahear atreae u*v * 

Figures 127 *■ 130 show the shear stress distributions measured by a 

split film probe for various angular positions and blowing rates. | 

The results r although subject to the errors discussed in Section 

3.3.4, 5.4.3. 2 indicated the expected trends. A negative value near 

the wall was not measured but was indicated by the results, in agree* 

ment with a positive wall shear | 

T ■ - p u*v' I 

f 

28 ' 

An example of the wall shear stress as measured by Englar is given 

in Figure 152. ) 

Just inboeurd of the velocity maximum a zero shear stress point was 
indicated although it was apparently at a greater y/y value than was 
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suggested for a curved wall jet by Wilson and Goldstein This may 
be an effect of the surface roughness of the trailing edge cylinder. 


A maximum value of the shear stress occurred close to y^^^ 
distributions then quickly reduced to a second zero value at 
approximately Outboard of this zero point the shear stres* 

became negative 1 consistent with the return of a poaitive velocity 
gradient. 


Figures 134, 13S give a conparieon between the measured dletri but lone, 
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those produced by CZRCON and those steaeured by Wilson and Goldstein 

They indicate the large difference in the distribution caused by the 

presence of the external flow and also tha inability of the eddy 

visooeity model used in CZAOGM^^ to predict tha shear straaa in the 
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out«r raglon of tha flow. Flquraa 134, 135 uaa tha raaulta ahown In 
rigura 133 for tha naaaurad u'v' faotorad to taka account of tha 
llmltad bandwidth of tha split film probaa, aa dlaouaaad In faction 
4.2.3, 


Tha Intarpratatiorv of thaaa shaar atraaa raaulta with raqard to tha 
propoaad vortax straam is mors oomplax than for tha normal atraasaa 
u'^, v'2. It would Beam raaaonabla howevar to axpact tha maximum 
Bhaar to occur in the vicinity of tha counter flowing atraam between 
each vortex (see Figure 151). Also, since tha fluctuating velocity 
at tha inner edge of the vortex stream imparts a shear force on the 
surrounding fluid of constant sign, then it is reasonable for the time 
averaged shaar streaa at that point to be finite and non-sero. Hence 
the zero shear stress may wall occur away from the velocity maximum. 

At the outer edge of the vortax straam however, the shear stress is of 
fluctuating sign. This suggests that the time averaged shear stress 
should be zeroj as was shown to be the case (Figure 127 - 130) . 

It la necessary to remember that the measured shear stresses (u^) 
are perhaps not a realistic representation of the actual shears 
experienced by the fluid since the measurements were time averaged with 
respect to a fixed co-ordinate system. 

5*4*3.4 The radiaJ static pressure distribution 

These results have already been discussed in Section 5.3.1, Thm 
unusual shape of the static pressure distributions is however now simply 
explained by the proposed discrete vortex flow field. The high 
pressure gradients that were measured represented the time averaged 
values of the pressure gradients appropriate to the high curvature of 
the streamlines of the flow around their local centres of curvature for 
each passing vortex, combined with the spacing of the vortices. 

The radii of the proposed streamlines around tha vortices is roughly 
an order of magnitude smaller than the surface radius (see Figure 141) . 
This explains the disageemsnt between the measured pressure gradients 
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and thoM produoad by tha radial monantum aquation. 

5.5 Ganaral Conoluaione Raqardinq tha Bxparimantal Raaulta 

Qsaat cara waa takan in obtaining and avaluating tha axparimantal 
raaulta. Tha unaxpactad tranda of tha radial atatic praaaura 
diatributiona lad to a naw undaratanding of a Coanda flow fiald which 
not only allowad a full intarpretation of tha turbulanca maaaurotnanta , 
but alao axplainad many of tha unuaual phanonana aaaociatad with tha 
parformanca of circulation control aarofoila. 

It ia hopad that tha propoaed flow fiald will lead to an improved 
theoretical prediction and the incorporation of auch a discrete 
vortex calculatio'i into an aerofoil prediction schme will now be 
diacuaaed in some detail in Chapter 6. 
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6. THE THEOMBTICAI. IMPLICATIONS OP THB PRB8BNT RESULTS AND THE 
PROPOSED PLOW FIELD 

6-1 T h« Dsfloitncap of th« M«thod of Dvorak 

Th* dlacustion of th# currant raiulti in Chaptar 5 Indioatad a numbar 
of daflclanclaa, othar than tha aaaunad flow fialdi in tha tlma 
avaragad paramatara uaad in tha trailing adga wall jat calculation of 
Dvorak and Kind^®, 


i) Radial ahaar atraaa diatribution 

Figura 134 givaa a jompariaon batwaan a typical ahaar atraaa diatri- 
bution predictad by Dvorak and Kind and the currant raaulta. Tha 
eddy viacoaity model uaad by Dvorak and Kind^® aquation (1.4) ia uaad 
to calculate the ahaar atraaa diatribution and it can clearly be aaan 
how the model fails to predict the negative values in the outer region. 
This negative shear stress region was also shown to exist by Jonas^^ 
but only close to the slot exit. 


11) Radial static preaaure gradient 

The meaaured static pressure gradients (Figures 139 - 142) indicate a 

significant departure from the usually assumed distribution given by an 

overall radial force balance. The values used in CIRCON^® are taken 

from tha potential flow calculation with a correction for excess jat 

momentum. Only tha longitudinal pressure gradient is calculated by 
36 ^ 

CIRCON and an example of the distribution is shown in Figure 153. 


ill) local streamline radius of curvature 

The experimental results (Figures 101, 103, 105) indicated that tha 
local streamline radius of curvature was considerably larger than tha 
usually assunad value of(R^ + y) . This would have an effect throughout 
tha wall jat calculation, most obviouely in tha continuity aquation. 
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and in tha calculation of the following tozina in the radial wonenturo 
aquation. 


£ iV 1 3(u*v» ) 

R 30 ' R ' R ' 36 


It should of course be noted that the use of such time averaged 
relationships has produced a calculation method that is not based on 
a true representation of the flow field. The above deficiencies exist 
'^^thin the present calculation scheme; however there are seme other 
aspects of the assumed flow field that CIRCON^® is unable to account 
for, such as; 

i> jet exit angle relative to the local surface 
ii) surface roughness 

Hi) variation of jet contraction geometry 
Iv) variation of slot lip geometry 

It is suggested that all of these may produce large changes in the 
performance of a given aerofoil. 


If a wall jet calculation scheme is to be based upon the solution of 
the momentum equation, closed by a suitable turbulence model, then the 
finite difference technique is considered to be the most suitable. 
However, if the proposed flow field is to be accurately modelled then 
an inordinately fine mesh grid would be necessary, creating problems 
in terms of computer time and storage. It is proposed that a more 
physically representative approach would be to incorporate a discrete 
vortex wall jet model into the overall viscous/potential flow inter- 
action technique. 


114 


6.2 Son* Conntnts upon th« Incorporation of a Diierata Vortex Modal 
Into an Aerofoil Calculation Schama 

Vlacoua/potentlal flow interaction tachnlquaa exist which enable the 
calculation of the boundary conditions around the trailing edge wall 
jet of a circulation control aerofoil. These boundary conditions are 
very important in this type of calculation since the conditions at the 
slot determine the initial vorticity and the lower surface boundary 
conditions usually define the wall jet separation pressure. The 
present results may be used to make some suggestions regarding the 
improvement of the existing discrete vortex wall jet calculation 
scheme (Smith et al^^' . 

6.2.1 Determination of the initi-^.l vortex strength 

Since it is known that apart from viscous effects the vorticity cannot 
be destroyed in a flow, merely dissipated, the determination of the 
initial vortex strength becomes of prime importance. It is this 
parameter that will control the extent of the Coanda attachment for a 
given surface geometry. 

It would seem reasonable that the jet; free stream velocity difference 
should be the major factor in determining the vorticity? the external 
and slot flow are assumed steady. A realistic value should be obtained 
by considering the velocity profiles and in particular the instantaneous 
value of (Ujjj - close to the slot. It can be seen from Figure 149, 

that the time averaged values of this velocity difference exhibit am 
exponential type decay, suggesting some similarity between this 
parameter and the proposed vortex decay of Smith et al in terms of 
energy dissipation. Before considering the actual vortex strength 
within a Coanda jet, two other factors must be considered? the frequency 
of the vortex shedding and the starting length of the flow. The 
frequency is considered in Section 6.2.2. 

iha question of starting length has long been a problem in all wall jet 
oalculaticxi methods, although problems are usually centred around the 
poaalble existence of a potential core and the extent of the slot lip 
wake, in the case of the log spiral trailing edge it is usual to assume 
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* origin to tho Coanda lurfaco In ordar to produoa tha raqulrad 

ilmilarlty. Aa shovm in rigura 146 tha propoaad flow fiald ia nora 
complax whan conaidaring atarting langth. 

From tha axparimental reaulca (Pigura 107) tha extant of the alot lip 
wa)ce appeared to be small, of tne order of 5 alot widths, and did not 
vary significantly with changes in alot lip thickness. This would 
suggest that a value for the wake length as a function of the blowing 
rate would be reasonable and not significantly in error, it is proposed 
that the wake length would reduce with increasing blowing due to the 
increased entrainment of the stronger jet. 

The more difficult problem is that of the transition from the assumed 
vortex sheet at the end of the lip wake to the fully developed vortex 
flow. More recently Jimemea^^ has considered this problem without the 
presence of the wall. Considering an Infinite uniform vortex sheet 
with strength ^U, then if the sheet is perturbed locally it will tend 
to roll up into a tight spiral vhose evolution and growth rate can be 
determined. 

A double armed spiral is formed with a tight almost circular, central 
core. The streamlines form the familiar 'cats-eye' type pattern 
(Fi^^^^® 154) and indicate that the la^ge vortical structures are 
elliptical in shape with a semi axis ratio of 1.76. The growth rate 
® free shear layer cc*^structed of a number of such structures has 
been shown by Jimenez^ ^ to be: 

X " (Uj^ + Oj) 

Ibis is in good agreement with the experimental observation* of Brown 
and Roshko^^. 


The distance over which the vortex sheet rolls up to form such discrete 
structures is not clear. Bowever it would be logical to suggest that, 
providing the perturbation was not a randomly occurring event but 
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contlnuou*, th.n th. .trarting l.ngth could ba approximated to on# 

wavelength of the vortex etream, eimllar to the photographic evidence 

62 

of Chandreuda et al 

The required perturbation could be due to either three dlmeneional 
effects within the turbulence structure or, as ie more likely in the 
case of the wall jet, to the strong entrainment of the outer flow at 

the Jet exit. 

It is clear th* these suggestions are reasonable for free shear layers. 
The effects of the presence of the wall and viscosity upon the formation 
of the structures is as yet undetermined. The wall will presumably 
have a more significait effect upon the growth ra^e than on the 
formation of the structures and this will be discus 1 in Section 6.2.3. 


Having determined some criteria for the starting length, the problem of 
relating the time averaged velocity profiles to the instantaneous 
vortex strengths is now considered. Assuming the simplified vortex 
array and with the notation as shown in Figure 155 it can be shown 
that, assuming the vartlces to exactly span the shear layer and that 
locally, the vortex strengths and size are constant, then 
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(6.3) 


If this ie Integrated to produce the time averaged velocity then: 

“x.r ■¥■ L r * 7 


-n 


Hence the ratio of instantaneous to time averaged velocity at the edge 
of the si^lifled vortex stream is given by, 
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Thii li taken at a point on the edge of the shear layer coincident 
with a vortex centreline. The parameter d is the wavelength of the 
stream, Figure 1S4 shows some results from this calculation for 
ratios of vortex strength for typical length scales and indicates the 
magnitude of the corrections in order to allow for time averaged flow. 

Using the simple model of Figure 155 it is obvious that the 
instantaneous value of V, the normal velocity at any vortex centre, 
providing it is not at the extremeties of the stream, will always be 
zero, since the normal components of the vortices tend to cancel. 

Thus it becomes clear that the dissipation of the vortex strength and 
changes in vortex spacing dictate the locus of the vortex stream. 

This will be further discussed in Section 6.2.3. 


6.2.2 Determination of the vortex shedding frequency 

The determination of the vortex shedding frequency is of vital 

importance to any discrete vortex theoretical method since it dictates 

the time step interval of the calculation. Ttie current experimental 

work suggested that the frequencies of the stream were above the 

normal operating bandwidth (20 - 40 KHz) of the hot wire anemometers. 

This would tend to be supported by other experimental evidence, such 
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as Brown and Roshko , Horne et al . it is also reasonable to 
suggest that, not including vortex pairing or dissipation, which will 
be discussed in the next section, the time interval between vortices 
remains constant for the wall jet. Now, since 


fr«,u.ncy - 


Streaming speed 
separation distance 


( 6 . 6 ) 


the distance between adjacent vortices must vary with changes in 
streaming speed as well as with their mutual interaction. Figure 157 
shows the variation of «m/2' the proposed streaming speed, with angle 
from the jet exit. Since is measured at the vortex centres, its 

variation will be alsmst independent of the else, strength and 
frequency of the vortex stream. 
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Ths shedding frequency is usually determined, such as In the casts of 
cylinders, blunt aerofoils and spoilers, from the Strouhal number. 


where d and U are assumed typical lengths and velocities (such as the 
free stream velocity and the dieuneter in the case of the cylinder) . 

For a given situation, the Strouhal number remains constant and the 
frequency then becomes a linear variable with velocity. 

For the case of the I *.qhly curved wall jet with an external stream, it 
is proposed that, 

d > f^ (slot lip thickness, jet boundary layer displacement 
thickness) 


U “ velocity difference at the slot 


The effect of slot height upon the value of d, given as 6 , the 

displacement thickness of the jet boundary layer, Is suggested in 
order to account for the presence of the jet boundary layers, as these 
are the source of the majority of the vorticity flux (-3u/3y) , and the 
possible existence of a constant velocity core. It is possible that a 
similar effect for the upstream boundary layers should be included. 
However the strong entrainment at the slot exit has been shown to alter 
the shape of the upstream boundary layer velocity profile and 
insufficient data is available on which to base any qualitative 
suggestions. Certainly, since the vorticity of a fluid element is 
defined as 


It) 


av 
“ "ST 


du 


( 6 . 8 ) 


the largest contribution will coma from the jet boundary layer/slot lip 
wake interaction i.nd the suggested scaling length is representative, 
rrea the current geometry where d = 0.25 x 10*’ m (0.010 inches) and 
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Msuming a conatant Strouhfcl nunbar of 0.2 (typical of vortax ahaddlng 
flows), tha fraquanoy variation shown in Flgura 158 can ba obtainad. 
Ttia fraquancias indicatad ara in tha corract ranja comparad with tha 
bandwidth of tha anenomatar aquipment and suggast a wavalangth 
(distance between adjacent vortices) given by 

U 

^ ^ _8 ^ streaming speed 
f frequency 

of the order of 1 2 mm (0.04 -► 0.08 inches) . 


Since the formation of the vortices is suggested to be largely 
dependent upon the interaction of the jet boundary layer and the slot 

lip wake, the effects of the proximity of the wall are expected to be 
small. 

On the gr owth and dissipation of the vortex stream 

The vortex stream which makes up the shear layer of the wall jet grows 
due to the individual growth of each vortex by entrainment. As the 
vortices roll up, they continue to entrain the outer fluid at a rate 
dependent upon the individual vortex strengths and the spacing between 
adjacent vortices. In his analysis, Jimenez^ ^ considered the entrain- 
ment process of each vortex. He showed that the vortices will continue 
to grow in the familiar 'cats-eye' shape but that the expansion of the 
individual vortex in the lateral direction saturates before the 
longitudinal axis. The limiting condition then becomes the point at 
which adjacent vortices touch. At that point it is suggested that the 
process of pairing occurs. Op to the present time this process has not 
been allowed for in discrete vortex wall jet calculations. 

The pairing process is a very complex cne. There are two possibilities 
to bo considered. The first possibility is that, at the very point 
that the vortices touch, the effect of having opposite velocity vectors 
at the same point causes a severe disruption of the local streamline 
jwttern and instant amalgamation of the two vortices occurs. Upon 
smalgamatlon, the nett vortlclty will be conserved while the wavelength 
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of th« local vortex stream is doubled. If the cross sectional area 
of the 'new' vortex is roughly double that of a previous vortex, then 
separation between vortices again exists and entrainment continues. 

The second posoibility usually considered is that, upon touching, the 
two vortices roll round each other for a fraction of a turn before 
amalgamatlonj the so called 'helical pairing'. This process is 
indicated by the results shown by Chandrsuda et al®^. It is interest- 
ing to note that Jimenez®^ also shows that the pairing and entrainment 
process can approach a limit cycle at which point the vortices are no 
longer generated downstream. This is interesting and could have 
significance in determining the extent of the wall jet attachment 
for a given blowing momentum. 

The effect of the presence of the wall is obviously important in 
determining the growth rate of the vortex stream. Compared with a 
free shear layer, the growth of the vortices is restricted due to the 
proximity of the solid surface and of course a certain amount of 
energy will be lost from the stream due to viscosity within the inner 
layer. The growth of a free shear layer derived by Jimenez®^ is given 
as equation (6.3) and indicated in Figure 159. The reduction due to 
the presence of the wall and the effect of the conservation of 
angular momentum in the curved wall case are clear. 

The single vortex stream can be simply modelled in potential flow and 
a wall can be included by adding the usual mirror images. Introduction 
of the additional boundary condition that u->0 as y -► o is more complex 
and is yet to be satisfactorily attempted. 

It is generally assumed that the vorticity will be dissipated by the 
smaller scale eddies within a flow and that large vortical structures 
will degrade to a smaller scale - the so called 'energy-cascade', 
although this is usually a strong three dimensional flow, in the 
proposed flow field it is obvious that no significant degradation of 
the flow occurs, it is therefore suggested that the energy is 
dissipated in three waysi- 
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I) by tha auall acala turbulanoa aurroundlng aach Individual 
vortax 

II) by the catastrophic pairing phenomenon 

ill) via the wall ^hear atreas, which exerts a moment on the 
vortices . 

It Immediately becomes obvious that the modelling of such a flow field 

is inherently extremely complicatedf and the estimation of the 

dissipation of the vortical energy is a key part of the technique. 
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The work of Jimenez would seem to indicate that the pairing process 
is the dominant factor in the growth of the shear layer and that 
inclusion of a limit on the size of adjacent vortices is necessary. 

Once the limit is reached then pairing should be accounted for with a 
coincident loss in vorticity for the new single structure. 

The assumption that the structures can be represented by simple point 
vortices is valid since, although the deduced theoretical streamline 
patterns tend to cylindrical cores, when viewed from a distance, a 
vortex with a cylindrical core will appear to be a simple point vortex. 
With this in mind it would suggest that the proposed pairing limit 
should be checked by estimating when two adjacent 'cats-eye' stream- 
lines touch. 
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The invlscid potential flow method of Smith et al , has already shown 
that a discrete vortex stream will remain attached to a highly curved 
surface. It was however necessary to include an exponential vortex 
strength decay to obtain reasonable agreement with experimentc.1 lift 
coefficients. It la hoped that incorporation of some of the above 
suggestions may allow ren»val of the somewhat arbitrary vortex decay 
rate without reducing the accuracy. 
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Sow further oonBldTatlona 

Th. .bov. oo™„t. h.v. b..n conc.rn.4 prl».rily with pr«Iuetn, «r. 
t..lnt;c t.olmiqu .1 for Incorporotim • diwcrot. vortn c.leuUtlon 
in . eutvod W.11 ),t «thod. M»h conrld.rir« th. tru. ..rofoll, 
■ooa further conaideratiom may need to be made. 

^^fect of longit udinal preeaure gradient 

Aa has been aho«n, the wall jet exists in the presence of extremely 
strong positive and negative longitudinal pressure gradients, it 
appeared from the turbulence results, (Section 5.4.3) that the 
initial negative pressure gradient stabilised, or in some cases, 
gradually reduced the measured turbulence levels. As soon as the 
suction peak was passed and the pressure gradient became positive, 
the turbulence levels increased. This may well indicate that little 
vortical energy is dissipated within the flow until after the jet has 

negotiated the suction peak, with the exception of that due to akin 
friction. 


^- 2 - 4.2 Ett«!t of the pre.ence of .hoot u..... 

In .cat fUsht CM It btccuM n.cs.ty to ohok. th. .lot exit to 
obMln .uffioi.nt j.t m... flow. In th... o.... th, J.t 1. und.r- 
•xp«d.d Md CM b.oc«. lodly .up.r.onlc .t th. J.t ..It, Llttl. 1, 
known of th. d.t.ll«i .tructur. of thi. flow, or of coiiipr...lbluty 
.ff.ct. in .xc.pt th.t M •.;pM.lon fM' .pp,„. to .ct. ftob 

th. .lot Up Md th.t .v.ntu.lly, ,t pr...ur. r.tio. .pproMhin, 2.5, 
the jet will detach. 


sine. th. ..loolty dlff.r.nc. .. d.fln«J In s.otlon 6.2.2 r.Min. of 
ord.r, th. .h.ddinp fcqumey .lio r«.ln. of .imilax 
■.smltud., aow.v.t, .Ine. th. .tTMUng iptmS of th. vortic 1. 
9 i».tly lncr....d, th. w.«.l.n,th of th. vottlc b.ooM. Ion,. Ihi. 
«ould t.nd to .nn..t th.t .t .CM point, th. vortic b.c« 
ln.ufflcl.ntly cion to ..oh other to nlntln th. .tuchnot, «d 
hone. th. flow ..pet.., Thl. My .xpUln why th. ohn,. In flow 
•ICUcn occar, to .uddenly ud .1», .i„o. th. lonMi ••v.Unjth. 


rsduoc •ntr«lnm«nt, why th« trantonlo parfomano* dacllnas. 

Many of thaaa pointa wiU ranaln unxaaolvad until furthar axparlnantal 
data baconaa avallabla. 
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7. COWCLUBIONS 

A progrm of raaaarch hat baan oonduotad on tha aarodynamioa of a 
ncnlnal 20% thicknanaiehord ratio oiroulatlon eontroUad aUlptlo 
aarofoU taction aqulppad with loading and trailing adga blowing alota. 

ovarall lifting parfornanca raaulta agraad wall with tha llmitad 
eunount of pravious work availabla. The detailed Invaatigation of tha 
thin, highly curved trailing edge wall jet demanded the davelopmant 
of new teat technlquee and equipment ae well as tha capability to 
obtain repeatable, steady, two-dlmenalonal flow conditions around tha 
model section. The results obtained have led to a new insight into 
the nature of Coanda flow which offers the possibility of a simpler 
theoretical solution to the complex wall jet flow. 

7.1 The Nature of the Coanda Flow 

The results from the detailed trailing edge wall jet Investigation led 
to the proposed discrete vortex flow field, fully discussed in 
Chapter 5. 

Each passing vortex produces a reduced static pressure between itself 
and the adjacent surface, drawing the vortex towards the surface. The 
balance between the ' cent*' if ugal ' force of the stream and the attach- 
ment of the vortices causes the typical wall jet properties of attach- 
ment and growth. The effect of the high frequency vortex stream is 
also to strongly entrain the outer fluid, be it moving or stationary, 
due to the influence of the velocity field of the vortices. 

The assumption that the initial vortex strength and shedding frequency 
(and hence the degree of attachement of the wall jet) are dependant 
vqpon both and the slot geometry has allowed a variety of observed 
effects on overall lifting performance of circulation control aerofoils 
(e.g. slot height, probe interference) to be simply explained. 
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7.2 Tti«orttioal Impltoationi 

Tha propoaad flow flald la bollevad to axplain tha Inability of tha 
axlatlng "tiroa avaragad” thaoratical mathoda (aithar Intagral or 
finita dlffaranea approximations) to satisfactorily pradiet 
circulation control aarofoil or avan highly curvad wall jat flows. 

Tha axisting potantial flow discrata vortax mathod of Smith at al^® 
providas an indication of tha simplicity of raplacing a wall Jat by a 
stream of discrata vortices. Tha present results would tend to 
suggest that a more successful method would include tha effects of 
viscosity, vorticity dissipation and pairing of adjacent vortices at 
some critical point in their development. The potantial for an 
extremely quick calculation scheme still exists however. 

7.3 Conclusions Regarding the Bxperimentil Techniques and 

Instrumentation 

7.3.1 Two-dimensionality of flow 

The attaining of two-dimensional flow, both over the aerofoil and on 
exit from the slot, was of vital importance to the validity of the 
experimental results. The use of some form of secondary blowing was 
shown to be necessary; however experience suggests that the fixed 
nature of the tip Jet sealing strips became a problem when the slot 
height was varied. It was also apparent that at higher blowing rates 
the small tip Jets were insufficient to fully compensate for the end 
plate effects. It is suggested that some form of wall blowing on the 
end plates would be a more universal solution providing sufficient air 
supplies were available. 

7.3.2 Anemometry 

The techniques developed for the hot wire and film prices were shown to 
provide repeatable, accurate results. By using both types of probe, 
.allowance could be made for the frequency response and sise limitations 
of the split film probes. 
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Ih# digital analyaii taohniqiwi ware shown to b« vary aultabla to thla 
typa of data aoqulaltion and axtandad tha applicability of the apllt 
film prebaa by allowing tha flow vaotor to ba at aoma unknown angla 
to tha aanaors, 

Mantlon should also ba mada of tha capaoltanca position control systam 
davalopad during this rssaaroh. it had applications for positioning 
any matal probe in close proxinit/ to any metallic portion of tha 
aerofoil and was unaffected by airflow, its accuracy and repeatability 
ware axcallant and similar systems may ba of significant value in a 
variety of other aerodynamic testing environments . 

7.3.3 Determination of aerofoil effective incidence 

The use of on-line mini-computers and simple approximate theories to 
assess the effective incidence of tha aerofoil proved to be tha 
quickest, most reliable technique developed thus far. The equipment 
also provided instant verification of the validity of the experimental 
pressure distributions and monitored the stability of a given flow 
condition over an extended period of time. 
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8* SUGGESTIONS FOR FURTHER RESEARCH 

Thii axp«rlm«ntal raaaareh hai Indioatad a naad for furthar Invaat- 
Igation of tha vortax straam that has baan ahown to axlst in tha 
trailing adga wall Jat. It would ba datirabla to invaatigate a 
number of diffarant aspacta. 

- tha vortex ahadding frequency. 

- the pairing proceaaea of the vortieaa. 

- the effacta of the geometry of tha slot on the vortex formation, 
the effect of tne severe longitudinal presamre gradients upon 
the vortex stream. 

~ the Influence of free stream and jet flow Mach number. 

Due to the inherently high shedding frequency of the vortices it Is 
suggested that an optical technique, either a stroboscopic achl^eren 
or holography for example, would yield the most uaeful results. 

It would also be of interest to investigate the effects of non- 
circular trailing edges to enable the influence of surface curvature 
on jet attachment to be determined. 

Prom a full scale application viewpoint, further information regarding 
the drag of circulation control aerofoils is required. This is an 

area still grossly deficient in experimental and reliable theoretical 
results . 

Taking into account the suggestions of Chapter 6 regarding the inclusion 
of a discrete vortex method into an aerofoil calculation schema, it 
should be possible to replace the somewhat arbitrary vortex strength 
decay rate from the method of Smith et al^^ by a viscous approximation 
of the dissipation of vortlcity and the pairing of adjacent vortices. 
This should provide an improved prediction method with more relevance 
to the real flow. 
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APPENDIX I i THE SCANIVALVB CONTROL PROGRAM 


AI . 1 Main Program 

The arrays and control paranatera are eatabliahed and a blockage 
correction applied to the input dynamic preaeure. 


1500 

1510 

1000 

C 

1001 

99 


COMMON TCONl/ TC0N2^ DVNP^ SP50<4>> SP<44>, SPAV 
COMMON X<44>> XNEW<44>, SPNEW<44) 

COMMON 2ER01, RLPHfl, CD 
DIMENSION IDftTE<j:> 

WRITE<4, 1500> 

FORMAT <" INPUT RUN NUMBER AND DATE I FORMAT 
READ<4. ieO0>IRUN^ <IDATE<I ), 1=1^ 3;) 

WRITE<4. 1510) 

FORMAT <' INPUT DVNP, TCONl. TC0N2 '> 

READ<4> lOODDVNP, TCONl, TC0N2 
F0RMRT<I4, 312) 


BLOCKAGE CORRECTION 
DVNP=DVNP+1. 0275 
FORMAT <3F 10. 4) 
RUN«FLOAT<IRUN) 
CONTINUE 


' > 


The scanlvalves are checked to their home positions and the program 
awaits an external evnnt {+5 volts d.c. on A/D channel 4) to begin 
sampling. The subroutine STEP, steps the Scanlvalve 1 port; the 
passed argument defines v .ch scanlvalve is to be stepped, -me 
required D/A channel is held at +9.5 volts for 40 ms and then at 
0 volts for 100 ms. The subroutine TDWELL allows a 250 ms delay, 
enabling the pressure at the transducer face to reach a steady state 
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C ZERO D/R CHANNELS 

CALL DAC<O^0. e> 

CALL DAC<l«e, e> 

C HONE SCAN 1 VALVES 

5 Ht1SIQl«ADC<0> 

XF<HMSIQ1. QE. 0. 5) GOTO 10 
CALL STEP<0> 

CALL TDMELL 
GOTO 5 

10 CONTINUE 

C S/Vl HOMED AT FIRST PORT 

6 HMSIG2»ADC<2> 

IF<HHSIQ2. GE. 0. 5> GOTO 11 
CALL STEP<1> 

CALL TDWELL 
GOTO 6 

11 CONTINUE 

C S/V2 HOMED AT FIRST PORT 

C WAIT FOR EXTERNAL EVENT TO TRIGGER SAMPLING 

50 STSIG=RDC<4> 

IF <STSIG. QE. 2. 0> GOTO 12 
GOTO 50 

12 CONTINUE 


Upcm the external event trigger, the program samples the four half 
chord pressures • These were not used in this work. The 44 static 
tappings around the centreline are then san^ilodj allowance is made 
for the fact that tapping^4 was unusable. The subroutine SCANCL 
8anf>les and stores the static pressures! an average of 50 samples of 
each pressure Is used. The Scanlvalves are then checked to have 
returned to the home position. 


oo 
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r 


13 

14 


ee 

ei 

20 

c 


2000 

30 

2001 

31 


li«0^ SCRNNINQ THE FOUR HRLF CHORO FRESSURES 

CflU STEP<e> 

CRU STEP<1> 

CRLL STEP<0> 

CRLL STEP<1> 

CRLL TDWELL 
CRLL SCRN50<1, Il> 

CALL TDWELL • 

CRLL SCRN50<3^I1> 

CRLL STEP<0> 

CRLL STEP<0> 

CRLL TDWELL 
CRLL SCRN50<1, Il> 

CALL STEP<1> 

CALL STEP<1) 

CRLL TDWELL 
CALL SCRN50<3, Il> 

NOW SCAN THE 44 TAPPINGS ON CENTRELINE CHORD 


DO 20 Jsl, 43 

L»L+1 

L±=L/2 

LPsPikl 4 

IF<L2. Ed. L> GOTO 13 

M“0 

N»1 

GOTO 14 

M«1 

N»3 

CRLL TDWELL 
CRLL ZERO<N> 

CALL STEP<M) 

CRLL TDWELL 

IF<L. NE. 4>G0T0 80 

L"L+1 

CRLL STEP<M> 

M“0 

N“1 

GOTO 14 

CONTINUE 

CRLL SCRNCL<N. J> 

CONTINUE 
CALL STEP<M> 

CONTINUE. 

HMSIG2-RDC<2> 

*• *• »> M 

OOTO^jr VALVES OUT OF SVHC',/> 

WR1TE<4> 2001) 

eOHTJHUE HOHEft IN $VNC',/> 
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I 

The Bubroutlna CALCSl lnt«rpolatt« for extra point* between the 
meawured preisurei uelng a Langranglan interpolation technique and 
then integrate! the preeaure eoefflcienta to give and Cq^. 

Subroutine EFFECT allows for eatimation of the effective incidence of 
the aerofoil from the measured pressure distributions emd is shown 
separately in Section AI.2. The main program then outputs the 
results. 


C RUN COMPLETED - PRESSURE COEFFICIENTS NOW EVALUATED 

CALL CALCSl 
CALL EFFECT 

C OUTPUT RESULTS 

WRITE<4. 40O0> 

4000 FORMAT<' DO VOU WANT OUTPUT? i^VESO 
READ*:4/ 4010>IDEC 

4010 F0RMAT<I2> 

WRITE<4* 2002><IDRTE<I>, 1=1. 3), RUN 

2002 NUMBER'/ F5. 1, //> 

2003 F0RMflT<' ' / 5X> ' INPUT PARAMETERS /> 

WRITE<4/ 2004>DVNP/ ALPHA 

2004 F0RMAT<' 5X, ' DYNAMIC PRESSURE ='/ FS. 3/ A 

* 5X, 'MODEL EFFECTIVE INCIDENCE = 'F8. 2/A 

* SX/'C MU TRAILING EDGE •= '/A> 

1F<IDEC. NE. DGOTO 75 
WRITEC4, 3O10> 

3010 F0RMRT<' MODEL STATIC PRESSURES') 

WRITE<4/ 2006)<SP<I)/ I«ri, 43) 

2005 F0RMRT<SF9. 3) 

WRITE <4/ 3050) 

3050 F0RMAT<' INTERPOLATED PRESSURES') 

WRITEC4, 3060><SPNEW<I)/ I-l, 43) 

3060 F0RMAT<6F9. 3) 

75 CONTINUE 

WRITE<4/ 2010)SPAV 

2010 F0RMAT<' '/ 5X> ' %»*<•♦<«♦••*• LIFT COEFFICIENT 

* 5X/ 'CL - '/ F10. 4 / /'> 

MRITE<4, 2O20)CD 

2020 FORMAT <6X/ '*#»♦•>♦<•»« DRAG COEFFICIENT *•••*•*'•;'/ A 5X/ 'CD«'/ F10. 4/ /> 

GOTO 99 

STOP 

END 
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AI . 2 Sub rout! n« IP FBCT 

Thii .ubroutin* initially caloulatai tha laadin? adga lift ooaffioiant 
(aa far aa x/c ■ 0.5) from tha axparimantal raaulta. 


C 


1Q0 


150 


210 

220 


200 


250 


1500 


SUBROUTINE EFFECT 

TO FIND EFFECTIVE INCIDENCE OF MODEL 
COMMON TCONl* TC0N2/ DVNP> SP50<4>, SP<44>, SPAV 
COMMON X<44># XNEW<44), SPNEW<44) 

COMMON ZEROl, RLPHfl> CD 
DIMENSION THETft<44).CP<44> 

SUM-0. 0 


DO 100 1-33/ 43 

KFflCT2-<XNEN<I)“X<I)> 

SUM-SUM+XFRCT2*<SPNEW< I )+SP< I > 0 

CONTINUE 

CONTINUE 

DO 150 I-l, 11 

XFnCT«<XNEW<I>“X<I )) 

SUM-SUM+XFfiCTH.<SPNEWCI>+SP<I>)^2 0 

CONTINUE 

DO 200 1-33/43 

Il-I+l 


IF<I. EQ. 43>G0T0 210 

GOTO 220 

Il-l 


CONTINUE 

XFfiCTl-CX< I1>-XNEM< I ) > 

SUM«SUM+XFRCTl»t*<SP< I1)+SPNEW< I > >/2 0 

CONTINUE 

DO 250 1-1/10 

Il-l+l 

XFflCTl«CX<Il>-XNEM<I>> 

SUM-SUM+XFflCTl»H<SP< ID+SPNEWC I ) >/2 0 
CONTINUE 

SPRVl-SUMXll. 712 

CL—SPRVl 

CLl-CL 

WRITE<4/ 1500>CL 

FORMRT<' CRLCULRTED LEADING EDGE CL -%F6 3> 


Tha axparimantal data ia than plottad on a C»T diaplay at tha tunnal 
conaola. 
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CONTINUE 
CALL CLRPLT 

CRLL SCALE<e, I, 24* -2. 0O> 

CALL PLINE<e^e#24.e> 

CALL PLJNE<e.i.e«-2> 

PLOT EXPERIMENTAL DATA 
DO 900 Z-23< 42 

CALL PLINE<X< I >« SP< 1 >/ XWEW< I >, SPNEW< I > > 

CRLL PLINE<XNEW<I>,SPNEW<I),X<I+1>,SP<1+1>> 

CONTINUE 

DO 950 I»l. 22 

CALL PLINE<X<I>, SP< I XNEWd >> SPNEW<I)> 

CALL PLINE<XNEW<I), SPNEU<I>, X<I+1), SP<1+1>> 
CONTINUE 


The operator inputs an estimated effective incidence and the program 
iterates around the overall theoretical lift coefficient (using the 
equations given in Section 3.1.1) until the theoretical and 
experimental leading edge lift coefficients (0 < x/c < 0.5) are 
tolerably close . 

360~ WRITE<4, 10O0) 

1000 FORMAT<' INPUT ALPHA EST" > 

RERD<4> 2000>ALPHA 
2000 F0RMAT<F6. 3> 

300 CONTINUE 

DO 310 1=1,43 

THETA<I >=ACOS< <2, 425>-i > 

IF<I. LT. 22>G0T0 310 X. ? 

I”ETA<I)=2. H.3. 14159~THETA<I> 

310 CONTINUE 

RLPHA=ALPHH/57. 296 

311 CONTINUE 

DO 320 1=1, 43 
DUMMV«CL/’<2. m3. 14159>t<l, 2) 

C AMMA=AS I N < DUMMV ) -ALPHA 

JSf "?■ jt'f > >*siN(nLPHR+enMMft > ) 

e*THETfla» 

320 CONTINUE 

TSUM*0. 0 
DO 600 1-33,42 
XFRCT=X<I+1>-X<1 > 

TSUM-TSUM-fXFRCTf*<CP<I>+CPCl+4 1 wo a 
600 CONTINUE I+l) )/2. a 

DO 650 1=1, 1- 

XFACT»X<l-fl)“ ,;<i > 

TCL— TSUM/11. 712 

goto 311 
610 CONTINUE 


non 
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Th. th,or6tlo«l pr...ur. di.tribution 1. th.n •up€rimpoi«d upon th. 
exp»rlro«ntal results and the operator evaluates the natch visually. 
A satisfactory incidence is transferred to the main program upon 
conclusion. 


PLOT THEORV 


3S0 


346 

iieo 

1260 


DO 340 1=1, 43 
CONTINUE 
C1“CP<I>- 05 
C2=CP<I>+. 05 
C3=X<I>+, 25 

C4=X<n-25 

CALL PLlNE<C3,CP<l>,C4,CPa>) 
CALL PLINE<X<1>, Cl, X<I>, C2) 
CONTINUE 
WRITE<4, 1100> 

FORMAT <' 0. K. =1" ) 

FEAD<4, 120O>1DEC 
F0RMAT<I2> 

IF<IDEC. NE, 1>60T0 312 
ALPHR-^ALPHA^S?. 296 

return 

END 
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^PENDIX II I HOT WIRE AWBMOMBTRY 

P*»^^vatlon of a Hot Wlr » Ambient Temperature Correction 


The tenqperature correction given in reference 44, equation 72 is, 

+ B(PU)"(1 - ^ At ) 

where 


a m 



overheat ratio 


At » t - t 
n m me 



and the subscripts are 


w ■ wire 

m ■ ambient fluid 
me - ambient fluid at calibration 
then, u can be given by 


a 


<’’w - ’’m) 


(T 


T ) 
me 



continued 
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AZZ.2 Derivation of the Equations to Oiva U and V from an X-array 
Hot Wire Probe 


Aeauming the coordinate system and conventions shown in figure 49, 
and that represents the effective cooling velocity on a wire, 
them 


« (Usin0 - Vcos0)^ + 


%long " VsinO)^ 


assuming 


where 


r .^ -2 ^22 

%ff * %long 

a ” wire direction sensitivity coefficient 

• 0.2 


Then, 

2 


- Vcos0)^ + w‘ 

i.e 0 ' 

- 90" 

2 

%ff • 

+ w‘ 

1 , i.e. 

0 ■ 0" 

**eff " 

{“V) - 


2„2 


2.2 


continued 
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Cro»i Wlr> A nalyxi Coaputar Program 

Dub to thB conplBxlty of tho oquatloni dsrlvBd In tho prtvlouB 
BBotion, thB antlyBlB of tha dual aanaor hot wira signala waa 
parformad in two parta. Program XRUN aanplad tha aignala, convartad 
tham to floating point nudbara and than wrote tha raaulta, including 
run paramatara, to file, on a hard di.c A aacond program, XSUMS, 
waa then uaad at a latar time, to recall the data and parfonn tha 
analyaia and time averaging. 

Ml. 3,1 Program XRUN 

initially the program defines the various file names to hold the run 

data and samples, and gathers the required mean and r.m.a. voltages 
and flow temperatures. 


continued over 
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or poo* ou*UTf 


f' 

Irf 


1000 


1010 


c 

c 

c: 

10 


THISJ ZQ A DUAL OHANNEL QAMP'LINfJ 
EXTERNAL XRUNHB 

COMMON /BLKl/:m,Jr<SlUO> rIDUFX(Sj,;‘) 

COMMON /BLKlVIDEFrlCMF »IBUr4.<512) 

COMMON /DLK3/ICIIAN1 r ICI IAN2 
COMMON /BLK4/El<2!5<i>»E2(2nA) 

COMMON /CAIN/IRUNr II. »>flCLAArX2CLAArXlDTAA»X'»riTAA 

«REm'; ‘ "^2 - "2 - BIOMaJ . TC^:.2 , 

REAL1C4 HTl (2) »HT2<2) »PRnGNMf2> 

REAL NlrN2 

DATA HTl/<f*RDKlXlAr<f»R DAT/ 

DATA HT2/4RDKlX2ArAf? DAT/ 

DATA PROONM/ARDKOXRUr ARN SAU/ 

CALL RD lA IN ( 1 1- LAO r ,T RUN r 60 ) 

OAMPLINC PROGRAM' »*0) 

XUrlrtNl^IGb. I L ( ) 

IFdOHANl ,LT,0)0TaP 

II" ( II- T: f CH ( 1 1 ri ) » NE ♦ 0 ) GTOI'- 

irn:ENTERacilANl»HTl»20)*LT*0)ST0P 

ic:han2:=igetc<) 

Ir(ic:(-IAN2.LT.0>E;T0P 

ir ( IFETCH ( I IT2 ) . NE , 0 ) STOP 

IP ( lENTER ( IP IAN2 r UT2 » 20 ) . LT * 0 ) STOP 

IRUN^O 

^^I■^ ''■NAI-VSIS PROORAM',-0) 

CALL PRINT ( INPUT TOTAL NUMDIIR OF RUNS'p'O) 

READ(7f 1000)NT0TAL 
FORMAT (I 3 > 

r*?An / EOl r A1 r N1 » OIGMAl r TCALX ' p ’0 ) 

RLAIK /r lOlOEOl »A1 »N1 » GIOMAl r TCALl 
format (5F0. 3) 

CALI PRINK' INPUT E02» A2rN2rCIGMA2rTCAL2' » *0) 

READ ( / » 1 01 0 ) E02 r A2 » N2 » 0 1 0MA2 r TCAL2 

CHAIN LOOP COMMENCE*:; 

LUN1»14 

__LUN2«15 

“Ti»iiTr — ‘ . ■■ ■■ - 

IRUN==«IRUNM 

CALL PRINT<' INPUT RUN NUMBER' p»0) 

READ<7rl000)TRUNl 


f 
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OP POOR QUALITY 


CrtLL PRINK' INPUT rRMJSt » T1 ' » '0> 

READ < 7 » 1 030 > EMt^ANl i FitNK 1 » T 1 
J030 |•ORMAT(:'^^H*^) 

CALL PRINT(' INPUT IL'ME:AN2rRNC2» T2' » *0) 

READ < 7 » 1 030 ) EHEANl’ r RMC2 r T2 

LUN3'«il2 

LUrM'«13 

CALL ACS I GN ( LUN J, r X 1 CLA A > 

CALL ASSIGN <LUN2» XI DTAA> 

CALL ASSIGN ( LUN3 » X2CLAA > 

CALL A«SJ:GN<LUN^rX:>DTAA) 

IREC1,”J. 

IREC3-1 

DEFINE FILE LUNl (3»:>5A» 'J» IRECl ) 

DEFINE FILE LUN2<20»Sl2i.Uy IREC2) 

DEF I NE F I LE LUN3 < 3 » 2lj A » Li » I REC3 ) 

DEFINE FILE LUN4 ( 20 » 1512 » U» IREC4 ) 

URI TE ( 7 » AOOO ) XICLAA r X2CLAA 
AOOO FORMAT (' ' » 1L5A1 * 5X » 15A1 > 

WR I TE < 7 » A 1 90 ) EO 1 ^ A 1 1 N 1. y S I GM A Ir TC AL t » I RUN 1 1 EHEANl ► 
♦RMSlrTlrNTOTAL 

A 1 90 FORMAT < 5FL5 ♦ 2 1 1 3 y 3FLi * 2 y 1 4 > 

C 

t: WRITE OALIDRATION CONSTANTS TO FILES XICL** » X2CLJ»:<c 

C 

WRITE CLDNl ' 1 >E01 y Al y N1 y SIGMAl t TOALl » IRUNl rEMEANl » 
J|!RMS,l rTl rNTOTAL 

WR I T £ ( LLIN3 ' 1 ) E02 » A2 y N2 y S I GMA2 y TO AL2 r I RUN 1 » EME AN2 » 
)KRMG2rT 2 rNTOTAL. 


The prograun aeunples the input signals upon a carriage return by the 
operator. The completion routine XRUNSB writes the samples to the 
previously defined data files. 

i 

C 
C 

c 
c 


PAUSE nrPE <CR> TO SAMPLE DATA' 

SAMPLE DATA 
ICMF'^0 

I CMFl ^ ^*UF r 51 20 r 1 0 r 25A0 r 0 r 2 » 1 » 2 » I CMF » 1 DEF r 25A t XRUNSDj, 

CALL SEfR<3rGHr5.0rICMFl) | 

CALL PRINK' WAITING RTS'r'O) 

CALL LWAIKICMFrO) 

CALL SETRL'lryr) 

CALL PRINK' FINISHED SAMPLING' » '0> 
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*I1ia latnplta ar« rscallad from th« fllaa and convartad to floating 
point niunbara. 


DO 100 >lrlO 
K«'J *:>-1 

inooeii;;'«:i:reaeiw ( ai 2 r .t bui- j » k » iohani ) 
I^K;orM^sIRr:ADW<K^2^ iBur2»Kr i:rHAN2> 
rci 11.0 

El < M > ::■ < rtn A < iBuri < m ) ) 2 . > *2 . s/ai 2 ► 

E 2 < M > < I •L.T 1 A < I BUI"2 < f1 ) > 5 1 2 . > <«2 . 5/a 1 2 . 
110 CONTINUE 

WRITE <I..UN2'J) El 
Wr<ITEa..lJN^,'J)E 2 
100 CONTINUE 

CALL, ^■•|^•INT(' BATA WR,TTTEN TO FILES' f* 0 ) 


Upon, completion, the files are closed. The names of the data files 
are updated (this provides each run with a dedicated series of file 
names) and the program chains back to itself to begin a new satiple 
point . 


ORIGINAL PAGE 18 
OP POOR QUALITY 


C(UI. (LUND 

CALL CLC)i;i:(LUN2> 

CALL CLCiLi:; (U)n;^> 

CALL C:L0CI.:(I,.UN4) 
IF(n*na«2A>(30T(3 2Q 
GOTO so 

no i;>o 

X 1. D T A A < 1 0 ) :l nr AA (1 0 > •• • 1 
XI CL A A (10) -=x 1 C:L AA (10) ' 1 
X2nr AA (10) »x2nr AA ( i o > * i 
X2CLAA ( 10 > '«X2C1..AA ( 1 0 ) ■■ • 1 
CONTINUE 

X 1 nr A A ( 7 > -x 1 OTA A ( 7 ) I " 1 
X1CLAA(7)~X1CLAA(7) » ‘1 
X2tiTA A ( 7 ) «X2nTAA ( 7 ) J * 1 
X2CLAA ( 7 ) - X2CLAA ( 7 ) i " 1 
I1»0 
GOTO ^0 

XinTAA(10>'=XinTAA(10)i':l, 

XI CLAA ( 1 0 ) ^^X 1 CL AA (10) !• • 1 
XanTAA ( 10 ) >==X2DTA A < 1 0 ) !M 
X2CL AA (10) ‘^X2CLAA ( 1 0 ) i " 1 
CALL IDELE r (ICHANl »I-iri ) 
CALL HiELET ( ICHAN:M-IT2 > 
CALL CLOSECaCHANl) 

CALL CLOSECdCHANS) 

CALL IFREECaCIIANl ) 

CALL IFREITCdCHAN:?) 

IRt:Cl*d 

IREC2‘=LI 

IREC:5"1 

IREfd«l 

IF ( IRON . I.‘.0 , N TOTAL ) STOP ' TH 
CALL C:HA IN < PROGNM , IRON > > 

STOP 
ENn 


ri-IATS ALL FOLKS! I 


The completion routine XRUNSB Is called every time 1 buffs, of data 
(256 samples per wire) Is collected by the CALL RTS Instruction. 


POOR 


XRUN«EUF'(3R WITH XRUN.GrtV 
GUDROUnNI" XRUNGI* 

COMMON /IJLK 1 / 1 BUr ( '.'J V20) a EKJI ■ 1< IB 1 2 > » I HUr 2 < tB 1 2 ) 
COMMON /BLIC’/IBIl'l'' » ECMF 

COMMON /bi,.k3/i;c:nan:i » eciian;* 

DATA ICALL/0/ 

I'CALL'-ICALLll 
DO 100 

M-2KA!|!( lOALL - 1 >i( ,t#2 - 1 ) 

IDUri<I)»IDUr <M) 

iBLinia > -i'DLir(Mii > 

100 CONTINUE 

IDLIM=aCALL>|:2-l 

i: NCODi;*" I un i ru < ii 12 fX duf 1. » 1 nuM •• i chan 1 ) 

INCODE-a: UR I TU ( :> 1 2 » I DI.IF2 P .T: III )M * :c Cl IAN2 > 

IDEF>=:IBIII"1 1 

RiniiiHi 

END 


All. 3.2 Program XSUMS 

This program performs the analysis of the results stored by program 
XRUN. 

Initially the data files created by XRUN are redefined and made 
accessible. The run constants are read in and the various required 
parameters calculated. The accumulators are all set to zero. 



IS4 


j 


omaiNM. MQi II 

or POOR QUALITY 


XflUMn.F'Cll^ NJW '7V 
CROSS WIRC rtNALY«n{ 

OOHMON/riLKl/NTirr AI,. » .U » i:^>^Ui(«iAATX10ArtA^X2DAAA»X2CAAA 
CCJHMnN/rH,.K2/l::i (2'JiA) riraC’l'iA) .ll:l (MA) »U;i(2IS6> 

cc)hmon/plt/uimr:an » u:;Mr::AN » m„AO r i< 

RC:AL>K4 NiRUN»N2RI,)M 
REAUl;ri E0MUSf!:SMI.,20 
REAL NlrM:,( 

REAL.t>l PR£)0NM<2> 

DATA FT<tlC)NM/ARriKOXBl)»ARM!:; fJAVX 
CALL RCIIAl'N(iFLAO»NTOTAL^AO) 

IF(irLAG,LT,OX3(JTQ 10 

LOGICAL*! XlDAAAaS) »Xl(;:AAA<ir.5) »X2DAAAa!5> »X2CAAA<13) 

DATA XIDAAA/'D' » 'K' r ' 1 ' t ' ! ' t 'X' i ' 1 S '11' » 'AS *B' » 'A' i 
*' . ' » 'D' > 'A' p 'T' T '0/ 

DATA XtCAAA/'D' p 'K' p '1 ' p ' t ' p 'X' p '1 ' i 'C' » 'A' » 'DS 'AS 
*' ♦ ' p 'D' p 'A' p 'T' p *0/ 

DATA X:iDAAA/'D' p 'K' p ' 1 ' p ' { ' p ' X ' p '2 ' p 'D ' » 'A'p 'B'p 'A'p 
*'.'p 'D'p'A'p'T'p'O/ 

DATA X2CAAA/'D' p 'K ' p ' 1 ' p ' { ' p ' X ' p ' 2' p 'C' p 'A ' p ' B' p 'A' p 
*','p'D'p'A'p'T'p’0/ 

CALL PRINK' DATA RECALL PROGRAM ' p *0) 

I2"=0 

n'=i!+i 

I2^I2H 

LUNl-14 

LUN2’=1S 

LUM3--1.2 

LUN^<^!3 

CALL ASOION < LlJNl p X.1.CAAA ) 

CAL L ASS I ON ( LUN2 p X 1. DAA A ) 

CALL ASSIGN < LUN3 p X2CAAA ) 

CALL ASSIGNL LUN4 p X2DAAA ) 

DEFINE FILE LHN.1. ( J, p 2K<f> p II p IRES 1 ) 

DEF I NE F I LE LLIN2 ( 20 p 5 1 2 p LJ , I REC2 ) 

DEE I NE F I I E LUNl^ ( J. p 25 A p U p I RLTIS ) 

DEFINE FILE LUN4(20pi512pl,lp IRIIC4) 

READ (LUN! ' 1 >E01, p A1 p M! p SIGMA 1 p TCALl p IRUNl pEMEAN.1 pRMSI pTI pNTOTAL 
READ ( UJNU ' 1 > E02 p A2 p N2 p SKIMA2 p TCAL2 p IRUN J p EMEAN2 p RMS2 p T2 p NTOTAL 

SET UP CONSTANTS * - 

TCALl«TCAI..!i27;5, 

TCAL2‘-TCAL2+27;?>. 

N1.RUN'«1,/NJ. 

N2RUN==:1,/N2 

E0JRUN:E01**2 

E02RUNi:.E02**2 

Tl--T,ti27;,;. 

T2'-T2i27;T, 

EPSLN ,1. s: ( r 1 • T CAI.. 1. ) /TCAL. ! 

EPSLN2 * < T2-TCAL2 > / TCAL2 
CORR.t«l ./< 1 .iEPSLNt/SIGMAi, ) 

C0RR2«1 ./( J ,iEI*SLN2/SlGMA2> 

A1RUN-A1*CC)RR1 

A2RUN:^:A2*CnRR2 
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ORIGINAL page 18 
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n:i>AVi-o . 

E.'lF'RS=^0, 

E2PRS-'^0 . 

E12PR=«0. 

U1MEAN«'0, 

LI2MILAN-S-0. 

URM£!1«0. 

URMf{2-0 » 

U,12:=0. 


The mean and r.m.s. values of the input signals are calculated and 
then used to obtain the attenuation due to the unknown input 
amplifier gain. 


C 


no 

100 


fn ART CALCS 
DC) 100 J-1,10 

READ ( LUN2 ' J > ( EK J) , 1 :=t 1 , 256 ) 
FJEAD ( UJN4 MX E2 ( I ) , I y 2V)4 > 
DO 110 M=--=l»25/> 

ElA'v» -ElAVH;;i (M) 
E2AV-C2AViE2(M> 

E 1 PRS^^Iil PRSiE KM) XcE KM) 

E2PRS:^ E2PRST)i:2 < M ) >f:E2 ( M ) 
COMTINUE 

UR I TE < LUM2 M X E K I ) » I 1 ► 25 A ) 
WR IT E ( U.JM 1 MX E2 Q ) y K: K 25 A > 
C0NTI)R)E 

Ell>RR:.'^Eir-Rr>/25A0» 
E2PRS^=E2PRS/25A0 . 
E1AV=^E1AV/25A0, 
E2Ay«=E2AV/25A0. 

E IPRSr^^SQRT < E 1 PRB-- ( E 1 A 1 A^ ) ) 
E2PRS-S0PT ( E2PRS - ( E2AV^?E2Ag ) ) 
ATTI^Nl-RMiil/ElPRS 
ATTEN2-RMS2/E2PR3 


The attenuations are used to scale the samples and the 
derived in All. 2 can be used to obtain U, V. 


equations 
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DO ;120 

REAti<LUN2'J)<Eia)».T-«1 ri>56) 

READ<LUN4'J) (Eaa) 

DC) 130 M*«lf2SA 
El(M)!«(El<N)-r;:iAV>>KATTI-Nl 
E2 ( M > = ( E2 ( M ) ~E2A V ) )):ATTEN2 

E£5ML1S==EXP<N1RUN:KAL0C)( < CEMEANH El CM) >*.t2- EOIRUNJ/AIRIIN) ) 
ESML2£5“EXP(N2RUN){{ALOR( ( (EMEAN2+E2CM) >*X!;»-E02RIJN)/A2RUN) > 

E 1 2PR>=E 1 2PR < E 1 ( M ) :f;i;:;2 < M > > 

ESML 1 S~ESHL 1 S^ESML 1 £) 

ESML2£>«ESML2S*E0ML2f; 

EP-^ C ESML 1 £)+ESJM(..2£J ) / 1 , 04 
EMi= ( ESML 1 K-E0ML2S ) /O * 9*4 

U1 (M ):= (SORT (EPiSQRT ( ADO (EP*)K2-EM!|cEM) > ) >/l .4142 
U2CM>-“-EM/(2.ScUl<M)) 

I UlMEAN:=mMEANHJl(M) 

j U2hIIAM ■i=LI2MEAN f U2 ( M ) 

130 CONTINUE 

WRITE < LUN2 ' J ) < Ul (I ) , I;=il , 256 ) 

WRTTECLUrM' J> (U2 ( I ) » I-l » 25/> ) 

120 CONTINUE 

U 1 M E A N - U 1 ) i E A N / 2 5<fi 0 . 

U2MEAN -U2MEAN/2H<r.O . 

E12PR^^E12PR/25A0. 


Prom the data, now stored as instantaneous U emd V components, the 
required turbulence parameters u'^, v*2, u'v' can be calculated. 


DC) 140 JsltrlO 

READ ( LUN2 ' J) < U1 ( I ) » :«^ » 25o > 
READ ( LUN4 MX U2 <1 ) » I - U 2 . 56 ) 
DO ISO M'=1»2S6 
UIPR -U1(M)“U1HEAN 
U2PR‘--'U2<M) -t)2MEAN 
URMS 1 -URMC 1 iU 1 PR»U 1 PR 
URMr;2 HIRMC2 I U2PR:|tU2PR 
U12«-U12iUlPR)|cU2PR 
ISO CONTINUE 

1^0 CONTINUE 

URMC :l, -‘'SORT ( URMSl /25A0 . ) 
URM52«S0RT a)RM«2/2S<f)0. ) 
U12:=U12/2BA0, 
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The results are output; the data files closed and the file names are 
updated ready to find the next set of stored data. 


IFarLAO.LT.OGOTCl 20 


WR .T TI-: < 7 y 1 000 U-0 1 * A 1 » N 1 ; £5 1 OM A 1 » TC AL l»i:02fA2yN2,S 

URrrE(7yll00) 




c;c)nti:nliii 


UR rriz^ 7 y 1 20 r t.;i. MIIAN y i;2M£ AN y URMSl y URMS2 » U 1 2 y El 2PR 

CALL. CLOSE (LLJN2) 

CALL CLOSE (LUN;3) 

CALL CLOSE (LUN4) 

IF(12J=:0.2<f))G0T0 30 
GOTO 40 
no SO y25 


X 1 DA A A ( :l 0 ) ."X 1 Li AAA (10) • • 1 
X 1 C A A A (1 0 ) '- X 1 C A AA (10) • • 1 
X2DAAA ( 10)=-X2riAAA ( 10) •• * 1 
X2CAAA (10) -X2CAAA ( 1 0 ) " J 
CONTTNUE 

X 1 DA A A ( y ) X 1 DAAA ( 9 ) f * 1 
X ICAAA ( 9 ) =-X 1 CA A A (9)1*1 
X2riAAA ( 9 ) "iX2l'iAAA. ( 9 ) 1- * 1 
X2CAAA ( 9 ) ••=X2CAAA ( 9 ) 1- * 1 
r 

GOTO 60 

X 1 DAAA ( 1 0 ) -0< 1 n AAA (10)1* 1 
XICAAA ( 10 ) -=X1CAAA ( 10 ) f- * 1 
X2DAAA ( 10 ) :^=X2riAAA (10)1*1 
X2CAAA ( 10 ) :=X2CAAA (10)1*1 
ir (Il.EQ.NTOTAL)STOP 
CALL CHA I N ( rR0(3NM y NTOT AL » 6 0 > 


> F0RMAT(/»10Xy »GXr 'DIGITAL HOT WIRE ANALYSIS ' yOXy 
>!c'*mt'y/y5Xy'INPUT CALIDRATION CONDITIONS TOR WXRF t'?/' 

*5X » ' ZERO VOLTS ' r 1 5X * FS , 3 y / y 
#SX y ' INTERCEPT ' y 1 AX y FS , T5 y / y 
#SX y ' SLOPE ' y 20X y FS , 3 y / # 

' OVERHEAT RATIO ' y 11 X y FS . 3 y / » 

' CALI DR AT I ON TEMPERATURE 'y2XyF6,ly/, 

'INPUT CALIDRATION CONDITIONS FOR WIRF 2'y/y 
'ZERO VOLTS' ylSXyFS. 3 y/y 
' INTERCEPT' y 1 AX yFS. 3 y/y 
'SLOPE' y20XyFSy3r/y 
' OVERHEAT RAT 1 0 ' y 1 1 X y FS » 3 1 / y 
*SXy 'CALIDRATION TEMPERATURE' yZXy FA. 1 y/y > 

FORMAT ( //yl AX y ' y OX y 'CROSS WIRE ANALYSIS' 

♦ y 1 3X y ' JKtJU*# ' y /5X y ' RUN ' y 5X y ' 111 MEAN ' y 5X f ' 112111:; AN ' 

*y3Xy 'URMSl ' ySXy 'URMS2' y4Xy 'U12' ySXy 'E12PR' ./) 

FORMA r ( 3X y 1 3 y 3X y FO » 2 y 3X » FO . 2 y 5X y FO ♦ 5 y 3X y FO . S y 2X y FO ♦ 5 y 2X y Ffl . S ) 
B ro!** 


#RXy 
ilfSX y 
tSXy 
*SXy 
JitSXy 
^KSXy 
*SXy 
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All. 4 The Split Film Analyslw Program 

The split film analysis program (TSIRUN) Is essentially very similar 
to XRUN, XSUMS except that the two functions are combined. The 
program does not create stored data files l.ut processes each batch as 
[ It Is taken. 

The program Initialises the spilt film constants and defines a 
temporary data storage file. The run constants are Input (E. , E-, 

/e j/e^'^ ) and the program then waits to start sampling. 


COMHON /ELK 1 / 1 HUF ( '.fi 120) 

n I ML’MEt ION i: III 0 ( 5 ) » ira g i ( » n:.jm ( 5 > 

DIMENSION [i: 1 < 2S AO ) y E2 ( 256 0 ) 

REAL. N1 » N2 . Nl. RUN » N2RUN » N » NUM 
WR;[TE(7rlOOO> 

READ (7» 2000 ) i;;o:i »E02 
A 1-41 ,77 
N1-0,3E)A 
SIGMA1==,5 
A2==9,U5/> 

N2:^.05A 

SIGMA2i--',5 

J:^l 

NTOTAL-^l 

LUN==3 

CALL AG:.;iriM(LUN^ 'DKl JGPLITl .DAT' ) 
DEFINE FILE LUN<H00v20rli* IREC) 

10 CONTINUE 

WRITEC7»1100) 

READ ( 7 » 2100 >EMEAN1 »RMG1 » EMEAN2 rRMG2 
WRJ.TE(7»3000) 

3000 FORMAK' IB INPUT O.K. 1™YES'> 
READ<7»3100)IBWCH 
3100 FORMAT (12) 

IFdGWCn.NE.DGOTO 10 
PAUBi:; 'TYPE <CR> TO SAMPLE' 
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The signals art sampled and converted to floating point numbers. 


ICMF»»0 

CALL K-TEJ ( 1 BUF p S 1 20 r t » 2S<f»0 0 » 2 » 1 r 2 r I CMF p I BEF p 2SA0 p ) 
ICMF1<=0 

CALL St;TR<;3pBllr5»Opi:CMF.1.) 

CALL LWAITaCMFpO) 

CALL S£TR<-1ppp) 

DO 100 I;=1»25A0 
M»l!K2-l 

E 1 a ) == < F L.T 1 A ( 1 DUr ( M > ) • 5 1 2 » > >|c2 , lVL=i 1 2 . 

E2 <1 ) « ( FLT 1 6 ( IBUF < M f- 1 ) ) -S 1 2 . ) 5JC? . 5/5 J. 2 * 

0 CONTTNUII 


The attenuation of the input signals is calculated. 


H1RIIM = 1 ,/NJ. 

EO J, 2 -' ( FOl ^ , 92 ) { ( EC2;1; . 92 > 

ElA0:=:0. 

E2Ak>:=0, 

E1FR0:::0. 

E2PR0^-0. 

DO 200 M=l»25('.0 

ElAV='.^ElAOfEl(M) 

E2AV~E2A0Hr2(M) 

E1PRS=E1PRS+E1(M)>KE1(M) 

E2PRJL«E2PR0+E2 ( M ) ^!cE2 < M ) 

CONTINUE 

P*25<40» 

ElPROirElpRG/S 

E2PR'>™E2PR0/S 

ElAVi«ElAV/£; 

E2AUwE2AV/0 

E1PR0-£;0RT<E1PR£{ - (E1AV>|!E1 AV) ) 
E2PR£;«sqrt < E2PRS E2AV*E2A V ) 
ATTENl-RMta/ElPRS 
ATTEN2'=RMG2/E2PRS 
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The input data la scaled by the attenuation and the values for U, V 
are calculated. 


DO 300 M‘=J.r2SA0 
El(M>»(P:nM>-El,AV)!HATTENl 
E2 ( M > !* ( E2 < M > • E2A0 ) * ATTEN2 

300 CONTINUE 

DO -100 M«1»25A0 

ElSa='(i:MEANliE;l<M) )5f!*2 

E200- ( EMI-:AN2 » E2 ( M ) ) >K*2 

IIN'“< (E1£!CM E2S0 - E012)/A.t )«!(:N1RUN 

E2<M)==(E10Q -E2S0)/<A2H<(UN5!:5>:CN2-1. ) > ) 

El < H > -£5QRT ( UN*5!<2-E2 ( M > >}!5li2 ) 

400 CONTINUE 

UAU^’^O , 0 
VAV’--0 ♦ 0 

DO 500 M“1»25A0 

UAV-UAOiEKM) 

uau-vaome;2<m> 

500 CONTINUE 

UAU===UAV/5 
VAV'^UAV/S 


The velocity fluctuations are determined and u'^, v'^ and uW simply 
obtained. 


UP-^0 . 0 
0 . 0 
UV-O ♦ 0 

DO AOO ,i=lr 25^*0 
Ul-EKM) UAV 
V1-^E2<M)-VAU 
UP«iUPiUU:!<:2 
UPi^UPi0l!|c>|f2 
UVHJUmU:Vl 
600 CONTINUE 

UPRr.j;;uRr(UP/‘;> 
vpr<--^aar<T(VP/o) 
UVPR = U0/!!5 
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Th« value of the traverelng gear counter ii obtained to give the 
radial location of the neaaurement point. 


IDIG<1)-»20 

IDIG<2>«M0 

IDIG<3>«"dO 

IP1'CI<4)»*100 

I <=1 


I N = I n I R ( 1 1 . • 1 ,i 00 1 4 ► • 1 77 » 1 » 1 VAR ) 

1 70 I nOR a r " 70 » • J, 77777 > r N y ;i VAR ) 

I Cl 1 0 1 ( 1 ) := I D T I U i , • 70 y ■ i A 0 » 1 y I V AIV > 
IF ( iniGi a ) . NE . iriiG a ) ) gotid i 
IN-ICilR< J. t »lA00J4y ■177y 1 » IVAR> 

1 72:.= I DOR ( J. y • 72 , • 1 77777 , 1 N , T VAR ) 
IDKM ( I ) == r DIR a y • 72 y - 1 AO y 1 y I VAR ) 
IF ( IDIGI ( I ) . NE» IDIG ( I ) ) GOTO I 
IN.‘=IDIR<1 y "lADOlAy M.77yl y TVAR) 
I74==ID0Ray *74y •i77777yIHyIVAR) 
IDIGI n ):=IDIR( 3. y “74y ■ .lAOy 1 y TVAR) 
IF < IDIGI (I ) . NF . IDIG ( I > )GnT0 X 
INUM:.-: I DIR (1 y • 72 y • 3 7 y :l y I VAR ) 
NUMCD^FLOATaNliJI) 
NUMa)«i5..-NUMn:) 

I~IH 


IF(I J:lo.A)GOTO 2 
GOTO 1 
CGNTINUF 


N'-NUM ( 3 ) M 0 ♦ -tNUM <2)1100. 
N~NH0000*WUM(5) 


:f:NUM(;5)H000»*NUM<4) 


The results are output and written to the temporary storage file and 
the program loops to enable the next probe position to be sampled. 
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WRITE < 7 r nOO > UAU » UPR » VA V » VPR » tIVpR , NTOTAL » N 

NTC)TAL«NTOTALil 

WRITE(7»1200) 

READ<7f2200)ItiEC: 
iFaDEc,Ea,2>Gtrrc) 30 
iF(iDEC,Ea*o)otrro 20 

WRI T E ( LUN M ) IJA V » UPR r VAU f VPR t UUPR f J » N 
J«JH 

CONTINUE ; 

GOTO 10 

*»^™F^rE0lfE02^ sampling AND PROCESSING PROGRAM',/ 

FORMAT (2Fci» -4)*” 

FORMAT ( ' INPUT EMEANl p RMS J, p EMEANl p RMS2 ' ) 

FORMAT ( 4FB « A ) 

»='^ES,0-NC, 2=EXIT and URITE<) 
r < tiX p 'UDAR’“ ' p I" 8 . 2 p 2X p ' UPRIME« ' p F8*3 p2X, ' VE<AR= * > 

*J;B;3;2X,'Vf.RI«E.'H-fU3./,SX,-U. = -,W.3,2 -“S, si. • . 

WRI TE ( LUN ' J ) UAU p UPR , VA V p Ul >R y UVPR , J y N 
U-'~0.123 

J-Jil 

WRITE(LUN' J)U 
CALL CLOSE (3) 


Subroutine TSISBl is a completion routine called by the CALL RTS 
routine upon completion of the sampling. 


TEJISDJ. WITH TSIRUN.SAU NOU 7V 
SUDROUTINE TS.TSD1 
COMMON /DLKl /IDUF ( 5120 ) 

COMMON /DLK2/IDEr p ICMJ y IC| IAN J, 
lNCODE:.^IURrrU( 102-40P IDIIFp 1 p ICHANl ) 
IDEF-IDEFTJ, 

return 




TABLE 1 


BASIC DATA SLOT HEIGHT - .021" 


Cy ■ 0 

.0065 

6 - 

o 

in 

y«/2 “ 

0250" 

1 

9 


y(in) 

u/u» 



v/u« 



(P I,P») 
Nm 2 









-507 

Knail 

0.009 

1.703 

0.36 


-0.035 



-647 


0.012 

1.712 

0.48 


-0.011 





0.014 

1.674 

0.56 


0.012 



y*O.016 


0.017 

1.600 

0.68 


0.027 



-424 


0.020 

1.532 

0.80 


0.038 





0.022 

1.450 

0.88 


0.033 



35 


0.029 

1,259 

1.16 


0.023 



213 


0.036 

1.030 

1.44 

0.049 

-0.014 

0,023 

0.267 

173 


0.049 

0.737 

i.ge" 

0.036 

-0.100 

0.016 

0.072 

-124 


0.062 

0.709 

2.48 

0.037 

-0.122 

0.011 

-0.154 

-265 


0.076 

0.731 

3.04 

0.039 

-0.131 

0.010 

-0.182 

-286 


0.089 

0.751 

3.56 

0.038 

-0.133 

0.011 

-0.199 

-290 


0.102 

0.776 

4.08 

0.038 

-0.136 

0.011 

-0.166 

-273 


0.116 

0.803 

4.64 

0.037 

-0.134 

0.012 

-0.205 

-270 


0.129 

0.823 

5.16 

0.037 

-0.133 

0.012 

-0.159 

-268 


0.142 

0.845 

5.68 

0.038 

-0.135 

0.011 

-0.180 

-264 


0.209 

0.942 

8.36 

0.033 

-0.125 

0.012 

-0.165 

-234 


0.276 

1.005 

11.04 

0.027 

-0.118 

0.011 

-0.154 

-220 


0.342 

1.028 

13.68 

0 .ri 8 

-0. 105 

0.010 

-0.099 





cu - 0.0065 


y(in) U/U 



V 2 " O.onV U. - 31.9 


1.306 

1.290 

1.275 

1.237 

1.171 

1.665 

0.818 

0.692 

0.619 

0.623 

0.639 

0.671 

0.682 

0.704 

0.815 

0.901 

0.938 


0.013 

0.010 

0.010 

0,011 

0.011 

0.011 

0.012 

0.013 

0.012 



0.356 

0.261 

0.145 

-0.092 

-0.147 

-0.089 

-0.079 

-0.073 

-0.137 

-0.119 

-0.209 

-0.195 


y-0.016 

-371 

-275 

-159 

-90 

-8 

-62 

-127 

-163 

-178 

-192 

-185 

-177 

-165 

-131 










CU ■ 0.0065 


y{in) U/U^ |/^2/U, |u'v'(mV’ 



1.076 

1.103 

1.117 

1.117 



0.716 

0.590 

0.549 

0.551 

0.574 

0.586 

0.618 

0.718 

0.815 


0.256 

0.344 

0.'401 

0.487 

0.573 

0.630 

0.631 

1.032 


1.777 

2.178 

2.550 

2.923 

3.324 

3.696 

4.069 

5.989 

7.908 


0.041 


1.404 0.041 


0.043 

0.040 


0.045 

0.046 

0.047 

0.045 

0.041 


- 0.044 

- 0.027 

- 0.009 

- 0.003 

0.009 

0.009 

0.015 

0.020 

0.024 

0.007 

- 0.006 


0.044 - 0.018 

0.045 - 0.020 


- 0.014 

- 0.006 


0.037 

0.066 


0.018 


0.338 

0.322 

0.301 

0.033 


0.011 - 0.034 
0.012 - 0.096 
0.010 - 0.056 
0.010 - 0.072 
0.011 - 0.062 
0.013 - 0.186 
0.014 - 0.150 


0.015 - 0.250 














CW • 0,0065 


y«/2 - 0.0413” 


yUn) U/U^ y^yn/2 


0.792 


1.012 

1.026 

0.996 

0.973 

0.945 

0.852 

0.759 

0.624 

0.527 

0.469 

0.451 

0.444 

0.465 

0.574 

0.705 


0.218 


0.412 

0.484 

0.533 

0.702 

0.872 

1.166 

1.501 

1.840 

2.155 

2.470 

2.809 

3.123 

3.438 

5.061 

6.663 


0.037 

0.036 

0.037 

0.041 


0.037 

0.044 

0.046 

0.052 

0.053 

0.038 


- 0.002 

0.007 

0.012 

0.031 

0.045 

0.052 

0.051 

0.044 


0.039 0.032 


0.022 

0.012 

0.016 

0.023 

0.076 

0.116 


0.027 

0.021 

0.015 

0.015 

0.013 

0.014 

0.014 




0...93 

0.220 

0.096 

0.193 

0.131 

0.151 

-0.256 


-337 I WAZX 


y-0.016 


0.013 -0.308 













cw " 0.0139 


y«y3 “ 0«0274" U_ - 31.9 mi‘^ 


y(in) U/U. V/U. ^2/u, 


0.009 

0.012 

0.014 


2.351 

2.310 


0.328 

0.438 


0.014 2.232 0.511 
0.017 2.166 0.620 
0.020 2.103 0.730 
0.022 1.994 0.803 
0.029 1.724 1.058 
0.036 1.392 1.314 


0.020 

0.022 

0.029 

0.036 

0.049 

0.062 

0.076 

0.069 

0.102 

0,116 


-0.066 

-0.O31 


0.031 

0.053 

0.053 

0.044 


0.054 


0.027 

0.969 1.788 0.035 -0.116 0.014 

0.934 2.263 0.031 -0.147 0.011 -0.103 

0.959 2.774 0,030 -0.154 0.012 -0.102 



0.209 


0.987 I 3.248 0.031 -0.154 


1.031 4.234 


0.030 


1.129 7.628 


0.276 1.185 10.073 0.017 



0.008 



-1127 I wall 


0.011 -0.089 


1.009 3.723 0.028 -0.160 0.016 -0.055 


■0.157 0.011 I -0.101 


0.129 1,050 4.708 0.030 -0.160 0.011 -0.089 

0.142 1.066 5.182 0.028 -0.157 0.012 -0.071 


0.025 I -0.154 0.009 -0.089 


-0.055 


























TABU8 7 


cvi ■ 0 

.0139 

9 m 

35 ° 

V2 ■ 

y(ln) 

fl/o. 

1T/V2 


v/u. 

1 ^ 2 / U ^ 

0.009 


0.186 


- 0.044 


0.016 

1.621 

0.331 


- 0.013 


0.022 

1.611 

0.455 


0.003 


0.029 

1.561 

0.600 


0.009 


0.036 

1.483 

0.745 

0.054 

0.016 

0.031 

0.042 

1.389 

0.870 

0.052 

0.019 

0.031 

0.049 

1.770 

1.014 

0.052 

0.016 

0.030 

0.062 

1.088 

1.284 

0.053 

0.013 

0.029 

0.076 

0.897 

1.573 

0.046 

- 0.006 

0.023 

0.089 

0.796 

1.843 

0.044 

- 0.019 

0.019 

0.102 

0 . .46 

2.112 

0.046 

- 0.025 

0.015 

0.116 

0.746 

2.402 

0.048 

- 0.025 

0.013 

0.129 

0.799 


0.050 

- 0.022 

0.013 

0.142 

0.809 

2.940 

0.050 

- 0.013 

0.012 

0.209 

0.881 

4.327 

0.045 

0,034 

0,014 

0.276 

0.931 

5.714 

0.037 

0.069 

0.013 

0.342 

0.944 

7.081 

0.028 

0.100 

0.012 


OT.IgiaHT ^ .OiV 


U _ ■ 31.9 ms 


“ 1 


k»’v*(m*i"^ 


(P 


Jia 


p«) 

"2 


0.971 

0.871 

0.842 

0.790 

0.403 

0.191 

-0.037 

-0.094 

0.092 

-0.113 

-0.161 

-0.208 

- 0.220 


-847 


-743 

-709 

-572 

-385 

-223 

-108 

-35 

56 

57 
9 

-154 

-223 

-283 

-284 

-221 

-182 















CM * 0.0139 


■ O ‘ O504 " u - 31.9 


I 


y(in) U/U^ V/U^ 1^2/0, 



0.209 

0.276 


0.342 


1.376 

1 . 33S 

1.226 

1.097 

0.931 

0.784 

0.687 

0.633 

0.621 


0.787 

0.846 


0.878 


0.154 

0.274 

0.377 

0.497 

0.616 

0.839 

1.062 

1.301 

1.524 

1.747 

1.986 

2.209 


3.579 

4.726 


5.856 


- 0.053 

- 0.013 


0.016 


0.034 0.026 


0.041 

0.056 

0.053 

0.046 

0.040 

0.044 


0.050 

0.047 

0.044 

0.031 

0.019 

0.009 

0.006 


0.043 0.009 


0.048 

0.045 


0.063 

0.107 


0.141 


0.026 

0.029 

0.031 

0.027 

0.023 

0.018 

0.016 

0.013 

0.015 

0.014 


0.013 




1.135 
0.895 
0.441 
- 0.157 
- 0.244 
- 0.359 
- 0.108 
- O. 244 
- 0.261 


























T»1LB 10 




0.0367" 




BiM man m .o9i" 

).0367" 

0. • 31 

.9 m*“^ 



uW^(mV^ 

(P “ P.) 




HALL 1 



-2003 




y>0.016 




-1396 




-1053 




-672 


0.031 

0.740 

-493 


0.028 

0.684 

-367 


0.026 

0.713 

-346 


0.018 

0.373 

-540 


0.011 

-0.030 

-759 


0.011 

-0.048 

-801 


0.010 

-0.054 

-775 


0.010 

■0.036 

-762 


0.010 

■0.054 

-749 


0.011 

-0.057 

-729 


0.008 

-0.048 

-661 


0.008 

-0.085 

-560 


0.006 

-0.077 




















CM • 0.0197 


V2 " 0 . 04B0 " . 31.9 m .'^ 


y(in) U/0. y/y ^^2 V/U, 1/7^’Vo, 



0.014 

0.017 

0.020 

0.022 

0.029 

0.036 

0.042 

0.049 

0.062 

0.076 

0.089 

0.102 


1.853 

1.843 

1.821 

1.771 


0.292 

0.354 

0.417 

0.458 

0.604 


1.687 0.750 


1.571 


0.875 


0.056 

0.053 


1.320 


1.292 


1*113 1.583 


1.009 

0.984 


1.854 

2.125 


0.116 0.987 2.417 


1.000 


1.006 

1.053 


0.055 

0.047 


0.009 

0.022 

0.034 

0.047 

0.047 

0.053 

0.053 

0.050 

0.053 


1.445 1.021 1 0.053 0,053 


0.041 

0.009 


0.040 - 0.016 

0.040 - 0.022 

0.044 - 0.019 

- 0.016 
- 0.009 
0.028 
0.066 
0.091 



0.033 

0.031 

0.033 

0.032 

0.025 

0.018 


0.012 



0.595 

0.180 


0.014 - 0.043 
0.011 - 0.024 
0.011 - 0.041 


- 0.087 


0.012 - 0.145 
0.010 - 0.157 
0.010 - 0.160 


-1149 












TABtB 12 


CM ■ 0.0264 


W 2 * 0-0297" • 31.19 


y{in) 5/U^ 


0.009 3.5^8 


0.016 


0.022 2.934 

0.029 2.545 


0.036 

0.042 


0.076 

0.089 

0.102 


0.129 

0.142 

0.209 


2.053 

1.577 


0.049 I 1.326 
0.062 I 1.241 


1.260 

1.273 

1.288 


0.116 1.307 


1.313 

1.329 

1.354 


0.303 


3.348 0.529 


0.276 1.332 


0.751 

0.976 

1.202 

1.424 

1.650 

2.098 

2.5.49 

3.007 

3.434 

3.906 

4.343 

4.781 

7.037 

9.293 


0.092 

0.083 

0.055 

0.029 

0.026 

0.026 

0.025 

0.025 

0.024 

0.024 

0.021 

0.017 




-0.050 

-0.006 

-0.006 

-0.025 

-0.085 

-0.163 


-0.232 0.020 




1.012 

1.159 

0.755 


-2187 


-2542 


-0.288 0.011 -0.042 -1091 

-0.307 0.011 -0.080 -1100 

-0.317 0.009 -0.086 -1111 

-0.326 0.010 -0.077 -1110 

-0.332 0.010 -0.075 -1066 

-0.339 0.011 -0.072 -1087 

-0.342 0.011 -0.068 -1078 

-0.357 0.008 -0.058 -982 

-0.357 0.007 -0.062 -869 

















TABLE 13 


CM ■ 0.0264 

0 ■ 

25® 

y(in) 

u/u. 

y/ir./a 


0.009 

2.900 



0.012 

2.922 

0.299 


0.014 

2.893 

0.348 


0.017 

2.843 

0.423 


0.020 

2.787 

0.498 


0.022 

2.765 

0.547 


0.029 

2.608 

0.721 


0.036 

2.370 

0.896 

0.025 

O.U42 

2.066 

l.p45 

0.042 

0.049 

1.765 

1.219 

0.055 

0.062 

1.476 

1.542 

0.033 

0.076 

1.376 

1.891 


0.089 

1.373 

2.214 


0.102 

1.370 

2.537 


0.116 

1.382 

2.886 


0.129 

1.389 

3.209 


0.142 

1.389 

. 32 


0.209 

1.401 

5.199 


0.276 

1.367 

6.866 



y „/2 - 




U ■ 31.9 ms''^ 



-0.053 
-0.053 
-0.038 
-0.050 
-0.053 
-0.075 
-0.103 
-0.160 
-0.207 
-0.219 
-0.223 
-0.229 
-0.229 
-0.229 
-0.219 
-O. 207 


0.014 

0.020 

0.023 

0.016 



-2227 WALL 


-2440 

y-0.016 
-1405 

-960 

-730 

0.169 -648 

0.449 -566 

6 -513 

2 -863 

- 0.011 -1212 

-1223 
-1191 
-1189 
I -1163 
-1126 
-1000 
-869 















y„/a ■ 


y(in) u/u^ y/yn/2 V/u^ 


0.089 
0.102 
0.116 
0.129 
0. 142 
0.209 
0.276 


2.451 

2.464 

2.470 

2.549 

2.542 

2.508 

2.445 

2.332 

2.185 

2.009 

1.774 


1.386 

1.354 

1.348 

1.348 

1.357 

1.370 

1.348 


0.174 

0.232 

0.270 

0.328 

0.386 

0.425 

0.560 

0.695 

0.811 

0.946 

1.197 


1.718 

1.969 

2.239 

2.490 

2.741 

4.035 

5.328 


0.036 

0.048 


0.009 

0.019 

0.016 

0.025 

0.031 

0.034 

0.034 

0.034 

0.028 

0.031 

0.003 


0.036 -0.041 

-0.066 
-0.078 
-0.075 
-0.075 
-0.069 
-0.044 
-0.025 


0.023 

0.027 

0.021 



0.681 

0.028 

-0.024 


-1987 I HALL 


y«0.016 


-163" 


-1412 

-1321 


-1052 

-1039 

-1037 














TABLE 15 


CU >0.0284 


V2 ■ 


yUn) 

u/u. 



V/B. 


0.009 


0.029 

0.036 

0.042 

0.049 

0.062 

0.076 

0.089 

0.102 

0.116 

0.129 

0.142 

0.209 

0.276 


2.082 

2.194 

2.223 

2.207 

2.154 

2.065 

1.994 

1.824 

1.596 

1.426 

1.323 

1.273 

1.257 

1.254 

1.263 

1.235 


0.134 

0.236 

0.328 

0.432 

0.537 

0.626 

0.730 

0.924 

1.133 

1.326 

1.520 

1.729 

1.923 

2.116 

3.115 

4.113 


0.066 

0.063 

0.064 

0.065 

0.065 

0.059 

0.048 

0.037 

0.034 

0.033 

0.030 

0.023 


0 

>.019 
>.034 
0.028 
0.031 
0.031 
0.034 
0.034 
0.016 
- 0.003 
- 0.016 
- 0.025 
- 0.025 
- 0.022 
0.013 
0.053 



.038 
.039 
.038 
0.040 
0.034 

n a 


0.022 

0.016 

0.013 

0.012 

0.012 

0.011 


1.740 
1.676 
1.338 
0.986 
0.623 
0.177 
0.046 
‘ 0.013 
‘ 0.096 
-O.I J 


-1656 

-1470 

-1144 

-1046 

-642 

-699 

-621 

-562 

-446 

-462 

-594 

-662 

-730 

-746 

-638 

-524 



























1 


Cm * 0.0264 


6 -65' 



y(in) U/u^ |*^*/U, 


0.0 X 7 

0.020 

0.022 

0,029 

0.036 

0.049 

0.062 

0.076 

0.089 

0.102 

0.116 

0.129 

0.142 

0.209 

0.276 


II 


0.793 

0.843 

0.950 

1.078 

1.056 

1.169 

1.150 

1.191 

1.194 

1.232 

1.132 

1.066 

1.019 

0.824 

0.740 


0.052 

0.070 


0.687 0.081 


0.098 
0.116 
0.127 
0.168 
0.209 
0.284 
0.359 
0. 440 
0.516 
0.591 
0.672 
0.747 
0.823 
1.211 
1.599 


0.117 

0.109 

0.079 

0.092 


-0.075 

-0.053 

-0.053 

- 0.022 

-0.025 

-0.031 

0.006 

0.028 

0.082 

0.107 

0.163 

0.197 

0.179 

0.238 

0.260 

0.282 

0.194 

0.219 


0.121 

0.137 

0.123 

0.107 

0.193 

0.166 

0.154 

0.127 

0.089 


34.722 

47.790 

29.487 

14.851 

14.870 

0.468 

5.981 

5.874 

6.624 

5.203 

1.543 




























































FIGURE 1: EXAIPLES OF THE APPLICATION OF COANDA EFFECT 
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FIGURE 3: A SKETCH OF THE X-WING STOPPED ROTOR AIRCRAR 
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FIGURE 19: SLOT GEMCTRY MID SEAL1H6 
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FIGURE 20 i A CHECK ON m Tyn.ntMrycmuAi i 
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FIGURE 25i 
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B. MOVING AWAY FROM SURFACE - U^ - 31 ms 



A CHECK ON THE VIBRATION OF A HOT WIRE PROBE 
DURING A RADIAL TRAVERSE FROM THE SURFACE 
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FIG 26 

(A) & (B) 




(b) VlIM THROUGH HORXIMO SECTION WINDOW 


FIGURE 26 THE CIRCULATION CONTROLLED AEROFOIL LOCATED IR THE 
7ft X 5ft working SECTION 
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(o) GENERAL LAYOUT OP CONTROL CONSOLES 



(d) DETAIL OF PROBES MOUNTED ON AEROFOIL TRAILING EDGE (••• Also Pig. 26(a)) 


FIGURE 26 CONTINUED 
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(•) VIEW tiOOKXNO UPSTREAM AT ABROPOXL TRAXXiING BDQBi SMOMlNO BLOMXNQ 
SLOT AND TRAVERSING GEAR SUPPORTING THE PROBES 

FIGURE 26 CONTINUED 
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FIGURE 27s TUNNEL REFERENCE PRESSURE DIFFERENCE CALIBRATIOM 
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FI6URE 51: REYNOLDS SHEAR STRESSES IN THE HIGH SPEED WORKING SECTIOH BEFORE CLEAWlilS GAUZE 



FIGURE 52: EFFEa OF TIP JETS UPON OVERALL LIR COEFFICIENT 
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FIGURE 33A: SCHEHATIC DIAGRAH OF BLOWING SUPPLY 



UPPER AIR SUPPLIES 



figure 33B: BLOWING SUPPLY SYSTEM INSTALLED UNDER 7ft x Sft wn»yr.li; ^FfTinn 
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FI6UREW1 EXAUPLES OF THE EFFECTIVE INCIDENCE ESTIHATION 
TECHWIWE FOR THE UHBldWN » 














FI6URE H2-. VELOCITY PROFILES DOWNSTREAH OF CAI rRRATimi wn 77 lF EXIT 
AS MEASURED BY A DUAL SENSOR HOT WIRE PRHRF 


mji 




15 FROM SLOT 
C^-0.0158 


SLOT HT.-0,040in 


LZP-O.OlOin 


(VoltM) 


y(in) 

U(ra§“^) 

0 0.1373 

45.1 

0 0.1107 

43.4 

X 0.0840 

41.2 

a 0.0707 

42. a 

D 0.0573 

53." 

V 0.044 

.72.2 

+ 0.0307 

88.1 



3 4 

log (high pass frequency) 


ItlSMlSl 


ElOS 
























FIGURE 48i 









FIGURE 51: 
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FIGURE 52 : Bjj)CK DIAGRAM OF THE DUAL SENSOR HOT WIRE PROBE SAMPLING AND ANALYSIS PROG 
























NOWW. CALIBRATIOH of A SPI IT piim p.«ftF 












EJi.,&Q 



o 













IZE UPON RADIAL STATIC PRESSURE HEASUREHBIT 

































5,0 



FIGURE 68: VARIATIOHS OF CENTRELINE LIFT COEFFICIENT WITH GEOMETRIC INCIDENCE 
TRAILING EDGE BLOWING ONLY 
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FIGURE 72: DEDUCED EFFECTIVE INCIDENCE CORRECTION 
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FOR TRAII IMG EDGE BLOWING ONLY 



FIGURE 7H: EXMfflES OF THE LIFTING PERFORtWNCE OF OTHER AEROFOILS COWPARED WITH 
THE PRESENT STUDY 
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FIGURE 79Bs EXPERIMENTAL CENTRELINE PRESSURE DISTRIBUTIONS 

TRAILING EDGE BLOMIN& 




VARIATION OF LIFT COEFFICIENT WITH TRAILING EDGE BLOWING FOR A FIXED 
LEADING EDGE BLOWING RATE = n.ni25 
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figure 86: VELOCITY PROFILES IH THE TRA1LIR6 EDGE HAU JET C» = 0.0159 









FI6URE88. VElflCm PROFILES IN THE TRAILING ED6E WALL JET C* = 0.0197 
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FIGURE 92: LOCAL FLOW AUGLE IN THE TRAILING EDGE LV':.^ JET C. = 0.0065 
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LINES OF CONSTANT 
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FIGURE 104: WEAH STREAHL1HE PAnERM IN THE TRAILING EDGE WALL JET C. = 0.0284 





















FIGURE 108! CHANGES IN ANGLE FROM SLOT TO SEPARATION WITH 

TRAILING EDGE BLOWING FOR A VARIETY OF AEROFOILS 
AND EVALUATION TECHNIOUES 
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FIGURE U4: L0H61TUD1HAL TURBULENCE INTENSITY IN THE TRAILING EDGE WALL JET = 0.0139, ccg = 
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FIGURE 126 : COnPARlSON OF DEVELOPMENT OF THE LONGITUDINAL AND 
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FieURE 130: REYNOLDS SHEAR STRESS DISTRIBUTION IN THE TRAILING EDGE WALL JET = 0.0284. «g = 




































FIGURE 138: TRAILIHG EDGE STATIC PRESSURE DISTRIBUTIONS INDICATING SEPARATION POINT 
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FIGURE 142Bs RADIAL STATIC PRESSURE DISTRIBUTION IN THE TRAILING 
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FIGURE W7! SKETCH OF OBSERVED HOT WIRE OUTPUT SIGNALS 
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FIGURE 153s A SAMPLE OF THE VARIATION OF RADIAL PRESSURE 


GRADIENT IN THE TRAILING EDGE WALL JET FROM THE 
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VARIATIOH OF THE PROPOSED VORTEX STREAMING SPEED WITH ANGLE FROM THE SLOT 
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FIGURE 158: THE PROPOSED VORTEX SHEDDING FREQUENCY AS A FUNCTION QF THE VELOCITY 
DIFFEREHCE AT THE SLOT FOR A VARIETY OF TYPICAL LENGTH SCALES 









